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PREFACE. 



The present volume contains all -thftt -Professor Boole 
wrote for the purpose of enlarging his Treatis6 on Differential 
•Equations, Had he lived to publish the second edition he 
Tvould doubtless have incorporated his more recent investi- 
gations with the original work, and it is therefore necessary 
•to explain why another plan has been a;dopted. 
. In some cases Professor Boole had indicated that certain 
portions of the original work were to be omitted and their 
places supplied from the manuscripts; but on examination 
it appeared that in subsequent passages of the work there 
were references and allusions to the portions thus marked to 
he omitted which would not apply to the substituted matter. 
/Thus in attempting to carry out the directions it would 
have been necessary to accept the jesponsibility of making 
many alterations, and consequently to incur the risk of fail- 
ing in the attempt to improve the original form. 

Moreover the Treatise had been for some time out of 
print, and thie long delay which must have been caused by 
the labour of reconstruction would have ptoduced serious 
inconvenience to students at Cambridge and elsewhere. iPro^ 
fessor Boole himself was always especially anxious to consult 
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VI FKEFACE. 

the advantage of students, and those who had the charge 
of his manuscripts were naturally incEned to adopt a course 
of which they believed he would himself have approved. 

The design of reconstructing the Treatise was therefore 
abandoned; and it was resolved that the original volume 
should be reprinted, and that the manuscripts should be 
collected and published separately. This plan has the ob^ 
vious recommendation of enabling those who are already 
familiar with the original work to turn their attention 
readily to the new investigations. It will be seen that 
xnauy of the Chapters of the present volnme may be re- 
garded ad independent essays or memoirs which lose notbing 
by being separated from the other volume; and indeed so 
indications had been left by Professor Boole of the place 
which such Chapters were to occupy in the enlarged edition. 

I have printed ftU the unpublished matter relating to 
Differential Equations which I found among Professor Boole's 
papers* In a few oases it will be seen that an investigation 
is incomplete; such investigations have however been in» 
eluded in the volume, because I was unwilling that anything 
should be lost which so great a mathematician liad written 
on a subject he had long and carefully studied. 

I trust that no serious error will be found in the voltmie, 
and that any laults which may be detected will be excused 
on account of tiie nature and difficulty of the task that had 
to be performed. Many of the manuscripts had not been 
finally revised ; some of them were very obscure and had 
to be . carefully and laboriously copied for the press. In 
general the equations were not numbered, and thus tmly 
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PREFACE. Til 

blanks occurred in place of references; this circumstance 
often csnsed great trouble and perplexitj: I Hope however 
that a satisfactory result has been finally attained. 

I may state for the benefit of those who are conversant 
with the first edition of the original work Ihat the theo- 
rem which in tiie paresent volume is cited as contained in 
Chap. IL Art. 1 will be found in Chap. IV. Art. 2 of 
the first edition: ihe change was made by the direction of 
Professor Boole's interleaved copy. It was judged conve- 
nient to number the Chapters in the present volume in con- 
tinuation of those in the original work. 

All additions of my own are enclosed within square 
brackets. The sheets have been read by the Bev. J. Sephton, 
IFellow of St John^s College^ as well as by myself, and the 
volume is much indebted to his care and accuracy. Obvious 
mistakes in the manuscripts were of course corrected; thus, 
for example, the table at the end of the volume was calcu- 
lated by Mr Sephton, because the table in the manuscript 
was rendered erroneous by the use of a wrong sign in a 
formula* 



I. TODHUNTEE. 



St John's Coujeoi^ 
November, 1865. 
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CHAPTER XIX. 

ADDITIONS TO CHAPTER II. 

1. [In Chapter ii. Art. 9, two methods axe given for 
solving the differential equation 

{ax + br/ + c)dx+ {ax + Vy + c) dy = 0.] 

But there exists another transformation by which the equa- 
tion may be reduced to, (because it may be constructed from), 
an equation in which the variables are separated. 

Assume as this equation 

(^y+ C)daj'+ (^V+ CV/ = P (1) 

and let cc' = a; + w^y, y' = a? + m^. 

It will be seen that in these equations united we have as 
many constants as in the original equation. Now on substi- 
tuting in the assumed equation the values of oi and y\ and 
comparing with the equation given, we deduce a system of 
relations equivalent to the following, viz.: 

The quantities m^, m^ are roots of the quadratic 
am'-(J + a)w + J' = 0. 

The quantities -4, -4', (7, G' are determined by the system 
of equations 

Am^ + A'm^ = a\ Gm^ + G'm^ = c , 
B.D.E. II. 1 



Digitized by 



Google 



2 ADDITIONS TO CHAPTER II. [CH. XIX. 

from which we find 





m^ — m^ 


m^ — m^ 





Now (1) gives on dividing by (^V+ 0') (^i^'+ C) and 
integrating 

^ log [A'x' + C') + 2 ^^S (^y + ^) = ^^''^*" 

or (-4 V + (7') ^' (^y' + a)^ = const., 

which on substitution and reduction gives 

1 
{[am^ - a') {x + m,y) + cm, - c]<--c-' _ ^^^^^ ^^^ 

{{am^ — a') (a? + m^y) + cm^ — c'}«^V^' 

2. Under certain circumstances the ffeneral solutions of 
differential equations of the first order faiL This happens in 
the above example if m, = tw^, the solution then reducing to 

, 1 = const. 

The theory of the deduction of the true limiting form of 
the solution in such cases requires a distinct statement. 

Let the supposed general solution be represented by 

G being the arbitrary constant and u a function of x, y, and 
constants which are not arbitrary. Suppose too that when 
one of these constants h assumes a particular value ^, the 
function u reduces to a constant v. Then we have 

u-^v G — v 

Now the second member being a function of an arbitrary 
constant is equivalent to an arbitrary constant and may be 



Digitized by 



Google 



ART. 2.] ADDITIONS TO CHAPTER lU 3 

replaced by C. The first member is a vanishing fraction, the 

limiting value of which is (;t» ), the brackets being used to 

denote that after the differentiation h is to be made equal to k. 
Hence the solution becomes 

fdu' 



(S)='^- 



In applying this theory to the reduction of the general 
solution (2) in the case in which m^^m^, it must be 
observed that the numerator of the first member is the same 
function of m , x, y, as the denominator is of tw,, jc, y ; or 
attending solely to their functional character with respect to 
Wj, m^, we may affirm that the numerator is the same function 
of m^ as the denominator is of m^. Representing these func- 
tions by 0(>wj, <l>{ni^ respectively, we have 

But nij, wtjj being roots of a quadratic equation may be 
represented in the form 

m^ = m + k, m^ = m — k^ 

the roots becoming equal when k = 0. Hence 

^^ <l>{m + k) ^ 

Therefore since 

d<l>{m + k) _ d^{m + k) dj> (m — k) _ d<f> [m — k) 

dk " dm ^ dk "^ dm ^ 

we have 

,, T.d4>[m-^k) ,' ^ j.d^(m-k) 
du '^("'-^) dk -'^("^ + ^) dk 
dk~ ld>(m-A:)|' 

1—2 
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ADDITIONS TO CHAPTER II. . [CH. XIX. 

[dkj 



mm) ^d^Jm) 
J. 

d 



fdu\ ^ ' dm dm 



= 2^log^H. 



C 
Thus the solution becomes on putting (7 for — , 

^^ dm ^^T^' ^^S ^^"^^ "■ ^'^ (^ + ^^) + <^ - ^'1 = <^- 

3. [The next Article seems to have been intended to ap- 
pear in the enlarged form of Cliap. ii.; but I cannot discover 
what precise position it would have occupied. I conjecture 
that *' the above demonstration" refers to Chap. II. Arts. 2, 3; 
and I have accordingly supplied a reference to equation (3) of 
Chap. IT. 

I had myself drawn Professor Boole's attention to Chap. Ii. 
Arts. 2, 3. The geometrical process of Chap. II. Art. 3, ap- 
pears to have been first given by D'Alembert in his Opiw- 
cules, Vol. IV. p. 255. D'Alembert calls it a demonstration; it 
seems to me only an illustration^ at least in the brief form of 
the text : and that such was Cauchy's opinion may perhaps 
be inferred from the elaborate investigation given by Moigno, 
to which Professor Boole refers in Art. 5 of the present 
Chapter. 

I had also drawn Professor Boole's attention to the state- 
ment at the end of Chap. il. Art. 12, that only one arbitrary 
constant was involved. Accordingly Article 5 of the present 
Chapter developes this statement, and Article 4 seems intended 
to bear on the same subject.] 

4. In the above demonstration the relation between y and 
X is regarded as one of pure magnitude, and the interpreta- 
tion of the differential equation becomes a limiting case of 
that of the equation oi finite differences (Eq. (3), Chap. ii.). 
But if. we represent x and y by the rectangular co-ordinates 
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ART. 5.] ADDITIONS TO CHAPTER H. 5 

of a moving point on a plane the differential equation may be 
interpreted directly. For supposing it reduced to the form 

we see that the direction of motion is constantly assigned as 
a function of the co-ordinates of position. The entire motion 
is therefore determinate as soon as the initial point is fixed. 
The result of the motion is a line or curve wholly continuous 
or subject to irregularities according to the nature of the func- 
tion /(a;, y). That the arbitrariness of oriein is geometri- 
cally equivalent to the appearance of a single arbitrary con- 
stant in the relation connecting x and y may be shewn thus. 

Let y = ^(«o»yo»^) 

be the relation between x and y indicated by the supposed 
motion, x^, y^ being the initial point of departure.. Then this 
point being on the line of motion, x^^ y^ are particular values 
of X and y, so that we have from the above equation 

w^hich establishes a relation between x^ and y^, and shews 
that there exists virtually but one arbitrary constant. 

5. It is proved in Art. 3, Chap, ii., that the constants 
ip^, y^, initial values of the variables x^ y in the solution of 
the differential equation of the first order, are necessarily 
equivalent to one arbitrary constant. I shall shew from th« 
form of the above solution that this a priori condition is 
actually satisfied. 

Developing the expression for y [see Eq. (30) of Chap, ii.] 
in ascending powers of a?, we have 

y = J. + J.a,+^ + ^ + &c., (32) 

in which ^,=s- ^'i"v^°v~"°r > 

^ 1.2... (/I — rj 

the summation extending from n = r to n = oo. Forming 
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6 ADDITIONS TO CHAPTER II, [CH. XIX. 

hence the differential coefficients of A^ with respect to cr^ and 
y^, and reducing by (28), we shall find 

whence in particular 

Eliminate between these equations j^(a?^, y^^ and we hare 

dx^ dy^ dy^ dx^ 

Therefore, by Prop, i., A^ is a function of ^^, so that the 
solution reduced to the form (32) contains but the single 
arbitrary constant A^. 

It remains to notice that the solution must be applied 
only under the conditions of convergency, i.e. under the con- 
dition that the ratio of the n*^ to the {n — 1)*^ term tends to a 
limit less than unity as n tends to infinity. For a discus- 
sion of the failing cases of this test see ' Finite Differences,' 
Chap. V. Generally it is desirable, in order to secure rapid 
convergency, to divide the interval x — x^ into separate equal 
portions, to each of which the general theorem of solution 
may be applied. If a? — a?^ be very small the theorem may 
be approximately represented by 

3^-yo=/(a?e,ya)(^-a?«). 

On these principles Cauchy has founded remarkable methods 
of solution, which deserve attention from the commentary on 
the limits of error on their application by which they are 
accompanied (Moigno, Vol. IL pp. 385—434). 
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CHAPTER XX. 

ADDITIONS TO CHAPTER VIL 

1. [This Article relates to Art. 2 of Chap, vii.] 

The sense in which (9) may be said to constitute the 
general solution of the differential equation is this. We 
obtain from it 

giving any particular value to G this will geometrically 
represent a curve consisting of two branches, and giving to 
G every possible value we obtain an infinite system of such 
curves, each consisting of two branches. The aggregate of 
branches thus obtained is evidently the same as the aggre- 
gate of curves given by the two primitives (5) and (6), un- 
restricted by any connexion between c^ and c^. In this sense 
then the solution (9) is general, that it includes all the parti- 
cular relations between y and x which are deducible from 
the original primitives (5) and (6). And it is only in this 
sense not general that it groups these relations together in a 
particular manner. 

To the expression of the complete primitive a certain 
variety of form inay be given without affecting its generality 
in the sense above affirmed. Thus, if to the solutions of the 
component differential equations we give the forms 

ye^-c^ = 0, logy + aa? - 0^= 0, 

we should have, by the same procedure, as the expression of 
the complete primitive, 

{jfe^ - c) (logy + oa? - c) = 0, 
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8 ADDITIONS TO CHAPTER VII. [CH. XX. 

an equation which may equally with (9) be regarded as the 
complete primitive of the differential equation given, and 
which in geometry re{)re8ents the same totality of branches of 
curves as (9), with this difference only, that they are differ- 
ently paired together. 

2. [This Article relates to Art. 3 of Chap, vii.] 

The question will here naturally arise, Since if F= c be 
a solution of one of the component differential equations, 
f[V) = c,in which /(F) denotes any function of F, is also a 
solution, by Chap. iv. Art. 3, why not give to the complete 
primitive the form 

{fAVd-c]{MV;i^c} {/n(FJ-cl=0, 

or the stricter form 

fAK)Aiv;) /.(F.) = o {F), 

in which /(F^), ^(FJ, .../»( FJ denote arbitrary functions 
of Fj, Fj,..., F„ respectively — stricter because the presence of 
arbitrary constants and functions in the previous form is a 
superfluous generality? It is replied that though the form 
just given is analytically more general than (15), it is not 
more general than (15) with such freedom as is permitted 
in the interpretation of the arbitrary constants. In a physi- 
cal or geometrical application we should not only be per- 
mitted to assign a particular value to the arbitrary constant 
in (15), so deducing what in reference to its source would 
then be termed a particular primitive, but to combine the re- 
sults of different determinations of c together, so as to obtain 
every form of solution which is implied either in the func- 
tional equation (F), or in its component primitives 

Kj = Cj , Kj = Cj . . ., »^» = o„. 

The same considerations justify us in speating of (15) as 
the complete primitive, and not as a complete primitive. 
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CHAPTER XXI. 

ADDITIONS TO CHAPTER VIII. 



1. [The Singular Solutions of Differential Equations of 
the First Order received great attention from Professor Boole, 
and the Chapter devoted to that subject is one of the most 
valuable and important in his work. He continued his re- 
searches after the publication of his first edition, and intended 
to reconstruct the Chapter with great improvements in the 
second edition. After carefully examining the manuscripts I 
Qamc to the conclusion that it would be very diflScult to re- 
write this portion of the work so as to connect the old matter 
with the new ; and thus it seemed best to reprint the original 
Chapter VIII. with corrections of obvious misprints, and to 
print the matter intended for the revised form in the present 
volume. The plan gives rise to some repetition ; but this 
seems unimportant, compared with the advantage of preserv- 
ing in the author's own language all that he left on an in- 
teresting and important point which he had carefully studied. 

2. It may be of service to the student to reproduce the 
substance of some remarks on his Chapter viil. which were 
sent to Professor Boole soon after the publication of his first 
edition ; for there is evidence in his manuscripts that he paid 
great attention to such remarks while engaged in the revision 
of his work, and thus the reason and the meaning of some of 
his additions and changes may be made more obvious. These 
remarks will occupy the next Article. 



3 

are 



. The two pages beginning with " And these conditions 
sufficient," and ending with "do not lead to conflicting 
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10 ADDITIONS TO CHAPTER VIII. [CH. XX I. 

results^'* forming part of Arts. 3 and 4 of Chapter viil., seem 
obscure and difficult* The following may perhaps be substi- 
tuted with advantage. 

The only ways in which 

dy __ df{x, c) J % _ ^/(^> c) ^/(^> c) ^^ 
dx dx dx^ dx dc dx 

can be equivalent when c is variable, are 
(1) when^ = 0, 

. (2) when ^) =00; 

in the latter case ^ = qo, and therefore ^- = 0, and this 
dx dy 

implies that the singular solution is of the form x = constant. 

Thus there can be no singular solutions except such as 

df (x c) 
are found from ^^—^j- — - = 0, and such as are found from 
etc • 

X = constant. 

Similarly, if the complete primitive be expressed in the 
form x = F{y, c), there can be no singular solutions except 

such as are found from — -r^ — = 0, and such as are found 

dc 

from y = constant. 

In Art. 8 of Chapter viii. we read, " We may pass over 
the case in which the above equation is satisfied independ- 
ently of c, because the relation obtained would involve x 

only, while it is a condition accompanying the use of ^ = oo 

that it leads to solutions involving y at least." It is ob- 
jected, Why may we pass over this case? Such a case might 
occur and furnish a solution, and then we should want to 
know the character of that solution. Take for example 

p = x''y; here if w is negative, -^ is infinite when a? = 0, and 
this is a singular solution. For the general solution is y = ce" + 1 , 
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ART. 3.] ADDITIONS TO CHAPTER VIII. 11 

and 80 a?=0 is not a case of it. The words — while it is 
a condition.,, at least — seem very diflScult, for by supposition 

we are now investigating what is furnished hy -^^' 



■CO , 



Professor Boole met the objection in substance thus : 

" It will be found that the rules in the book are correct in 
this case. What is implied in the Chapter, though not stated 

with suflScient clearness, is that if -^ = oo leads to a solution 

dy 

which does not involve y in its expression, nothing is to be 
inferred whether it is singular or not. Then the proper test is 



dx \p) 

In this example we have 



-J--=zco gives a;" = 00 ; no inference ; 



^(i) = oo gives ar-'"*V = «'- 
Hence a; =• 0, provided n is between and — 1, or y = 0. 



Consider these separately : 

First. Let w be between and— 1, and a? = 0. This is 
by the test a singular solution. Substituting it in the com- 
plete primitive we get y = c, which confirms this. 

Second. Let y = 0. This satisfies the difierential equa- 
tion; but firom the fact that it comes from -3- f-j = Qo we 
have no inference ; from the fact that it does not come from 
-^ = GO we have the inference that it is a particular integral : it 
corresponds to c = 0. 
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12 ADDITIONS TO CHAPTER VIII. [CH. XXI. 

There remains the case of a: = when n is between — 1 
and — GO . As this does not satisfy -r- ( - j = oo , we infer that 
it is a particular integral. To prove this we have 

When a? = this gives, since 1 + n is negative, 

c = QO orc = — 00, 

according as y is positive or negative. This is like Ex. 2 of 
Chap. viii. Art. 8." 

The remark made by Professor Boole in the above replv, 

that «y -p = 00 leads to a solution which does not involve j 

nothing is to be inferred., As important. It corrects the state- 
ment put too strongly in Chap. viii. Art. 7, " All we can affirm 

is that if -^ = 00 gives a solution at all it will be a singular 

solution." 

From Art. 8 onwards it seems assumed that a solution for 

which -^ = is always to count as a singular solution, even if 

it should coincide with a particular integral. This does not 
seem to have been quite tne view of the former part of Chap- 
ter VIII. : see Arts. 5 and 6 of the Chapter. 

In Ex. 3 of Art. 9 we read, " the second is obviously a 
singular solution." This means that since we have a solu- 
tion which makes -^ infinite, we conclude that it is a singular 
dt/ ' 

solution. 

So in Ex. 5 of Art. 11 we read, " is evidently a singular 
solution," when it seems better to say, " and is therefore a 
singular solution." 

4. The additional matter relating to Chapter viil. begins 
with another example which was to be placed at the close of 
Alt. 3 of that Chapter.] 
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AKT. 4.] ADDITIONS TO CHAPTER VIIT. 13 

Ex. The differential equation 

has for its complete primitive 

Here ^ = -=£_-l ^^ ^ — 

_. rfy Voj* 4- V* — wi" die Vaj^ + v' — «»• 
Hence -j^ = p=^====:, -tt = ^^ . 

Both -T^ and -r vanish then if 
dc dc 

This therefore is the singular solution and it satisfies both 
the tests, as both x and y are contained in its expression. 

Of the partial tests 

^=0, ^ = 00, ^ = 00, 
dc dx dy 

the first is not satisfied, the last two are satisfied. 

The determination of c as a function of x by the solution 

of the equation -^ y ' = is equivalent to determining 

what particular primitive has contact with the envelope at 
that point of the latter which corresponds to a given value 
of X. 

One important remark yet remains. The elimination of c 
between a primitive y^'fi^t o) and the derived equation 

--^ = 0, does not necessarily lead to a singular solution in the 
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14 ADDITIONS TO CHAPTER VIII. [CH. XXI. 

sense above explained. For It is possible that the derived 
equation 

dc 

may neither on the one hand enable us to determine c as a 
function of .a?, so leading to a singular solution; nor, on the 
other hand, as an absolute constant, so leading to a particular 
primitive. Thus the particular primitive 

being given, the condition -^ gives 

e«=0, 

whence c is + oo if aj be negative, and — oo if a; be positive. 
It is a dependent constant. The resulting solution y = 
does not then represent an envelope of the curves of particu- 
lar primitives, nor strictly one of those curves. It represents 
a curve formed of branches from two of them. It is most 
fitly characterized as a particular primitive marked by a sin- 
gularity in the mode of its derivation from the complete pri- 
mitive. 

All the foregoing observations and conclusions may be 
extended to the case of solutions derived from the condition 

dc 

5. We have seen that the equation ;^ = may be satisfied 

by an absolutely constant value of c, so leading to a particu- 
lar primitive and not a singular solution. In this case 
>|r (a: + A, c) as well as -^ (a;, c) would vanish, and the nume- 
rator of (9), instead of being the difl^erence of a finite and an 
infinite quantity, would be the difference of two infinite and 
equal quantities. [See Chap. viii. Art. 8.] It would not there- 
fore be infinite. Hence we conclude that -^ would not become 

ay 

infinite for a particular primitive in the strict sense of that 
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ART. 5.] ADDITIONS TO CHAPTER VIII. 15 

term, i. e. for a solution derived from the complete primitive . 
by giving to c an absolutely constant value. 

This is one point of contrast between the conditions 

^^ = ^=00. 
do ^ dy 

There is another not less important. As the numerator 
of (9) may become infinite not only when i/r (a?, c) = 0, but 
also when '^ (a;, c) = infinite, we see that a relation between 

y end a? which makes -^ infinite will not necessarily satisfy 

the differential equation. On the other hand, it is not a par- 
ticular primitive in the strict sense of that term. 

Wo* 

Exactly in the same way the condition -r- = 0, as relating 
to the complete primitive, leads to the condition 

as relating to the differential equation, with the same points of 
difference in the respective applications. 

Ex. Let -^^zmy"^ , and suppose m a positive constant 
greater than 1. 

Here ^=(w-l).V^ 

which becomes infinite when ^ = 0. As jthis involves y and 
satisfies the differential equation it is a singular solution. 

To confirm this conclusion we may refer to the complete 
primitive 

which does not give y = for any particular value of c. 

Now let m be a positive constant less than 1. We have 
still ;^ = 00 when y = ; but this value oiy no longer satis- 
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16 ADDITIONS TO CHAPTER VIII. [CH. XXI. 

fies the differential equation. It is not a solution at all, nor 

would it result from the application of the condition ^- = 

to the complete primitive. The distinction of character of 
the two tests is here made manifest. 

6. We may express the most important results of the 
foregoing investigations in the following theorem. 

Theorem. Every solution of a differential equation of the 
first order which is derived from the complete primitive by 
giving to c a variable value will, if it involve y m its expres- 
sion, satisfy the condition 

dp 
and if it involve a?, the relation 



L 
dx \p, 



©- 



00. 



But relations satisfying these conditions will not neces- 
sarily be solutions of the differential equation. 

In applying this theorem the following points must be 
carefully attended to. 

1st. No conclusion can be drawn from the satisfying of 

the condition -^ = oo when the relation in question does not 

contain y in its expression, nor from the satisfying of 

d 



dx \pj 



QO 



when the relation in question does not involve x in its ex- 
pression. For these conditions being respectively derived 

ui/ dx 

from -f-^0 and ^=0 are subject to the same limitations 

in their application. 

2ndly. It may be that ^ or ^ (-j assumes for a particu- 
lar relation between x and y the indefinite form - , In this 
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ART. 6.] ADDITIONS TO CHAPTER VIII. 17 

case we must seek by the development of its terms or by 
other known modes its true limiting value or values. Finite 
values will indicate particular primitives, infinite values sin- 
gular solutions, and when such values emerge together out of 
the same relation between the variables, the solution will be 
a particular primitive possessing the geometrical properties of 
a singular solution, its locus will be a particular cuKve en- 
veloping other curves of the same family. 

See Examples 2 and 3 of Chap. vni. Art. lU 
We have seen that the conditions 
dp d f \\ 



dy ' dx 



0" 



indicate in general the existence of a relation between o and 
a; or c and y. And when that relation is such as to enable us 
to determine c as a continuous function of one of the vari- 
ables, the corresponding solution of the differential equation 
is singular, and is geometrically represented by an envelope 
of the curves of primitives. But it may be, as we have seen 
in a particular example, that the relation does not determine c 
as a function oixox y\ but according to the language already 
used, c is a dependent constant, or in some other way different 
from the constant of an ordinary particular primitive. Let 
us examine in particular instances the kind of singularity 
which may hence arise. 



Ex.l. Given p = ?^I^. 

^ X 



Here ^ = ^(^ + ^^8^)' 

This becomes infinite if a? = ; but this not involving y 
must be rejected. A^ain, it becomes infinite if y » 0, and 
this proves to be a solution of the differential equation, the 
limiting value of the indeterminate function in the second 
member being (Todhunter's Differential CalculitSy Chap, x,). 
Now the complete primitive is y = €**, discussed in Art. 4. 
The constant c is there shewn to be dependent, the solu- 

B. D.E. II. 2 
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18 ADDITIONS TO CHAPTER VIII. [CH. XX [. 

tion y = emerging from the complete primitive by making 
c = — 00 if ic be positive, and c = oo if a? be negative. 

Ex.2. Given (gy-aTyg+/logy=0. 

Here ^^xy±y{a?-^^offy)^^^ 

therefore / = — ^^^-^^ — ^^^ + 



dy 2 (ic*-4logy)i' 

and this is made infinite by y=0 and by a^ — 4logt/ = 0, 
i.e. by 

y = 0, y = €*. 

Both satisfy the differential equation. 

Now the complete primitive is 

y = e^'*^'. 

We see at once therefore that the second of the above solu- 
tions is singular. The first however is deducible from the 
complete primitive by making c = oo or c = — oo, irrespec- 
tively of the sign or value of a;, provided only that x he 
finite ; not so however if x be infinite. The value of c is not 
therefore in the most absolute sense independent of that of x. 
If from the complete primitive we seek the singular solution 

by the condition -^ = 0, we get the two equations 

6«^=0, aj-2c = 0. 

The second of these determines c as a function of a?, and 
leads to the second of the solutions obtained above. The first, 
though it does not determine o as a function of a?, still ex- 
presses a relation between c and a?, which is the ground of 
the fulfilment of the condition 

ay 
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ART. 6.] ADDITIONS TO CHAPTER VIII. 19 

We may further notice a peculiarity arising from this rela- 
tion. Supposing X finite and the solution y = a particular 
integral, it presents the singularity that it is the only case in 
which two particular integrals agree. We might in any com- 
plete primitive, by changing c into c', get two values of c for 
the same particular integral, but then it would be for evert/ 
particular integral. 

One negative character seems indeed to mark all the cases 
in which a solution involving y in its expression satisfies the 

condition -^ = oo . It is that such solutions do not emerge 

from the complete primitive by the attributing of a single and 

absolutely constant value to c. The relation which makes -^ 

ay 
infinite satisfies the difierential equation only because it satis- 
fies the condition ^ = 0, and this implies a connexion be- 
tween c and a?, which is the ground of a real though it may 
be unimportant singularity in the solution itself. 

At this point, then, the question arises, whether the term 
singular solution shall be confined to that class of solutions, 
the loci of which represent the envelopes of curves of primi- 
tives, or shall be extended to all solutions which, satisfying the 

condition -^ = qo, indicate the existence of a relation be- 
dy ' 

tween c and a;, and possess an actual singularity arising from 
this source. While the all but universal consent of mathe- 
maticians is in favour of the former course, it is to be remem- 
bered that the question is solely one of definition. Not such 
is the question how singular solutions of the envelope species, 
or as would more generally be said true singular solutions, 
are to be distinguished from all other solutions. This we 
now propose to consider. The question is not an isolated one. 
It stands in close relation to a series of properties of singular 
solutions which admit of an orderly development. 



2—2 

Digitized by VjOOQ IC 



20 ADDITIONS TO CHAPTER VIII. [CH. XXI. 



Discrimination of singular solutions of the envelope species. 

7. A negative test, which in the great majority of cases 
suflSces for the present object, is suggested by the following 
consideration. 

The differential equation determining -^ as a function of 

X and y determines also ^, ^, ad inf., and th& know- 
ledge of these enables us to construct in a developed form 
the complete primitive. See Chap. ii. Art. 12. 

dii d tj 
The values of ■-?- , --A , &c. ad inf, as derived from the 

differential equation, are the same as those derived from the 
complete primitive. 

But a solution deduced from the condition ^ = oo is only 

dti 
constructed so as to yield the same value of -^ as the given 

differential equation does. If it be of the envelope species, 
the curve it represents has in general no continuous contact 
with the curve of any particular primitive. It will not there- 
fore generally yield the same values for -^%, -y^, &c. as 

the differential equation does. It will not therefore generally 
satisfy the differential equations of an order higher than the 
first, which would be derived from the given equation by dif- 
ferentiation. Hence we have the following Proposition. 

Proposition. If a relation which makes -^ infinite satisfy 

ike given differential equation of the first order, hut do not 
satisfy all the higher differential eqv/itions obtained from it, 
such solution will be singular and of the envelope species. 
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Ex. 1. By comparison with Its complete primitive we saw 

in Art.. 5 that ^^my** has for a singular solution y = 

when m Is a constant greater than 1. 

We will first suppose m a fractional quantity greater than 
1, and endeavour to deduce the character of the solution with- 
out making use of the complete primitive. 

From the solution we have 

But from the differential equation 

g = (m-l)/-g = »n(m-l)y^ 
and generally 

^ = w(m-.l)...(w-r + %"'. 

Hence, when r is less than m, the substitution of y = gives 

as Iwfore. But if r is greater than m, it gives 

We conclude that the solution is of the envelope species. 
Secondly, suppose m a positive Integer greater than 1. 

In this case we find, when r is less than w, the same series 
of values as before ; but ioxr =m we have 

g = m(m-l)...l, 

and this also shews the solution to be of the envelope species. 
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Ex. 2. The differential equation 

is satisfied by 

aj» + y« = 4. 

Is this a singular solution or a particular integral ? 

From the solution we find 

dy _ X i^y _ 4 

From the differential equation we shall have 

substituting in which the value of -j- , obtained from the 
proposed solution, we find 

d^y (a^' + y'r _ 8 

Now this differing from the value before obtained, we con- 
clude that the solution is singular and of the envelope species. 

And this result is verified by comparing the solution with 
the complete primitive 

As the test above exemplified is merely negative, it is in- 
sufficient. For it is conceivable that an enveloping curve 
should have an infinite order of contact with each ot the curves 
which it envelopes, and this is also possible. Any test found- 
ed upon a comparison of the values of differential coefficients, 
any test therefore furnished by the Differential Calculus, would 
be insufficient for the discrimination of such cases. 

Ex. 3. Given -j-^y (log yY. 
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Here -f- = ^ gives y = 0, and this satisfies the diflferen- ' 
tial equation. 

From this solution we find 



y = *' ^ = 0, &c.arft«/. 



From the differential equation we have 

5=y{aogyr + 2(logy)'), 

which consists of y multiplied by a rational and entire func- 
tion of logy. It is easy to see that all the higher differential 
coefficients of y hence derived will possess the same character. 
And all such vanish with y. 

We can therefore neither affirm nor deny that the proposed 
solution is of the envelope species. 

8. Before demonstrating a general Eule for the discrimi- 
nation of solutions of this character, we shall notice certain of 
their properties which serve to indicate in what direction the 
Kule IS to be sought. [See Chap. viii. Art. 14.] 

As the exact differential equation differs from the sup- 
posed given differential equation by having acquired a factor 
which the singular solution makes infinite, so the given dif- 
ferential equation may be said to differ from the correspond- 
ing exact one by containing a factor which the singular solu- 
tion makes to vanish. If we knew that factor, we could by 
rejecting it reduce the given differential equation to a form in 
which it would no longer be satisfied by the singular solution. 
Now Poisson has shewn on a particular assumption, which 
does not however affect the principle of the demonstration, 
th^t this factor can be found when the singular solution is 
known. His demonstration is in substance as follows. 

Let us represent the given singular solution of the dif- 
ferential equation by 

tt = 0, 

u being a given function of x and y. Then introducing u and 
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X instead of y and x as variables, the differential equation 
after transformation will assume the form 

Now this equation being satisfied by t« = and the first 
member vanishing, the second must also. Poisson now 
assumes, and the assumption must be carefully noted, this 
second member to be capable of being developed in ascending 
positive powers of u. Supposing it so developed, the diffe- 
rential equation becomes 

g = ^i*- + i?^ + &c., 

in which A^ B,... are functions of a?, and a, JS,.., ascending 
positive indices. 

Hence if w = be a singular solution we have, putting p 

n du 
^''di^ 

But this demands that there shf)nld be at least one nega- 
tive power of u in the development in the second member. 
Therefore a— 1, the lowest index, must be negative. There- 
fore a being already positive must lie between and 1. 

We may give therefore to the transformed differential 
equation the form 

« being a positive fraction, and Q aot vanishing with u. 
Hence, dividing by w", 

-or -J- v}-^ = ft 

1 — a ^x 
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a differential equation which is not satisfied by u = 0, since 
u = gives 1*^"" = 0, and the first member vanishes while the 
second member does not vanish. In its present form then 
the equation is not satisfied by w = 0. We see also that the 
property of being satisfied by m = has been lost not in 
reality through a transformation, but through the rejection of 
an algebraic factor w* from the transformed equation. It has 
been shewn in the treatment of Clairaut's equation, how in 
the ascent by differentiation to an equation of a higher order 
a somewhat analogous effect is proauced, the singular solu- 
tion emerging out of a factor of that higher equation. 

If we inquire what is essential in Poisson's demonstration, 
we shall find it to consist in that the transformed equation is 
of the form 

£-«^. 

in which while Q neither vanishes nor becomes infinite when 
u=0, the functions 



Z7and 






both vanish with «. The question whether U is of the form 
tt* as Poisson supposes, or is not, is wholly immaterial. 
This will fully appear from the demonstration of the follow- 
ing theorem, which is in effect Poisson's freed firom arbitrary 
assumptions. 

S. Proposition. If u = be a solutian of a differential 
equation of the first order between y and a?, and 

dU j^ y. 

represent the form which that equation assumes when u and x 
are assumed as variables instead of r/ and x, then if f(x^ u) be 
resolved into two factors Q, Uy of which Q neither vanishes 
nor becomes infinite when i* = 0, while the functions U and 

-jj b@th vanish when w = 0, then the differential equation can 

be reduced to a form in which it shall <cease to be satisfied by 



i 
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In the statement of this proposition x is supposed to be 
constant in the integration relative to u. 

The differential equation after the transformation which 
introduces u and x as variables becomes 



Let 






SO that V is in general a ftinction of x and w, the form of 
which is known by integration when that of U is given. 
And again, transform the differential equation by making v 
and X the variables instead of u and x. We have 

/dv\ dv dv du 



(dv\ __ dv dv du 
dx) dx du dx* 



dv 
in which ^ is the differential coeflScient of v with respect to 

a;, on the above hypothesis as to the constitution of v as a 

function of x and m, while f;!^ ) is the differential coefficient 

on the hypothesis that v is reduced to a function of x alone 
by the conversion of u into a function of x, 

the above equation becomes 

(l)-S-«- 

Now if u = give r = for all values of x, it will there- 
fore give 

/dv\ 



©">■ 



and further, 

dv^_ d p ^_/v 

dx dxj^ If" ^ 
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since we are permitted to make w = before effecting the . 
differentiation with respect to x. Hence the equation re- 
duces to 

0=(2. 

And this is not satisfied, since by hypothesis Q does not 
vanish with u, 

-ff^O, the transformed differen- 
tial equation will no longer have w = for a solution. 

Cor. Assuming Q = l, which does not violate the hypo- 
thesis respecting Q, and gives 

U=f{x, u), 
we see that if 

be satisfied by w = 0, And if at the same -time « = gives 

du 



f 



o/(^»w) 



= 0, 



the differential equation can be transformed so as to cease to 
admit of the solution w = 0. 

It is obvious however that it is best to assume Q so as to 
make the subsequent integration for determining v the sim- 
plest possible. 

It is manifest that a salution which can thus be made to 
cease to satisfy the differential equation cannot be a particular 
primitive. For the complete primitive of the transformed 
differential equation which it does not satisfy is convertible 
into the complete primitive of the original differential equa- 
tion which it does satisfy, merely by writing therein for v its 
expression as a function of x and t/. It cannot therefore be a 
case of the compkte primitive m an^ sense. It must be 
a singular solution of the envelope species. 

The converse proposition still remains to be proved. 
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10. Proposition. If u = Ohe a singular solution of the 
envelope species of a differential equation of the first order, and 
if hy assuming u and x as the variables^ the differential equa- 
tion is reduced to the form 

du j.f . 

then will 

du 



JoJ 



become when u = 0. 
Let the complete primitive be represented by 

F{X, U) = Cy 

then, since 

dF(x, u) dF(x, u) du_ 
dx du dx'^ ^ 

we have if for brevity we represent J?" (a?, u) by F^ 

dF 
du __ dx ^ 

du 



therefore 



I 



du 



f[x, u) 




Now w = being a singular solution, F{x, 0) is not a con- 
stant ; for if it were, the complete primitive would, on giving 
to C the constant value in question, yield w = as a particu- 
lar primitive. And this would e(jually be the case whether 
that constant were finite or infinite in value. We see then 

that F{x, 0) must be a ftmction of a?, and therefore — \ * 

must either be a function of x, or a finite constant differing 
from ; the latter if F{xy 0) be of the form oa? + J, the former 
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if it be not of that form. Therefore the value of \-^ — - 

ax 

when u = 0, since in this we are permitted to make w = 

before differentiating with respect to a?, will be a function of 

Xy or a finite constant differing from 0. 

Now it is manifest that in general 



[■ 



1 dF, o-r<^^j 



dF du ^, 

where H is some value intennediate between the greatest and 
least values which -nv assumes within the limits of integra- 
tes; 
tion. When these limits are, as in the above case, infinitesi- 
mal, we have 

1 

dx 
„ C du \ CdF ^ 

dx Jo 

dx 

But we have seen that — \ ^ does not vanish. Hence 

dx 

its reciprocal, the first factor of the ri^ht-hand member of the 

above equation, does not become infinite. Again, 

F{x,u)''F{x,0) 

vanishing when w = 0, we have 

when u is made infinitesimal as was to be shewn. 
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It will be observed that the previous general expression 
for I 27 \ becomes infinite if w = is a particular integral. 

JET/ /\\ 

For then, F{xy 0) being a constant, — ^ ' ^ vanishes, while 

F{x, u) —Fix J 0) does not vanish so long as u differs bj 
however small a quantity from 0. 

These propositions form the ground of the following Rule 
for the discrimination of singular, solutions of the envelope 
species from all others. 

11. Rule. The proposed solution being represented by 
tt = 0, let the differential equation, transformed by making u 
and X the variables, be 

^+/(«,«)=o. 



dx 

Determine as a function of x and u the integral 

'du 






in which TJ is either equal to /(a?, w), or to /(a?, u) deprived 
of any factor which neither vanishes nor becomes infinite 
when w = 0. K that integral tend to with u the solution is 
singular. 

Ex. 1. Determine whether y = is a singular solution or 
particular integral of the differential equation 

Here, since w =s y, no preliminary transformation is needed. 
"We have I ^ "^ • 

which tends to with y. Hence the solution is singular. 

To verify this we observe that the complete primitive is 

1 
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and this cannot be reduced to y — O by giving any constant 
value to c. 

We have seen in Art. 7 that the test which is founded upon 
the comparison of differential coeflScients does not suffice to 
characterize the above solution. 

* 

Ex.2. The equation ^ = ?L^M ia satisfied by y = 0. 
Is this solution singular or particular? 

Here also no transformation is required. We have, reject- 
ing the fector - which neither vanishes nor becomes infinite 
when y = 0, 

^^= log logy -log logo 

and this being infinite, however small y may be, may properly 
be said to tend to infinity as y tends to 0. The solution is 
therefore particular. 

It will perhaps appear at first sight as if in the above ex- 
ample we ought to write 

when y is made equal to 0. But the course of the demonstra- 
tion shews that the value of the definite integral must be first 
obtained on the hypothesis that u (in this case replaced by y) 
is finite, and then the limiting value which its exp^sion 
approaches to, as u approaches to 0, be sought. And in this 
case, since for all finite values of u however small the integral 
is infinite, its limiting value is infinite. 

The complete primitive in the above case is 

y = €«^, 

and the nature of the solution y = has already been dis- 
cussed in Art. 4. 
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History of the Theory of Singular Solutions. 

12. It is remarkable that while the theory of enveloping 
curves and surfaces was at once founded and developea by 
Leibnitz in 1692 — 4*, the corresponding theory of the singular 
solutions of differential equations has been of very slow growth. 
The existence of these solutions was first recognised in 1715 
by Brook Taylor; it was scarcely more than recognised by 
Clairaut in 1734. Euler, in a special memoir, entitled Expo- 
sition de quehues Paradoxes dans le Calcul Integraly published 
in the Memoirs of the Academv of Berlin for 1756, first made 
them a direct object of investigation ; but the foundations of 
their true theory were only laid in 1768 in his Institution's 
Calculi Integralis. Laplace, Lagrange, Legendre, Poisson, 
Cauchy, and De Morgan have in various ways developed and 
extended that theory; but there has been so remarkable a 
want of unity and connexion in this long series of researches, 
that important portions of the theory appearing in a too 
isolated form have been neglected, forgotten, and rediscovered. 
I purpose here to give a brief account of what seems most cha- 
racteristic, rather than of what is most original in their several 
researches ; for the germs of nearly all subsequent discoveries 
on the subject are to be found in the great work of Euler. 

Taylor and Clairaut appear to have been led by accident to 
the noticing of singular solutions ; the former while directly 
occupied on the solution of differential equations, the latter 
while discussing a remarkable class of problems relating to 
the connecting properties of different branches of the same 
curve. Taylor gave them the name singular, while Clairaut, 
and Euler too in his memoir, regarded them as a species of 
paradox, not merely from their non-inclusion in the general 
integi'al, but from the mode of their discovery through a 
process of differentiation. The memoir of Euler, though it 
sheds no light on the real nature of these solutions, contains 

• Ada Eruditonmy 1692, p. 168 ; 1694, p. 311. Opera, Tom. iii. pp. 264, 
296. 

Methodus IncremerUorum, p. 26. 

M4m<nra de CAcadhnie dee Sdmeee, 1784, p. 209. 
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an interesting theorem concerning their connexion with the 
form of the differential equation, viz. If this equation can 
be brought to the form 

Vdz^Z{Pdx+Qdy), 

in which « is a function of x and y, and Zof z^ then will 

be a singular solution. Jn his Instituttones Calculi Inteffralts, 
Tom. I. p. 393, however, Euler gives a rule which is the 
counterpart of that of Cauchy. [See Chap. viii. Art. 12.] 
He shews .that if w = be a particular integral, and if the 
differential equation be reduced to the form 
du ,f . 



then I , ^^ ^ = QO , 

^0 ' 



du 

<f> (ic, u) 



The limits of integration are here supplied. The reasoning, 
which is not fuU^ developed, is the following. From the 
transformed equation we have 

^^ ^^ 
<tx = 



Hence 



X 



.o.j 



du 



^ (a;, u) ' 
X ^ : 1 f du 



G ^ ^ CV - 



If this be satisfied by a solution involving x and y, and if 
that solution be a particular integral, then on putting for x 
its value in terms of u and integrating, the above equation 
will be satisfied by giving some particular constant value to 
(7. But if the supposed particular integral be u = 0, then x 
and u being independent, we may perform the integration 
with respect to w as if a; were constant. The resulting equa- 
tion cannot be free from x unless C be infinite, and then it, 

B. D. E. II. 3 
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win evidently not be satisfied unless j—- r be infinite. 

J <f){x, u) 

We infer then that this is a necessary condition in order that 
w = may be»a particular integral. 

This is Euler's fundamental theorem, and from this, by 
means of an hypothesis agreeing with that of Polsson con- 
cerning the form of the transformed diflferential equation, he 
arrives at the condition 

dp 

[In the passage to which Professor Boole refers, Euler 
does not undertake to discuss the nature of any solution, 
but only of a solution of the form x = constant. On his 
page 408 Euler proceeds to discuss the nature of any solu- 
tion. Professor Boole seems to me to attribute too much 
to Euler. For the convenience of those who wish to ex- 
amine the original, I will give the reference to the passages 
in the later editions oH^vX^r' ^Institutiones CahuU Integralis : 
Vol. I. pages 343 and 355 of the -edition of 1792; Vol. i. 
pages 342 and 354 of the edition of 1824.] 

Laplace in the Memoirs of the French Academy for 1772, 
p. 343, established the tests 

dy~ ' dx \p) "" 

and shewed their respective uses. He established also the 
test which consists in the comparison of differential coefiicients, 
and he supposes it universal. He adopts the hypothesis of 
his predecessors as to the forms of expansion, but with some 
recognition of its insuflSciency. 

Lagrange in the Memoirs of the Academy of Berlin for 
1774, p. 197, and 1779, p. 121, appears first to have developed 
the theory of singular solutions in its two forms of derivation 
from the complete primitive and derivation from the differen- 
tial equation, and to have established the essential connexion 
of these. But supposing the differential equation to be ex- 
pressible in the rational form 
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and employing the differential coefficients of F{x^y^p) in- 
stead of those of p he was led to sacrifice rigour to symme- 
try. One of his results has often since been adopted as a 
test of singular solutions. It may be thus stated. 

Pbop. a singular solution makes the general value of 

^y . ' 

-^ , deduced from the difierential equation in its rational and 

d^ 


integral expression, to assume the form - . 

[The demonstration is given in Chap. vili. Art. 14.] 

• This ambiguity of value of -j^ is evidently but an expres- 
sion of the fact that the contact of a curve of the complete 
primitive and that of the singular solution is not in general 
of the second order. 

The result given in equation (5) of Chap. vill. Art. 14 has 
also been adopted as the test of singular solutions. 

The researches of Poisson and Cauchy have already been 
noticed. It is certainly remarkable that the final test to 
which Cauchy's analysis led should be essentially the same as 
that which had been discovered by Euler so long before. 

Professor De Morgan has thrown an important light upon 
the nature of the conditions 

dp ^ dp ^ 

dy^ ^ dx^ '* 

which are fulfilled by all singular solutions in the expression 

of which X and y are both involved. He has shewn that any 

relation between x and y which satisfies these conditions will 

d\ 
satisfy the difierential equation unless it make ^, as derived 

from the differential equation, infinite ; that it may satisfy the 

d^v . 
differential equation even if it make -t4 infinite ; lastly, that 

3—2 
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if it do not satisfy the differential equation, the curve it 
represents is a locus of points of infinite curvature, usually 
cusps, in the curves of complete primitives. 

The proof is as follows : 

Let p = il> (a?, y) 

be the differential equation. Then the -proposed conditions 
are 

^<^(^>y) ^^ # (a?> y) ^ CO 

dy ^ dx ' 

therefore by differentiation, 

J^.d^dy^^ ^ + .^^ = 
dxdy dji^ d£ ' dm^ dxdy dx ' 

whence we have 

d^4> d^4> 
dy _ dxdy _ dx* 

dx d^^ ~ d^<l> * 

^^- dxdy 

These are two equivalent expressions for the same value of 
-^. The question now is, under what circumstances this 

value of ^ will satisfy the differential equation. 
Now from that equation we have by differentiation 



whence 



d^y __ d<f) d^ dy 
dx^ dx dy dx^ 

d^y d(l> 
dy _ dx^ dx 
dx d<f> 

dy 
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If then -tK be finite we have, since -j? and -^ are both 
cur ax ay 

infinite^ 

dy _^ dx * 

dx d<l> 
dy 
and this by the rule for the evaluation of fractions of the 
form § is equivalent to the value in either of its forms before 

obtained for -f^. Hence, any relation which satisfies the 

d^v 
given conditions and makes ^ finite, will satisfy the difie- 

rential equation, 

d^v 
And the same result holds even if ^ be infinite, provided 

that ^ -^ -5? vanish. 
dx^ dy 

Lastly, as when this result does not hold, the failure is due 

d\ 
to the infinite value of -t4 , we see that the line in which the 

locus of the proposed relation intersects the curves of primi- 
tives will be a locus of their points of infinite curvature. 

[Transactions of the Cambridge Philosophical Society^ 
Vol. IX. Part II.] 

Legendre's Memoir of 1790 throws but little light upon 
the subject of this Chapter. But it exhibits the theory of 
the singular solutions of differential equations of the higher 
orders, both ordinaiy and partial, in a form of great beauty, 
and will be noticed in the proper places. 
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CHAPTEE XXIL 

ADDITIONS TO CHAPTEE IX. 



1. Br successive application of the second theorem of 
Chap. IX. Art. 13, a linear equation of the n^ order may be 
reduced to one of the (w — r)*** order, if r distinct integrals of 
what the given equation deprived of its second term- would 
be are known. 

The reduction may however be eflFected immediately by 
the method of the variation of parameters. In this and in 
most general investigations connected with differential equa- 
tions great advantages in point of brevity and of the power of 
expression are gained by the employment of the symbol of 
summation 2, and of the language of determinants. I shall 
exemplify this here. 

Suppose the given equation to be 

^~+^^^"^^^^- +A=X (1), 

and let y^, y^^.*-yr be r particular values of y, satisfying the 
equation 

^ + ^^d^+^«^ :^^- = ^ (^)- 

Thus y = c^y^ + cjj^>..'\-c^y^ 

is a solution of the latter equation including these particular 
solutions. We shall represent this by 

.y = Sc^i (3), 

and regarding the quantities c^, c^,.,.c^, represented here by 
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Ci as variable parameters, shall seek to determine them so that 
the above value of y may satisfy the equation given. 

These r parameters, enabling us to satisfy r — 1 arbitrary 
conditions, besides satisfying the diflferential equation, we may 
choose these so that 

dt/ cPi/ dr^y 

may be the same inform as if c^, Cg, . . . c^ were constant. Now 
from (3) 

whence 

Tx-^^'^'di' 

provided that the condition 

be satisfied. Differentiating the first of these equations, we 
find in the same way that 

Tx'^^'^'dx'' 
provided that the condition 

^ dyi dci ^ ^ 
dx dx 

be satisfied. And thus continuing we see that the system of 
r equations 

y-^'- 1=^1' P-^£^' «■ 

will hold true provided that the r — 1 conditions 
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te satisfied. In each of these equations the symbol S is to 
he interpreted by giving to i the successive values 1, 2,... r, 
and taking the sum of the results. 

Difierentiating the last of the equations (4), we have 

As we cannot impose the condition that the last term of 
this equation shall vanish, let z represent its unknown value, 
then 

g-^t*- <«>• 

Now the system of equations (5), together with 
^Q^ dci_ 

constitute a system of r simple algebraic equations deter- 
mining by solution the r quantities 

dc^ dc^ dCf 
dx^ dx '"' dx 

in terms of their coefficients and of z, and therefore in terms 
of X and z, since the coefficients are known as functions of 
X. It is evident also that as the second members of all the 
equations but one vanish, and the second member of that is 
z, the values so determined will be of the form 

Xj,Xj,... X^ being known functions of x. Thus the r un- 

dc dc 
known quantities j-^v j^ ^® made to depend upon only 

one omknown quantity, viz. z. It remains then to deter- 
mine z. 

For this purpose we must complete the expression of the 
differential coefficients of y, and substitute in the given dif- 
ferential equation, and then seek to satisfy that equation. 
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Now differentiating (6) we have 

on substituting for -j^ the value X^z as above determined. 

We observe that the coefficient of z is here a known function 
of X. If we differentiate this equation and in the result sub- 
stitute as above for -j^ , we shall have a result of the form 
ax 

L and M being known functions of x. Ultimately tten 
we have 

Thus, while y and the differential coefficients of y up to the 
(r — 1)"* are of the same form as if c,, c^,... c^ were constant, 
the succeeding ones differ in containing an additional portion 
consisting of z, and differential coefficients of z multiplied by 
known functions of x. The result of substitution of these 
values in the given differential equation will therefore consist 
also of two classes of terms, viz. terms under the sign of 
summation, which will be the same in form as if c^, c^^...c^ 
were constant, and terms involving the differential coefficients 
of z up to the (n — r)^y with multipliers which are known 
functions of x. We shall in fact have 

Sc (^^4. A ^-^ + A ^^ 4-a\ 
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Now yi being by hypothesis an integral of (2), the first 
line of the above equation vanishes, and there remains the 
linear equation of the {n — r)*** order 

^^+^^^' +^=^' 

Supposing z hence determined, we have in general 

Ci= IXizdx^ 
and hence 

y = yAX^z dx-V yJX^z dx +7/JX^dx, 

and as z will have n — r distinct values, each involving an ar- 
bitrary constant, the above equation will furnish n — r distinct 
values of y, each involving an arbitrary constant. It is to be 
observed that no arbitrary constant need be added in the inte- 
gration of the terms Xizdx, for the effect of such addition 
would only be to reproduce the known integrals c^?^^. In 
this way, however, the equation would represent the general 
integral of the differential equation given. 

2. Let us examine the form of the result in the particular 
case in which r = w — 1 . 

Here we have 

dx'^^^^'dx'^ 
from w = Oto77i = n— 2, then 



dz 
dx' 
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Accordingly the differential equation for z will be 



43 



h^i^^^'-A'-x.. 



dz 
dx 



■{>)■ 



Now the equations for determining 

dc^ dc^ dc^^ 
dx^ dx^'" dx 

become on putting XiZ for ^ , and writing for brevity t/\ for 
*,y,"forgs&c, 



dx 



dx 



da? 



yr^x,+yr'x,... + 2^,_r"^^, = 0. 

Whence, by the theory of determinants, 



y 1 dM r _ ^ ^^ 



^1 Jf dj^^t"-^) ' ^--^ "" iV/ d^^,^/"-'^^ 



Jlf standing for the determinant 



Viy y%^' 
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Now the determinant is ultimately a function of x ; and 
such indeed that 






i^&) 



For 



dM^ ^ dM dyi ^ dM dyl ^ dM dyr^ 
dx dyt dx dyi dx '" dy^'^^ dx 

^ dM , , _, dM „ , _j dM ,^1, ,-v 

Now M being homogeneous and of the first degree with 
respect to the quantities y^, y^, y^^, we have 

V dM ,^ 

Hence 2 -r- yl is what M becomes when in its expression 

Viy ys, —yH-i are changed into y^, y,', ... y«_„ therefore it is 
what M becomes when two of its rows of elements become 
identical ; therefore it vanishes. In like manner all the other 
sums in (8) vanish excepting the last, for yj^*"*^....^^^^*^*^ is 
not a row of elements of the determinant M. Thus we have 

Hence 

Thus the equation (7) becomes 

dz fl dM ^ .\ TT. 

therefore z = -^e"-^^!'^ flfe/^i*^ Xdx. 
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Hence, since 

<?(Cj -p- _ 1 AM 

, fl dM J 

whence ''*=}Md^''^' 

we have y = ty,j^ ^^ zda, ■ 

z being given above. 
In the case of X= 0, we have 



whence 
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CHAPTER XXIII. 

ADDITIONS TO CHAPTER X. 



1. The theory of singular solutions of differential equations 
of the higher orders has been presented in the most complete 
form which it has yet received by Legendre. {MSmoires de 
V Academie Royale des Sciences^ 1790, p. 218.) He determines 
first the possible forms of these solutions considered as emerg- 
ing from the complete primitive by the variations of its arbi- 
trary constants, and secondly the theory of their derivation 
from the differential equation itself. I shall follow the same 
order, and shall in the end endeavour to point out in what 
respect Legendre's theory may be regarded as complete, and 
in what respect it is imperfect. 

Suppose the differential equation to be of the n^ order, 
and let it when solved with respect to the highest differential 
coefficient of y be represented by 

y»=^(^,y,yi,y., -.-yn-i) (1), 

in which, for brevity, 

t,=^ v=^ v=^ 

Let also its complete primitive, solved with respect to y, 
be represented by 

y=f{x, a„ a„ ... aj (2), 

flj, a,, ... a» being the arbitrary constants of the solution. If 
we differentiate (2) with respect to x, regarding a^, a,, ... a^ 
no longer as constants but as functions of a?, so to be deter- 
mined as to leave the expressions for y^, y^, ... y^ as functions 
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of aj, a^, ... «„ the same as before, we shall have, on repre- 
senting the second member of (2) by/, 



whence 



__ df df da^ df da^ df da^ 

^^^ dx da^ dx da^ dx 5a,» dx ' 



y^ dx' 



provided that 



^f_da^,^^JS +_^^=0 
da^ dx da^ dx da^ dx 

Differentiating on the same hypothesis the first of these two 
equations, we find in the same way 






provided that 



^/_ da^ ^ d^f da, I f.f <?<^._o^ 



dxda^ dx dxda, dx 



dxda, dx 



And continuing thus, it results that the system 

_df _^ _^ 
^'~^' ^*~dx"-"^''~ (hr 



(3), 



will be satisfied, i.e., that 3^1, ^gj • • • Vn will have the same ex- 
pressions when flj, a2,...a„ are variable as they have when 
these are constant, provided that the law of their variation be 
determined by the conditions 



df da^ df da^ 
da^ dx da^ ^ 



^df da^^^\ 
da^ dx 



d^f da^ ^ d^f da, , d^f da^^^ 

da^ dx dx da^ dx dx da^ dx dx 

dj da^ ^ drf da^ , dj da^_^ 

da^ dx'^ ^ dx da^ da?*"* dx da^ dx*"'^ dx -^ 



(4). 
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In this system the coefficients of 

da? ' dx^ '" dx 

are known functions of a?, a^, a^, . • . a^ when the form b{ f is 
known. 

Eliminating 

da^ da^ da^ 

dx^ da? '*"<&? ' 

we have a relation between a?, aj, a,, . . . a^J and this relation, 
with the given complete primitive and the first n — 1 of the 
derived and reduced equations, viz., with 

y-h Vi-dx' ^»""da?'*"^'^*~dl»*-'' 

will enable us to eliminate a^, a,, . . . a^, and to obtain a rela- 
tion of the form 



'^ \^y y, 



dy dv QUq rs^ 



da?' db»'"*da:"-V 



This IS a differential equation of the [n — 1)*** order. It dif- 
fers in its origin from the given differential equation, in that 
a new relation between a?, a^, a,, . . • «„ has been employed in 
place of the n^ equation, derived by differentiation from the 
complete primitive, for the elimination of the constants. 

The differential equation of the {n — Vf^ order thus obtained 
has an integral expressing y in terms of ar, and w — 1 arbi- 
trary constants. This is the most general form of a singular 
solution of the differential equation. 

It is possible that the elimination of a^ja^, ...a» may 
lead to a resulting differential equation which, instead of 
being of the order ti — 1, is of the order w — 2, w — 3, &c. 
The complete integral of such equation would be a singular 
solution of the differential equation. These possible types of 
solutions are distinguished by Legendre according to the 
number of arbitrary constants which they contain. A soIut 
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tion containing w — 1 arbitrary constants is called by iim a 
singular solution of the first order ; one containing n — 2 ar- 
bitrary constants a singular solution of the second order, and 
so on. 

Adopting this language we might term the complete primi** 
tive a singular solution of the order 0. 

Lastly, any relation between x and y, which satisfies the 
given differential equation, will constitute a particular case, 
either -of the complete primitive or of one of the general 
forms of singular solutions above defined. In the case of 
differential equations of the first order it is seen that no arbi- 
trary constant caja appear in the expression of the singular 
solution* 

Ex. The equation 
has for its complete primitive 



y 



+ 5aj + a" + y. 



.(6), 



required its singular solution. 

Proceeding as above, we find on the hypothesis of a and h 
being variable parameters, the safne formal expressions for 

-^ , -T^ as if those parameters were constant, viz. 



•('), 



provided that the variation of a and h be such as to satisfy 
the conditions 



(f + 2a)|+(«+.i,|>0 



da 



das 



»0 



(8). 



B. D.E. II. 
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Eliminating hence ^ and ^ , we have 
a? 

And from this, the complete primitive, Ai the first of the 
derived equations (7) eliminating a and b, we find 

This is the difierential equation of the first order, by the 
•solution of which the most general form of the singular solu- 
tions of the given differential equation will be determined. 

Eeducing it to the form 

and integrating, we find 

This then is the general expression for the lingular solu- 
tions of the given deferential equation. We see that it in- 
volves in its expression one arbitrary constant. 

The differential equation (9) may properly be termed a sin- 
gular first- integral of the given differential equation. The 
singular first integral (9) has itself also a singular solution, 
viz. 

but this is not a solution of the original differential equation. 
Nor have we any right to expect that it should be so. A 
singular solution of a differential equation of the first order 
does not necessarily satisfy the differential equations of higher 
orders derived firoih that equation. Chapter xxil. Art. 7. 

•2. -Itremains to establish the theory of the derivation of 
the singular solution from the ; differential equation without 
the mediation of the complete primitive. 
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Resuming the differential equation in its reduced form (1), 
and representing its second member by ^, suppose an infini- 
tesimal variation given to the arbitrary constants of its com- 
plete primitive, and let the symbol S be used to denote the 
corresponding derived variations of y, y^, ...y^. Then we 
have 

and so on. Hence, substituting and transposing, 

Let us consider the real nature of this equation. 

If a value of y, suppose y = '^(a:), satisfy ihe given differ- 
ential equation, that vsilue substituted in the coefficients 

00 dip o 

of the above equation will convert them into functions of a:, 
and the equation itself will become a linear differential equa- 
tion, the solution of which will determine Su as a function of 
X. If the differential equation (10) be really, as it is appa- 
rently, of the w*^ degree, Sy will have n arbitrary constants, 
oCj, • . • o^, and will be of the form 

JF^, Pj, . • . P^ being functions of x. Hence 

We see thus 'that the given solution y^'^{x) will be a 
particular case of this general inte^al involving n constants. 
JLt will therefore be a particular mtegral of the proposed* 

4—2 
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If, owing to the constitution of its coefficients, the differential 
equation (10) be of the degree n — 1, we shall have 

and y = '^{x) will then be a particular case of a solution 
involving n — 1 arbitrary constants. It will therefore be a 
singular solution of the first order. Even so, if the differen- 
tial equation (10) be of the degree n — 2, y = -^(a?) will be a 
singular solution of the second order. And generally, if 
the differential equation be of the r*^ degree, y^-^^x) will be 
a singular solution of the order w — r. 

Besuming the equation (10) it is evident that it cannot 
be of the degree w — 1, unless -j^ ht infinite. For, dividing 

by -T-^, we have 
^ dyn.1 






in which the first term does not vanish unless ~- be infi- 

nite. This then is the necessary condition for a singular 
solution of the first order. For one of the second order we 
must have in like manner 



( 

and so on. 



rfyn-i dy^ 



It follows hence that to find the singular solutions of a 
differential equation of the n^ order, we ought to differentiate 

the equation, regarding y> ^> ^ > &c« ^ Tarying through 

the variation of the arbitrary constants, to form in this way 
a linear differential equation for Sy, to examine the conditions 
under which this equation reduces to the (n — I)***, or to a lower 
degree, and to examine whether the most general relation be- 



Digitized by 



Google 



asm 



ART. 3.] ADDITIONS TO CHAPTER X. 53 

tween x and y which satisfies such condition, satisfies also the 
given differential equation,. K so it may be regarded as a 
singular solution. 

Besuming the last Example, viz. 
and operating with 8 we have 

wbich rednces to a linear differential equation of the first 
order for determining hf, provided that we have 

Eliminating -^ from the given equation by means of this 
there results 

and we find on differentiating this that it does constitute a 
solution of the given equation. It is therefore a singular 
first integral of that equation. We see that it agrees with 
the result obtained under the same name in the previous 
Article, and the rest of the solution need not be repeated. 

3. Upon Legendre's theory, and upon its results, the fol- 
lowing observations may be made. 

1st. We learn from it that there may exist ^ different 
general forms of the solution of a differential equation of the 
fC^ order, viz. the complete primitive involving n arbitrary 



Digitized by 



Google 



54 ADDITIONS TO CHAPTER X, [CH. XXIII. 

constants, and general forms of singular solutions containing 
fewer than w.jarbitrary constants.. A solution y=:'\{r{x) of 
unknown origin being given, we construct a differential equa- 
tion for determining oy, and, solving it, form the expression 
for y 4- 8y, and from the .number of infinitesimal arbitrary 
constants it contains,- determine the nature of that general 
value of y of which the given value is a particulw case. 
Now we are not to infer from this that the form of y + Sy will 
be the same as the general value of y in question. But we 
may infer that it will be a form to which that general value 
is reducible. And the actual reduction will be effected by 
expressing the general solution (as is always possible) in a 
form permitting its expansion in ascending powers of the 
arbitrary constants, and in the expansion making these con- 
stants infinitOrSimal^ and rejecting all powers of them above 
the first. In fact, if 

be any general form of solution which, when we assign to 
a^, a^,,.,af particular Values (e.g. make them vanish) re- 
duces to 

y^^[x), 
theii we shall have 

the brackets denoting that after differentiation we make 
a^y «2» ••• ^r vanish. 

This is that limiting form of the solution which Legendre's 
method enables us to construct by the solution of a linear dif- 
ferential equation; and the ground of the safficlencjfr of his 
method consists in this, that the infinitesimal quantities 

oa^y Sas, ... Sa^y 

which are in fact the arbitrary constants of that solution, are 
equal in number to the arbitrary constants of the general 
unlimited solution, the nature of which is thus made known. 
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2ndly. Legendre's tests for differential equations of the 
higher orders are in kind and effect analogous to the tests 

dy ' dxp 

for differential equations of the first order, Thej enable us 
to decide whether a solution possesses singularity, not whether 
it possesses the envelope species of singularity. The comple- 
tion of Legendre's theory would consist in the discovery of 
those ftirther tests dependent upon integration which corre- 
spond to the test of £uler and Cauchy for difierential equa- 
tions of the first order. 
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[CH. XXLV. 



CHAPTEE XXIV. 



ADDITIONS TO CHAPTEE XIV, 



[Art. 1 was intended to follow Chap. xrv. Art. 2.] 

1. As the condition of dependence of functions of two 
variables is of fiindamental importance in connexion with the 
theory of ordinary differential equations, so the generalized 
condition of dependence of functions of any number of vari- 
ables forms a fundamental part of the theory of partial differ- 
ential equations. This is contained in the following proposi- 
tion. 

Prop. I. If Wj, u^,..*Un are functions of w^, ^j,*..^n> 
but are as such so related that some one of them is expressi- 
ble as a function of the others, or more generally that there 
exists among them some identical equation of the form 

F{u^, w„ ... wj = 0, (I), 

BO that as functions of aj^, a?,, ...a?^ they are not mutually 
independent, then, adopting the notation of determinants, the 
condition 



du^ du^ 



du. 



dx^' 


dx^ 


* dx^ 


du. 


du^ 


du^ 


dx,' 


dx^' dx^ 


du. 


du^ 


du^ 


dx,' 


dx^' 


*^n 



= 0, 



.(2), 



is identically satisfied. Conversely, if the above condition be 
identically satisfied, the functions w^, w^, ... w„ are not mutu- 
ally independent in the sense above explained. 
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First let^ it be noticed that the Proposition is hut a general- 
ization of that of Chap. ii. Supposing U and w to be two 
functions of x and y. the condition of their dependence is 
affirmed to be 

[ dU dU 
dx* dy 

du ^ 
dx* dy 

£• e. it is the result of eliminating da?, dy^ from the equations 

dV^ dUj ^ 
^dx + ^^dy^O, 



= 0, 



and therefore it is 



dy 
du J du J 

dU du dU du ^ 
dx dy dy dx" ' 



B» expressed in Chap* ii. 

We proceed to the genetal demonstration* 

Let the first member of (1), considered as a function of 
«j, «„ ... w^ be represented for brevity by F; then differen- 
tiating, we have 

dF^ dFj . dF J ■ 

from which it follows that if du^^ du^,... du^^ are equal to 0, 
then is du^ equal to ; or, since Wj> w,, ..« u^ are functions of 
Xij x^j ••• x^j that if 



du^ y . du^ J, . du. J 

du^ J , du^ J , du, J 



^dx, + ^dx. 



dx^ 



dx. 



+^*'.=o 



(8), 
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then is 

^i'^+t^-'+t^-" <"• ' 

Thus the last n equations, linear with respect to 

are not independent, and therefore by the theory of linear 
equations the determinant of the system vanishes identically. 
Now this is expressed by the condition (2). 

It remains to prove the converse, viz. that if the condition 
(2) be identically satisfied, the functions «j, w„ ... «^ will not 
be mutually independent. 

First, the n — 1 fiinctions w^, w,, ... w^ axe either mutually 
independent or not mutually independent. 

If not, then the n functions w^, w,, ... u^ are not mutually* 
independent, and the Proposition to oe proved is granted. 

If Wj, w,, ... w^j are mutually independent as functions of 
a?j, x^y ... x^, they may be made to take the place of w — 1 of. 
these quantities, e.g, o?^, a:,, ... x^^ in the expressing of n^, 
i. e. we may, by means of the expressions for Wj, w,, ... w^^, 
eliminate 6om that of w» the quantities oJj, 0?^, ... a?^^, and so 
express t«« as a function of Wj, u^,..,u^^ anda;^. Suppose^ 
this done, then the system (3), (4) will be converted into 

c?t*j = 0, dt^^Oy ^M^i = 0, 

-j-^^du. + -j-^du^.,.+'j—^du^. + j-^dx^r==0^ 
du^ * du^ „ . . du^^ *^* dx^ 

Now, the deterministnt (2) vanishing, the equations of the, 
linear system (3), (4) are not independent ; therefore those of 
the transformed system, as written above,, are not independ- 
ent ; therefore the last equation of that system must be a 
consequence of the others which manifestly are independent 
But from the form of that last equation we Bee that such can- 
not be the cade' unless we have 
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which implies that u^ is a function of w^, i«,, ••. u^^ merely. 
Hence the functions w^, i«„ ... ti^ are not mdependent, as was 
to be shewn. 

The first member of the equation of condition (2) is com- 
monly called the functional determinant of u^, t^,, ... i^« with 
respect to a?j, «,,..• a:^,. The proposition may therefore be 
expressed as follows. 

The condition of dependence or independence of any sys- 
tem of functions of as many variables is the vanishing or 
non-vanishing of the functional determinant of the system. 

On account of the ^at importance of this proposition it is 
desirable to illustrate it by an example. 

Ex. Are the functions 

a? + 2y + «, aj-2y + 3«, 2ary-a»+4ya- 2«* 
mutually independent or not? 

The equation of condition is 

1, 2, 1 

1, -2, 8 =0, 

2y — i8?, 2a?+4«, — aj + 4y — 4i5 
that is, 

-4(-a; + 4y-4«) + 8(2y-«)^2(2iB + 4«) = 0, 

which is identically satisfied. Hence the functions are de- 
pendent. In fact, representing them by Uj v^ to, we have 

4tr =s w* — «". 

[Art. 2 was intended to follow Chap. xiv. Art. 4.] 

2. As it has beto shewn that a primitive 

t. = *W (1) 

leads to a linear partial differential equation of the form 

i^ +.<2j=i? (2), 
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provided that w = a, t? = 5, are integrals of the system of ordi- 
nary differential equations 

dx_dy dz .„. 

p-^-Q^M ^^^' 

it is evident that we shall obtain a solution of the partial dif- 
ferential equation (2) by constructinff the system of ordinary 
differential equations (3), deducing their general integrals 

and then constructing from these the primitive (1). 

^ But the question arises, Will this be the most general solu- 
tion of the partial differential equation given? 

That it will be so, may be shewn by means of the general 
proposition. See Art. 1. 

For let ti? = represent any solution whatever of the given 
partial differential equation. Differentiating this with respect 
to X and y, we have 

dw dw^ dw dw _^^ 

dx dz^" ^ dy .dz^" ^ 

substituting the values of ^ and q formed from this in the 
given equation, we have 

which must be identically satisfied. 

In like manner, w = a,. « = J being solutions of the same 
equation, we find 

which must be identically, satisfied* 
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Eliminating P, Q, R from these three equations, it results 
that the functional determinant ofw^ u, t;, with respect to x, y, z, 
Tirill identically vanish. Hence w is a function of u and v, 
and the equation w = is a particular case of 

which is thus shewn to be the general integral of the given 
equation. 

We are thus led to the following general Rule. 

EuLE. To integrate the equation Pp-\- Qq=^R we must 
form the system of ordinary differential equations 

dx _dy _dz 

deduce their general integrals in the form 

u^a^ v^h^ 

and construct the equation 

F(u, v) = 0. 

This will be the general solution sought. 

[Art. 3 was intended to follow Chap. xiv. Art. 5.] 

3. The above theory maj be extended to linear partial dif- 
ferential equations of the first order, without regard to the 
number of the variables. ........ . 

First, the theory of the genesis of such equations is ex- 
pressed in the following proposition. 

Pbop» a primitive equatioix of the form 

F{u,,u,,...u:)^0...... (1), 

in which w^, Wj, . . . t«„ are any given functions of the vari- 
ables 0, dependent, and x., x^^ .., x^ indej)endent, will satisfy 
the linear partial diflferential equation obtained by eliminating 
dZf dx^j cw?2,...<&?^from 

rfwj=0, dfwa = 0, dw»«=0, 
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expressed as total differential equations with respect to the 
primitive variables, and the equation 

Of this important proposition I propose to give two distinct 
proofs. 

1st proof. Forming the total differential of the given 
equation we, have, on representing its first member hy -F, 

dF. dF ^ , dF ^ ^ 

Now this cannot be true for all forms of the function unless 
we have the separate conditions 

Au^^% d!M, = 0, cjIm^ = 0. 

Strictly to prove this, suppose i^^, i^„ . . . J?; to be any n 
distinct and independent ftmotions of ^>W-^**-w«, and as 
such, distinct and independent forms of j^ Then the above 
equation gives 

dF^^ kF.^ . cZF- 
dF^^ dF^j dF,j ^ 



dF^j ■ dF^j dF^ ^ 

Now j;, F^,:..F^ being independent, their ftlnctiobal 
detennmant with resr^ct to w^^m,, •..«„, does not vanish, 
ihis again is the condition necessary and sufficient that the 
above system of linear equations maybe independent; and 
this lastly being the case, their only possible solution will b^ 

^2*1 = 0, d!wj = 0, d!M^ = 0, . 

as was to be shewn. ': 
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These equations in their developed expression 
du. , . du. J , du. y , du. y 

^'^^ :....,+^&=o. 



^'l 



enable nb to detennine the ratios of dx^^ dx^, . . • db«, dz in 
the form 

dx^_^ _dx^_dz .. ^ 

:^ ~ js; - ~x~:b ^^^' 

where XjV-X,,.*.^., JS are functions of the, original vari- 
ables. And now, forming the equation 

p^dx^-{rjp^dx^... +p^dx^-dz = 0, 
and eliminating the' differentials, we find 

for the partial differential equation sought. 

2nd proof. Differentiating the given primitive with respect 
to ojj, as contained explicitly m tlie functions w, , Wj, . . . w^, and 
also implicitly in the same through z, we have, on represent- 
ing the first member of the equation by Fj 

dF(du^ du\-dF^(di^^ da\ 
duXdx^^'d»J^duXdx^'^^'dz)''\ ' 

dF (du^ du,\_ 

or 

dF du, dFdu^ ^d^^Ai^.dF^ ^q 
du^ dx^ du^ dx^ '" du^ dx^ dz ^^ ' 

since 

dF du^ . dF du^ , dF du„ dF 
du^ dz du^ dz '** du^ dz dz * 
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Differentiating thus with respect to the remaining inde- 
pendent variables, we obtain finally the system 



dFdu^dF du, 



dF du^ ' dF 



vi^^j dx^ "^ du^ dx^'" du^ dx^ dz -^* * 

dJF^du^ dFdu^ .3Fdu^ dF 
du^dx^ du^ dx^"* da^dx^ da*^" ' 



d£^da^ dJ^da^ .dFdu^dF 
du^dx^ du^ dx^"* du^dx^ dz^^ ' 

from which, in combination with the equation 

d^du^ dF^dv^ ,^^_^F^ 
du^ dg du^ dz *" du^ dz dz ^ 



we can eliininate 






dF 


dF 


dF dF 


du,' 


du,* du^' ds 


The result will be 








du, 


du^ 


du. 




tfo,' 


dxr" 


•• dx, '-P* 




k- 


du^ 
dx,"-' 


"dx,'^' 




%' 


^ 


dK 1 
•&■'-* 



or, converting rows into columns, 

du, du, du^ du^ 



dx,' 


dxr: 


'"dxj 


dz 


du„ 
dx,' 


du, 

dxr' 


du^ 
"'dx,> 


du, 
dz 


Pt>. 


?*>••' 


•"Pni 


-I 



= 0, 



= 0. 
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which is the determinant form of the result affirmed in the 
proposition. 

The second of the above forms of demonstration seems to 
be preferable to the first, in that it rests only upon the consi- 
deration of the one general form of the function X I have, 
however, given the two proofs, chiefly in order to illustrate 
an important remark, viz. that, in nearly all general re- 
searches connected with partial differential equations of the 
first order, two modes of procedure, the one involving the 
use of differentials, the other that of differential coefficients, 
may be employed, and that between the forms to which these 
respective modes give rise, a certain law of reciprocity will be 
found to exist. 

The theory of the solution of the partial differential equa- 
tion 

follows immediately from that of its genesis. If we repre- 
sent by 

the integrals of the system .of ordinary differential equations 
(2) a solution of the given partial differential equation will 
be represented by (1). That this will be also the most gene- 
ral solution may be shewn by the argument of Art. 1. For 
if II? = represent any solution, then since 

dw dw ^ dw , dw ^ 

we fiind 

^dw ^ ^dw ^ dw j.dw . 

from which, in combination with the corresponding equations, 



^dx^ ^dx^ .,. d^^ . dz ' 
B.D.E. II. 5 
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eliminating X, X„ ... X^, 5 we obtain a result which ex- 
presses that the functional determinant of t(?, w,,... Mj, with 
respect to the original variables is virtually 0. Whence w is 
a function of w^, Wj, ... w„, and the proposed solution is in- 
cluded in the one to which the above method of solution 
leads. 

That method may therefore be stated in the following Bule. 
EuLE. To integrate the linear partial differential equaiion 

form the system of ordinary differential equaiions 
dx^^dx^ _ dx^ _ dz 

and deduce their general integrals 

then, 

mil be the general integral s<yughU 



[The general observations were intended to follow Chap* 
XIV. Art. 6.] 

General observations, 

4. The relation which exists between a proposed linear 
partial differential equation and its auxiliary system of ordi- 
nary differential equations should be carefully studied. While 
it is pro})er to say as above that the general integral of the 
one requires the knowledge of all the integrals of the other, 
it is also proper to describe that general integral simply as 
the most general form under which an integral of the auxi- 
liary system can appear. If 

are integrals of that system, then 
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is the one general form of an inte^l of that sj^stem, and 
due regard being had to the arbitrarineas of F^ this is equi- 
valent to 

5. The form which the auxiliary system assumes wben 
the given partial differential equation is deficient in any of its 
terms should be noticed. 

If Xj = 0, the auxiliary equation 

becomes, on clearing of fractions, 

dx^ = 0. 
And thus, if ^, ^, ... ^ vanish, the given equation being 

"^^^•••■*'"^^«:^ •••■*■ ^-S^,""^' 
the auxiliary system will be 

da?j = 0, £foj=:0,...dav = 0, 

and the integrals of this system being of the form 

the general solution of the given equation will be 

F{x^, Xr, w,^i, wj =0. 

This conclusion would follow also from the principle laid 
down in Chap. xiv. Art. 2. 

Linear partial diflferential equations in which the absolute 
term is wanting, and which are therefore of the form 



^»4'^^*^."*"^^"^"^' 



5—2 
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may be termed liomogeneous. As In this case one of the 
auxiliary equations is 

the general integral will be 

F{u,,u^, ^«-i,«)=0, 

Wp Wj, w^i being found by the integration of the remain- 
ing auxiliary equations 

When Xj, X,, X^ do not contain z, the solution is best 

exhibited in the form 

z==<f>{u„u^, O- 

6. Every linear partial differential equation can be converted 
into a hom6genepus one containing one additional variable. 
For it is shewn in Art. 3, that if w = be any integral of 

then is 

^ du ^ ^ du du ^ ^du 

a homogeneous equation with a new variable. 

From the general integral of this equation, that of the 
former one may be deduced by making w = 0. 

7. The solution of partial differential equations is some- 
times facilitated by introducing a new system of independent 
variables. The actual transformation is greatly facilitated 
by the following symbolical theorem. 

Theorem. If the partial differential equation 
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be expressed, sjmbolically in the form 
in which 

then, if yi, y,, y» he a new system of independent vari- 
ables given in expression as functions of the old ones, the 
transformed equation will be 

For, regarding z as a function ott/^yt/^, y^, we have 

^ — ^ ^ _L — ^3 a, — ^y» 
dxi "" di/i dx^ dy, dx^ *** dy^ dx^ ' 



dz ^ dz dy^ dz_ dy^ dz dy^ ^ 

^n^dy^dx^ dy^dx^*'* dy^ dx^' 

whence, substituting in the given equation we find, as the 



total coefficient of -j- , the expression 

or symbolically, Ay^; and so on for the other coefficients. 
The result then is 

It remains only after calculation of Ayj, Ay^, Ay„, as 

functions of a?^, a?^, a?^, to express these fuactions and X 

intermsofyj,y„......y.. 

[It appears from the manuscript that an example was to 
have been supplied here.] 
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[The next Article may be considered supplementary to 
Chap. XIV. Art. 10.] 

Singular Bohtions of partial Differential Equations. 

8. Legendre*s theory developed in Chap, xxiil. for ordi- 
nary, may be applied also without essential change to partial, 
differential equations. Eegarding the independent variable 
z as receiving an infinitesimal change iz through infinitesimal 
change, not in the values of the independent variables 

but in the values of the arbitrary constants of the complete 
or in the forms of the arbitrary functions of the general mte- 
CTal, and performing upon the given equation the operation 
denoted by S, we shall obtain a linear partial differential 
equation for determining the general .value of iz corresponding 
to any particular given value of z. If that linear equation be 
of a lower order than the differential equation given, then the 
equation expressing the value oi z-\-iz will be a limiting 
form of a solution less complete or less general than the com- 
plete or general solution of the differential equation given, 
and the given solution, formed by making the infinitesimal 
constants in the limiting form actually 0, will be singular. 

Conversely, to deduce singular solutions without the know- 
ledge of the complete or the general integral, we ought to 
construct the equations of condition for the reduction of the 
equation determining hz to a lower order than the eauation 
given, and the most general solution of the differential equa- 
tions of condition so formed, will be the most general expres- 
sion for the singular solutions of the differential equation 
given. 

Ex. [px - jy)* J + 4»ia?' (« — ay) = 0, 

in which 

dz _ dz 
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Eepresenting the firat member of the equation by F^ we 
have, on operating by S, 

dp dxc dq dy cfe ^"" ' 

and the conditions 

dF ^ dF ^ 

necessary to reduce the equation for hz to a lower order give 

{px- qy) {px-^gj/) ^0. 

From these we find 

p = 3m*aj*y*, q = m^xy'^j 

definite and simultaneous values of j? and j, which being sub- 
stituted in the given equation lead to 

z = 2w*a?*y*, 

and this, as it gives the same values of p and q as those 
obtained before, will necessarily satisfy the given equation. 
It is therefore a solution, and &om the nature of the analysis, 
a singular one. 

Legendre shews that this singular solution is also dedu- 
cible from the general integral of the given partial differen- 
tial equation. That integral is the result of the elimination 
of a from the two equations 

{^ {a)Y — 2ax<l> (a) + a^r — mxy = 0, 

{^ (a) - ax] if> (a) - 2xif> (a) + ;5 = 0. 

To deduce the singular solution he supposes ^ (a) to be not 
simply a function of a, but a function of a and of one or 
both of the independent variables. He expresses the varia- 
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tion p{ 6{a) derived from this new source by S, and operating 
on the nrst equation wiith B, finds 

therefore ' ^ (^) = <^* 

Substituting this in the equations of the general integral, and 
eliminating a, we find 

^s before. 

Legendre states his theory of the derivation of the singular 
solutions of partial differential equations from the equations 
themselves with great brevity, but still as a general theory. 
And there is nothing in the statement that carries with it any 
apparent restriction upon either the order or the degree of the 
equations given. Until however we are in possession of a 
perfect theory of the genesis of partial differential equations 
we shall not be entitled to say that Legendre's theory of 
their singular solutions is a perfect one; for until then we 
cannot even define, in a perfectly general way, the nature of 
the operation denoted by 6. 
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[The next three Chapters all relate to the subject of partial 
differential equations bt the first order. The manuscripts do 
not appear to have received their final revision from Professor 
Boole. It is certain that he intended the contents of Chapter 

XXV. to form a part of the new edition; and it is highly 
probable, although not certain, that the contents of Chapter 

XXVI. and Chapter xxvii. were' also to be included. 

The three Chapters are mainly derived from two memoirs 
by Professor Boole, published in the Philosophical Trans- 
actions. 

The first memoir is entitled On Simultaneous Differential 
Equations of the First Order in which the Number of the 
Variables exceeds by more than one the Number of the jEqua- 
tions: it occupies pages 437*. .454 of the Philosophical Trans- 
actions for 1862. 

The second memoir is entitled On the Differential Equa- 
tions of Dynamics, A seqttel to a Paper on Simultaneous 
Differential Equations: it occupies pages 485... 501 of tlie 
Philosophical Transactions for 1863. 

The first memoir was finished before Professor Boole had 
seen Jacobi's researches, which are cited at the beginning 
of Chapter xxvi; these researches indeed could only just 
have been published. In his second memoir Professor Boole 
describes Jacobi's methods, refers to his own already pub- 
lished, and points out the nature of the connexion between 
them.] 



Digitized by 



Google 



( 74 ) . [CH. XXV. 



CHAPTEE XXV. 

ON SYSTEMS OP SIMULTANEOUS LINEAR PARTIAL DIFFEREN- 
TIAL EQUATIONS OF THE FIRST ORDER, AND ON ASSO- 
CIATED SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS. 

1. The term stmultaneoua is here applied to a sjstem of 
partial differential equations, to signify that in that system 
there is bat one dependent yariable, the general expression 
of which, as a function of the independent variables satisfy- 
ing all the equations at once, is the object of search. All 
linear partial differential equations of the first order being re- 
ducible to the homogeneous form, we shall presuppose this 
reduction here. Under this form indeed the problem actually 
presents itself in Geometry, in the theory of partial differential 
equations of the second order, and in Theoretical Dynamics. 

We are sometimes led, in connexion with the same class 
of inquiries, to systems of ordinary differential equations 
marked by the peculiarity that the number of the variables 
exceeds by more than one the number of the equations. Such 
systems are intimately connected with the former — stand 
to them indeed in a similar relation to that which the 
Lagrangean auxiliary system bears to the single partial dif- 
ferential equation from which it arises. The theory which 
explains this connexion, and grounds upon it the method of 
solution of both systems will form the subject of the present 
Chapter. 

Connexion of the Systems. 

2. Prop. I. The solution of a system of simultaneous 
linear partial differential equations of the first order may be 
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made to depend upon that of a system of ordinary differential 
equations of the hrst order in which the number of the vari- 
aoles exceeds by more than one the number of the equations. 

The system of partial differential equations being reduced 
to the homogeneous form, Chap. xxiv. Art. 6, let n be the 

number of the equations, a?j, ar^, a?,^ the independent 

rariables, and P the dependent variable. 

Then from the n given equations determining 

dP dP dP 

dx^^ dx^^ dx^^ 

we obtain an equivalent system of equations which, by trans- 
position of its terms to one side, assumes the reduced form 



^1 "rfaj^i "^n+/ 



dx^ " dx, 






dP 

dx,^ 



+ A 



dP 
dP 



= 



' dx. 



= 



H+r 



dP, . dP , ^ dP 



---^^ 



.(1). 



Multiplying these equations by the arbitrary constants 

\> \> \ 

respectively, and adding the results, we have 

dP 



dP dP 



+ (KJ.^ + \A^ -^-KA^ 



dx^ 



dP 



+ M„+ \A„ + MJ ^- = 



(2). 
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a Single partial differential equation which, on account of the 

arbitrarihesa of \,,X,, \_, is equivalent to the system 

from which it was formed. ^ 

Of this equation the Lagrangean auxiliary system will be 
Xj \, ' \^ * 



\-^ii + X,-4,i + \^-4^ 

^n^ 



.(3), 



whence, eliminating X,,Xa, \„ we have the system of 

ordinary differential equations 






.(4), 



These equations being included in the previous system (3), 
any integrals 

tt = a, t? = i, tt? = c, &c, 

of them will be integrals of it. Therefore u, v, i£?,... will be 
values of P satisfying the partial differential equation (2). 
For they will be the only values which can satisfy it inde- 
pendently of \, \,, \^. Hence they will satisfy the 

equivalent system (1), and the general integral of that system 
will be 

F{yb, v, w, ...) =0 (5), 

the form of F being arbitrary. 

Thus the relation of the system (4) to the system (1) is the 
Bame as the relation of the auxiliary system of a single linear 
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partial diiTerential equation to that equation. And the ground 
of this relation is seen to be the same in both cases. The 
one form necessitates the other. 

3. Instead of employing the above mode of deducing 
the auxiliary system, we might employ the following which 
is practically more convenient. 

Since any value P which satisfies the partial difierential 
equations determines P=c as an integral ol the ordinary sys- 
tem, the latter must be consistent with (3?P=0 in its de- 
veloped form 

^y . ^. dP dP dP 

• Ehmmatmg ^, ^, ^^ 

by means of the n given equations (1), we have 

^^(^•+i-^n^i-A*»^* -Ai^.) 

•^-JZ—idx^-A^^dx^-A^dx^ -Aa^n) 



Whence, equating to the respective coeflScients of 
dP dP dP 

we have the system (4). 

* 
. In the same way we can pass from the system of ordi- 
-nary to that of partial differential equations. From the equa- 
tion dP=^ 0, in its developed form, we must eliminate a number 
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of differentials cJr,, eZr,, ••• eq\wl to that of ike given eqwi- 
iions, and then equate to tke coefficienta of the remaining 

differentials. 

4 Lastly, the formal connexion of tte two systems should 
be noticed. The partial differential equations being given in the 
reduced form (1), the ordinary system may be constructed as 

follows : For any differential coefficient, as -^ — , in any 

column after the first, write the corresponding differential 
dx^^^ subtract from this the sum of da?^, da?,, dx^^ mul- 
tiplied respectively by the descending coefficients of that 
column, and equate the result to 0. The system of equations 
thus successively formed will be the auxiliary system sought. 

The transition from the ordinary to the partial system may 
be effected by the same rule, substituting only differentials 
for differential coefficients, 

[It appears from the manuscript that an example was to 
have been supplied here.] 

Up to this point the theory of systems of partial differen- 
tial equations is in analogy with that of single equations. 
But here a difference arises. We do not know beforehand 
what number of integrals a system of ordinary differential 
equations, in which the number of variables exceeds by more 
than one the number of the equations, admits. 

The theorjr which removes this difficulty will be developed 
in the following sections. It will be shewn that a system of 
linear partial differential equations which admits of solution 
by the assigning to the dependent variable a value which 
satisfies all the equations in common, must either itself satisfy 
a certain condition, or be capable of being developed into a 
new but equivalent system which will satisfy that condition. 
It will be shewn that when that condition is satisfied, the 
auxiliary system of ordinary, is capable of expression as a 
system of exact differential equations determimng the inte- 
grals sought 
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It will be faund convenient to express by a single symbol 
the aggregate of the, operations to which the dependent vari- 
able is subject in the expression of a partial differential equa- 
tion. Thus the equation 

dz dz , dz - 

may be expressed in the form 

A^j = 
if we assume 

d d , d 



^ = dt^''dx-^ydy' 



Under this convention the following proposition is to be 
understood. 



5. Prop. II. If AP = 0, A'P=0 represent any two 
homogeneous linear partial differential equations of the first 
order,, then will 

(AA'-A'A)P=0 

also be a homogeneous linear partial differential equation of 
the first order, and it will be satisfied by all the common 
integrals of the equations from which it is derived 

First, the equation will be linear. For, let a?, y represent 
any two variables whatever, or the same variable repeated, 
out of the set a?i .... a?^, and let -4, B represent any functions 
of the variables aj^ .... ir». Then A may be represented by a 

series of terms of the form A -7- , and A' by a series of terms 

of the form B -7- . Hence (A A' — A' A) P can be expressed 
by a series of terms of the form 
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which, on effecting the differentiations, becomes 

dx dy dy dx^ 

tlie terms containing the second differential coefficients of P 
mutually destroying each other. Hence the equation 

(AA'-A'A)P=0 

will he a homogeneous linear partial differential equation of 
the first order. 

The constitution of the coefficients of this equation is easily 
determined. For suppose the given equations to be 

SO that 

^'^•k *^i.- ^■=^.i. +^4- 

then the equation 

(AA'-A'A)P=0 

may be written in the form 

--■(^^^^.s -<). 

and, since terras involving second differential coefficients of 
P will disappear, this becomes 



(^^^-^'^.)S+(^^«-^'^«)S 



'2 



+ (A5,-A'^jf^ = 0, 
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We see from this that the form of the result is the same as 
tf the ^ or A' from either equation operated only on the coeffi-- 
dents in the other equation. 

Secondly, the above equation will be satisfied by all the 
common integrals of the equations from which it is derived. 

For, let ^ ;= c be a common integral of 

AP=0 and AT=0, 
then 

A^ = 0, A'^ = 0. 

Performing on these the respective operations A' and A, 
operations which involve only differentiation together with 
algebraic processes, we have 

A'A^ = 0, AA'^ = 0, 

whence, by subtraction, 

AA'^-A'A^ = 0, 

or (AA'-A'A)<^ = 0, 

jfrom which it appears that j> is also an integral of the 
equation 

(AA'-A'A)P=0, 

as was to be shewn. 

6, Prop. III. If by the above processes of reduction and 
derivation we convert a system of partial differential equa- 
tions into a new system, such that if expressed in the form 

A,P=0, A^=0, A^P=0, 

the condition 

(AA-AA)^=o 

shall for each pair of equations be identically satisfied, then 
the system of ordinary differential equations corresponding 
to this new system will admit of reduction to the form of 
exact differential equations, the integration of which will 
enable us to construct the general value of P satisfying the 
system given. 

B. D. E. IT. 6 
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1st. Suppose the given system of n equations reduced to 
tte form (1), marked by the peculiarity that n of the differen- 
tial coefficients appear only in successive equations and with 
the coefficient unity. Then taking any two of those equa- 
tions (we select the first two), we have 

A _ ^ _L J ^ . A ^ 

^»-^/^"^ "^^"^' 

^«=^,+^"^ +'^'^' 

from the forms of which we see that the derived equation 

(AA-A.A.)P = 

cannot contain either 

dP ^dP 
dx^ dxc% ' 

It can only, as appears from Art. 5, contain the differential 
coefficients 

dP dP 



and must "be of the form 

^dx^, "^^ "d^^ 

It cannot therefore be an algebraic consequence of any of the 
equations of the system (1) from which it was derived. It is, 
unless by the vanishing oi B^y,...B^ it present itself as an 
identity, a new equation algebraically independent. Com- 
bining this with the former ones, we have a system of n + 1 
equations admitting of the same reduction as to form fol- 
lowed by the same subsequent process of derivation. And 
the result of each of these completed steps is to convert the 
system into one containing one equation more than before; 
but containing in each of its equations one term fewer than 
before. The process must then end either in the genesis of a 
system of partial differential equations such that the further 
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application of the process of derivation of Prop. II. shall only 
lead to identities, or in the emerging of the system 

dP ^ dP ^ dP ^ 

dx^ dx^ dx^ 

The latter supposition would imply that P is a constant 
The consequences of the former we proceed to examine. 

The final system of linear partial differential equations 
will be of the same type (1) as the original system, but will 
differ from that system in that n will be increased, and r 
diminished by the same amount. We shall therefore simply 
state the form (1), only under the condition 

(AA-AA)J'=0, 

and with the altered values of n and r. 

First, then, the common integrals of the new system will 
be the same as those of the original system. This is evident 
from Prop. II. 

Secondly. If we write 

m = n + r, 

the first equation of the system (1) will be 

and the auxiliary Lagrangean system of this will have m -- 1 
independent integrals 

among which the n — 1 known integrals (Chap. xxiv. Art. 5) 



6—2 
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are Included. And the general value of P satisfying the 
above first equation will be 

The assumption P=iX^ would not satisfy the said equation, 
for it would lead, on substitution, to 1 = 0. Hence we infer 

that while the functions w^, t^,, u^^ are independent with 

respect to each other, they are also independent with respect 
toa?j, so that the m functions u^, w^, , w^^, a;,, are mu- 
tually independent in the sense explained in Chap, xxiv. 

Let us now transform the equations of the system (1) after 

the first by introducing w^, u^, w^^, x^ as independent 

variables. Those equations being 

A,P^O, A,P=0, 

the result of the transformation will be (Chap. xxiv. Art. 7) 
/A ^^P . fK \^P . /A ^ ^P . /A \^P ^ 



But P = Xj^ being an integral of each of the equations of 
the system (1) except the first, as appears from their forms, 
we have 

\^i = 0, A„a;^ = 0, 

thus the last terms In the transformed system vanish. Further, 
the coefficients of the remaining terms reduce to functions of 

«^i, Wg, t^^_i merely. For, considering the coefficient A^w^, 

we have 

• (AA-AA)^^. = o, 

Digitized by VjOOQ IC 



ART. 6.] EQUATIONS OF THE FIRST ORDER. 85 

which, since AjUj = reduces to 

AAt^i = 0. 

Hence A^w must be a solution of A^Ps 0, and therefore a 

function of w,, w,, w^^. And so for the others. It 

restilts therefore that the transformed system is 

, . V dP . , . V dP , . . dF 



Wj, ttg, w^j being the a^twjH independent variables of the 

system. 

But the transformation having involved no loss of gene- 
rality, for a new system of m independent variables was 
simply substituted for an old one, the condition 

(AA-A.AOi'^O, 

satisfied before, will continue to be satisfied in the new sys- 
tem represented symbolically in the form 

A^=0, A.P=0, A,P=0. 

Any common integrals of this system will also be common 
integrals of the previous system. For as functions of 



thev will satisfy tTie first equation of that system, and they 
Hvill satisfy the other equations, because the present system is 
but a transformation of those. The converse is equally mani^ 
fest. 

Thus a system of n partial difierential equations contain- 
ing m independent variables and satisfying a certain condi-^ 
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tion, has in virtue of that condition been converted into a 
system of n — 1 equations between m — 1 independent vari- 
ables, and satisfying the same condition. This then is con- 
vertible into a similarly constituted system of n — 2 equations 
containing w — 2 independent variables, and so on till we 
arrive at a final single partial differential equation containing 
m — n + 1 independent variables. This equation has m — n, 
that is, r integrals, and these are the common integrals of the 
system (1). 

But the system of ordinary differential equations corre- 
sponding to (1) is in number r, and is satisfied by all the 
common integrals of that system. Hence these differential 
equations must admit of reduction to the exact form. 

7. We may deduce from the above investigation the fol- 
lowing Kule. 

To integrate a system, of simultaneous linear partial diffe- 
rential equations of the first order. 

BuLE. Beduce the equations to the homogeneous form 
(I), express the result symbolically by 

A,P = 0, A^=0, A^ = 0, 

and examine whether the condition 

(AA-AA)P=0 

is identically satisfied for every pair of equations of the sys- 
tem. If it be so, the equations of the auxiliary system, 
Prop. I., will be reducible to the exact form, and their inte- 
grals being 

U — Qy v = if w = c, 

the complete value of P will be F{u, t;, w, ...), the form of F 
being arbitrary. 

If the condition be not identically satisfied, its application 
will give rise to one or more new partial differential equa- 
tions. Combine any one of these with the previous reduced 
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system, and again reduce in the same way. With the new 
reduced system proceed as before, and continue this methoci 
of reduction and derivation until either a system of partial 
differential equations arises between every two of which the 
above condition is identically satisfied, or, whidh is the only 
possible alternative, the system 

appears. In the former case the system of ordinary equations 
corresponding to the final system of partial differential equa- 
tions will admit of reduction to the exact form, and the gene- 
ral value of F will emerge from their integrals as above. In 
the latter case the given system can only be satisfied by sup- 
posing P a constant. 

Ultimately then the determination of P depends on the 
solution of a system of ordinary differential equations reduci- 
ble to the exact form. This does not mean that each equation 
of the system is reducible to the exact form, but that the 
equations may be combined together so as to form an equal 
niunber of equivalent equations of the exact form. Generally 
when we know this combination to be possible it is easy 
to effect it, and best to endeavour to do so. We might how- 
ever employ the method of the variation of parameters as fol- 
lows. Supposing p the number of differential equations make 
all but ^ 4- 1 of the variables constant, integrate the reduced 
system, and then seek to satisfy the unreduced system by the 
same series of integrals with the arbitrary constants as new 
variables. The successive integrations and transformations 
of this method would amount to the same thin^ as those 
upon which the second part of the demonstration of Prop. iii. 
rests*. 

Lastly, given a system of ordinary differential equations 
containmg a superfluous number of variables without know- 
ing how many integrals they admit, we must, supposing 
P= c to be any integral, construct the corresponding system 

* It was thas indeed that the aathor wai first led to that theory. 
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of homogeneous partial differential equations satisfied by P, 
and apply to them the foregoing Eule. 

8. Ex. Required the integrals of the simultaneous par- 
tial differential equations 

j-p Jjy jj-p 

Representing these in the form AjP=0, AaP = 0, it will 
be found that the equation 

becomes, after rejecting an algebraic factor, 
dP ,dP ^ 

and the three equations prepared in the manner explained in 
the Rule will be found to be 

dP^ dP . 

— +aj^=0 
dt dz 

No other equations are derivable from these. We conclude 
that there is but one final integral. 

To obtain it, eliminate 

dP dP dP 
dx* dy ' dt 
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from the above sjBtem combined with 

dP, dP, dPj dP. . 

and equate to the coefficient of -7- in the result. We 

find 

<fe — (^ + Za?) dx — ydy — xdt = 0, 

the integral of which is 

An arbitrary function of the first member of this equa- 
tion is the general value of P. 

[It appears from the manuscript that another example was 
to have been added here.] 
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CHAPTER XXVI. 

HOMOGENEOUS SYSTEMS OF LINEAR PARTIAL DIFFERENTLAL 
EQUATIONS. 

1. The theory of homogeneous systems of linear partial 
differential equations in which when expressed in the sym- 
bolic form 

A,P=0, A,P=0, A,P = (1), 

the condition 

(AA-AA)i'=0 (2) 

is for all combinations represented by i and J satisfied in 
virtue of the constitution of the symbols Ai, Ay, forms the 
subject of important researches by Jacobi {Nova Meihodus... 
Crelle's Journal, Vol. LX. p. 1). The following are the most 
important of his results. 

1st. An integral of any one equation of the system beinff 
found, other integrals of the same system may be obtained 
without integration, by a process of derivation founded upon 
the condition (2). 

Let <f> be an integral of the first equation of the system. 
Then is the equation 

identically satisfied. 

Also the condition (2) being satisfied in virtue of the con- 
stitution of the symbols, we have 

(A,Ay-AA)<A = 0; 
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and in particular, making « = 1, and separating the terms,. 

which reduces by a prior equation to 

It appears from this that A,^, if it do not reduce to a con- 
stant, is an integral of the first equation Aj<^=0, and, if it 
prove to be not a mere function of ^i, a new integral. 

This process may be repeated upon the new integral with 
a similar alternation of results. It will be evident from this 
that if we confine our attention to the two equations 

A,P=0, A^=0, 

and suppose, as before, ^ to be an integral of the first, then 
will 

\<f>> \{\<f>), A.{A.(A,<^)},... 

or, as these may be expressed, 

Aa<^, A,'<^, A/<^, ... 

be also integrals of the first equation ; and this process of 
derivation may be continued until we arrive at an integral 
Aj^<f> which is not independent, but is expressible as a func- 
tion of prior integrals 

\<f>> \'<l> A,*^^. 

and, sooner or later, such a result must present itself, since 
the number of independent integrals is finite. 

It is further seen that the most general symbolic form of an 
integral derivable from the root integral <f> is 

A,«A3^ A^**^, 

Of fif /A, being positive integers. 

The above remarkable theorem was in some degree antici- 
pated by the researches of Poisson, 
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2ndly. Jacobi shews how by the aid of such derived in- 
tegrals of the first equation of the system a common integral 
of the first and second equation may be found, and how from 
this integral and its derived series a common integral of the 
first three equations of the system may be found, and so on, 
until a common integral of the entire system has been as it 
were built up out of previous integrals of less general appli- 
cation. 

Let 6, ^\ <f>\ <;^<'*-i> represent a series, of independent 

integrals of the equation AjP=0, of which ^ is the root in- 
tegral, and the rest are derived from it by successive applica- 
tions of the operation denoted by A^, so that 

^' = A.0 ^(^-l) = A/-l^; 

also let A/<^ be not a new integral but a function of 

^,f, ^^-'^• 

Now <^, (j), <^(^-i) being particular integrals of AjP=0, 

the function F{<f>^^\ <^(^-'i)) will also be an integral of 

the same equation irrespectively of its form. Let us inquire 
whether the form of the function can be so determined as 
to render it also an integral of the second equation AaP= 0. 

We have then to satisfy the equation 

A.i^(^,<A', <^<^-i))=0. 

By the principles of the DiflFerential Calculus this equation 
assumes the form 

A # ^F .A • / dF ... -. dF - 

But A^ = f, A^' = ^" ..A,^<.''-«=si^*'-i)j 

lastly, A,^(^-*> may bv hypothesis be expressed in the folrin 
/(^, ^', ^*'~^*). Thus the equation to be satisfied ia . 



Digitized by 



Google 



X 



ART. 1.] PARTIAL DIFFERENTIAL EQUATIONS. 93 

A.' ^^ A. A" ^^ , ■ Alc-l) ^^ 

+/(^.f. "^'^■'')5^=o. 

a linear partial differential equation of whicli the auxiliary 
system is 

Now the integration of this system may be made to depend 
upon that of an ordinary differential equation of the (ji — 1)*** 
degree between the two variables ^('*-^) and ^. 

For we have 



I y 



d^ <f> ' d^ <f> 

Differentiating the last equation with respect to ^, and attend- 
ing to the former ones, we shall be able to express ^^ in 

terms of the variables <^, <f>, ^^-"^K Proceeding with 

this in the same way and continuing the process we shall be 
able to express the series of differential coefficients 

~l4~' d^^ ' d<f>^'^^ 

in terms of <f>y <f>\ <^<^"^>. From these fju-l equations, 

eliminating ^', ^", ^<'*-^, we shall have a final equation 

between 

*^'^ ' "#"' rf^ ' 

ttat is, a 'differential equation of the (ji — 1)* order between <}> 
and ^<^i). 
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The complete integral of this equation will be of the fonn 

<A^"?=/(<^.Ci.<5«, C/x-i)- 

DiflFerentiating this /a — 2 times in succession with respect to 
A, and continually substituting for the differential coefficients 
of ^(^-1) their values as before assigned in terms of 

^,f, r-«, 

we shall have a svstem of /* — 1 equations connecting the 

above variables with the constants c^, Cj, c^-i. Finally, 

solving these equations with respect to the constants, we shall 
possess the integrals required in the form 

FA<l>,<f>', r-«) = <a. 

and each of these will be a common integral of the first two 
equations of the given system (1). 

[On the back of a page of the manuscript the following 
paragraph occurs, which seems to have been intended as a 
simplification of the preceding argument which begins with 
"The complete integral."] 

Suppose that a first integral of the equation can be found. 
Its form will be 



<*.*-'.^ ^)= 



c. 



Substitute in this for the differential coefficients of <^^'*-^> 
their values before assigned in terms of (f>, (f>\ <f>\...<f)^^^\ and 
we have an integral of the system (3), and therefore a com- 
mon integral of the first two equations of the system (1). 

[We now return to the place at which we inserted a para- 
graph.] 

Just in the same way Jacobi deduces a common integral of 
the first three equations of the system (1). For representing 
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any one of the first members of the above system by yjr^ and 
deriving thence the new independent integrals A^y^^ A*yjr,... 
he substitutes an arbitrary function of these for F in the 
equation 

A3P=0. 

It is evident that the solution of the partial differential 
equation so found will again be reducible to that of an 
ordinary differential equation between two variables. And 
so the process is carried on till all the equations are satis- 
fied. 

2. The above remarkable process was developed by Jacobi 
in connexion with the theory of non-linear partial differential 
equations of the first order. In that particular connexion it 
admits of certain reductions tending to diminish the order of 
the differential equations to be integrated. But these do not 
affect the. general principle of the method. It was in this 
special form that the theory of the solution of simultaneous 
linear partial differential equations originated. Jacobi does 
not consider the theory of equations in which the condition 
(2) is not satisfied ; but the language in which he refers to 
the condition shews that he had speculated upon the general 
problem — and it is difficult to conceive that he should have 
meditated upon it and not arrived at its complete solution. 

[The manuscript here gives the first two words of 
the passage from Jacobi's memoir which is quoted in the 
Philosophical Transactions for 1863, page 486.] 
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CHAPTER XXVn. 

OF NON-LINEAR PARTIAL DIFFERENTIAL EQUATIOfNS OF THE 
FIRST ORDER. 

1. In treating the present subject we shall first consider 
that class of non-linear partial differential equations of the 
first order which involves two independent variables, and 
then proceed to the general theory. The reason for this 
procedure is that the particular theory, though of course in- 
cluded in the general one, rests upon a somewhat simpler 
basis, and it was in fact developed by the labours of 
Lagrange and Charpit long before the general theory was 
known. The latter we owe to the independent researches 
of Cauchy and Jaoobi. 

[Here the manuscript refers to the matter contained in 
Chap. XIV. Arts. 7 to 12 inclusive; and then passes on to 
the general theory.] 



General Theory. 
2. Given an equation of the form 

Z—^{X^, 0?,, ... Xny «!, a^, ... an)j 

the number of arbitrary constants a^, a^, ... a„ involved being 
equal to the number of the independent variables x^, x^, ... x^^ 
we obtain by differentiation and elimination of the constants 
a partial differential equation of the first order. Of this the 
proposed equation is said to constitute a complete primitive. 
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The form of the above process which it seema best, as 
throwing light upon the inverse problem of deducing the 
complete primitive from the partial differential equation, to 
employ, is the following. Let the given primitive, solved 
with respect to one of the arbitrary constants a^, be presented 
in the form . . . 

f{x^, ... x^y z, a^, ... aj =a, (1). 

Differentiating with respect to each of the independent vari- 
ables we have a system of n equations of the forms 



/i (a?,, ... a?^, «,i?i, a», ... «») =0 ' 



j^ (a;,, ... x^, z,p^, a„ ... aO =0 ^ 



:(2). 



These n equations enable us first to eliminate the w — 1 

constants a,, a^, and so deduce the partial differential 

equation sought in the form 

F^{x^,... a?^, «,i?i,-..i?»)=0 ...(3); 

secondly to determine the n — 1 constants as functions of 
iCj, ... x^^ ^9JPii '" Pn ^^ ^^^ forms 



F^{X^, ... Xn,Z,J)„ ...^n)=a/ 



Ki^iy ... ir», «,i?i, ...pn)=a^^ 



w. 



As the system formed of these w — 1 equations, together 
with the previous one, is merely another form of the system 
(2) obtained by directly differentiating the primitive, it follows 
that if frofn these equations we deduce the values otp^, ...i?« 
as functions of a?^, ... a?„, a^, ... a^, and substitate them in the 
equation 

dz=2>jdx^+p,^dx^+ ...+p^dxn (5), 

they will render that equation integrable, and its integral 
will be the complete primitive (1), the constant d^ being re-, 
gained by integration. 

B.D.E. IT. 7 
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Examining the system (3), (4) we see that the first mem- 
bers of all the equations which it contains are functions of 
ojj, ... aJn> ^> i'l* •••i^«> while the second members are con- 
stants. The question then arises, What mutual connexion 
exists among these functions in virtue of which they yield 
values of Pi, •..i>«> which render the equation (5) inte- 
grable? 

The answer to this question must involve the entire theory 
of the solution of partial differential equations of the first 
order, so far as relates to the determination of a complete 
primitive. Given a partial differential equation of the form 
(3) it is evident that if we can construct a system of associated 
equations (4) possessing the character above described, the 
final value of z obtained by integration of (5) will both 
satisfy the given equation and contain the requisite number 
of arbitrary constants. It does not follow from this that 
it will be the only complete primitive, but it will be a 
complete primitive. 

3. The relation sought is expressed in the following 
Proposition : 

Proposition. If 

represent any two out of a system of n independent equations 
such that the values ofp^^ ••• ^n) thence determined would make 
the equation 

dz -p^dx^ + jpjdfej + . . . +i?H^H 

integrahle, then the first members of these equations hein^ 
represented for simplicity by F and ^, the condition 

^UdF ^ dF\d^ dFfd^^ d^S\ ^ 

the summation extending to all values of i^ from Iton inclusive^ 
will be satisfied identically. 
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Beciprocally, if the above condition he satisfied identically 
fir each binary combination of functions in the proposed system 
of equations^ and if these functions be independent^ then the 
values ofp^ ,... p^^ as functions ofx^ , . . . a?„ , ^f, which they yield, 
will make the equation 



dz -Pydx^ •^V%'^^% + ••• +?H^i% 



integrabU. 



It will be convenient to begin with the particnlar case in 
which the proposed equations do not explicitly contain z, the 
particular pair to be considered being represented by 



n: 



Differentiating with respect to a?„ and regarding ^j, ...pn as 
functions of the independent variables, we have 



dXf, dp^ dxi '" dp^ dxi ' 

rf^ , rf^ ^ d^ dpn^Q 

dXi dp^ dXi *'* dp^ dxi ' 



(6), 



to which we may give the form 

dF^_^ dFdpj •^ 
dXi ^ dp I dXi ' 

d^ _ _^ d^dpj 

dxi "" ' dpj dXi^ 



(7), 



the summation with respect to^ extending from y = 1 to j = n 
inclusive, 

From the first of equations (7) multiplied by j- subtract 

7—2 
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the second multiplied by -r- , and sum the result with respect 
to i from t = 1 to z = w inclusive. We have 

= -SS /^— — ^'-— — -^"^ (8) 

* ' \dpj' dpi dx^ dpi dpj dxj 



The expression under the double sign of summation in the 
second member vanishes when t=y; we may therefore re- 
strict the summation to unequal values of i and j. Now 
as for any particular combination of values, e. g. 2, 3, there 
would exist in the completed member both the terms cor- 
responding to t = 2, J = 3, and those corresponding to / = 2, 
V=3, it is evident that if we employ the symbol 2y to denote 
summation with respect to different combinations of i and^/, 
the second member of the last equation may be expressed in 
the form 

^ fdF J4> dpf dF d^ dp^ 
*^\dpj dpi dxi dpi dpj dXi 

dFd^^_^dFd^dpA 
dpi dpj dxj dpj dpi dxJ ' 

\\dpi dpj dpf dpj \dXj dxJ) 

so that the equation (8) becomes 

^ /dFd^ dF^\ 
* \(iXi dpi dpi dxJ 

^^tU——''——\{^^--^ f9). 

^{\dpi dpj dpi dpil\dXj dxJ) ^ ^* 

The number of terms of which the second member ex- 
pjresses the sum is thus — ^- — - , and it will be observed that 
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as to any particular term it makes no difference in what order 
the numerical values of i and j are assigned to these quan- 
tities; e.g. whether for the combination 2, 3 we make i = 2, 
y=3, or i = 3, j = 2; but we must confine ourselves to one 
order. 

Now when the equation 

is integrable in the manner here supposed, we have for all 
combinations of i andj, 

dx^ dx^ * 

All the terms in the second member of (9) therefore Ta^ish, 
and we have 



\dxi dp^ dpi dxj 



This is the direct form of the Proposition under the parti- 
culfur limitation supposed. 

. As jP, ^ represent, under the same limitation, any two of 
the first members of the n equations (3), (4), which determine 

Pi9'"Pnf there will exist — - equations like the above. 

It is usual to employ for brevity the notation 

^ (dFd^_dFd^\^^^^ 

^\dxi dpi dpi dxJ '■ •'* 

and this being done the above system of equations expresses 

71 (ft *~* 1 ) 

the — ^^ — -functions of the form \FiF^ as linear homogeneous 

functions of the — ^-^ — - quantities of the form ;^*""^^* 

It is hence that the vanishing of the latter series of quantities 
secures the vanishing of the former. 

The converse truth will therefore be established by shewing 
that the ^ quantities of the form -^ " ^ ar^j when 
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i^j, ^,,...jP^ are independent with respect to ^j, p,,...^^, 
expressible as linear homogeneons functions of the — - 

functions [i^<j;]. 

To avoid complexity of expression I shall establish this for 
the particular case of n == 3, and shall shew that the reasoning 
is general. 

The functions jP^, JF^, 2^, being independent with respect 
to p^iP%y Pz7 ^^ determinant 

dF\ dF\ dF^ 

dp^ ' dp^ ' dp^ 

dF^ dF^ dF^ 

^' ^' ^ 

dF, .dF, dF, 

d'Pi ' dPi ' dPs 

does not vanish. This determinant we shall denote by A, 

In (9) writing for ^and ^ first F and jP,, secondly JP, and 
JPj, thirdly F^ and F^^ we have on changing signs the system 

•^*J \djp, dp, dp, dpj\dx, dxj 

\.^xP^ - \^p^ ^^^ dp, dpj \dx, dxj'^'" J 

dF dW 

Multiply the first equation by j-* , the second by ^— ^ , the 

third by -j-^ and add. Then 
^ dp^ 

-g i^^J-f tvj -f K^j - A (|;-^j. 
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whence as A does not vanish we have, on dividing by it, the 
function ^ - j^ expressed as a linear homogeneous function 

of [i^^J, [Vj%d [i^,i^J. 

(IW flW fJW 
In like manner multiplying the equations by J-* , ^-3, -^ 

at>^ dp^ ap^ 



respectively, and dividing by A, we obtain -^ — -P^ as a 
similar linear homogeneous function, and lastly, multiplying 
V ^> ^S ^', and proceeding as before, we obtain 

^* — ^ as a similar linear homogeneous function. 

From all which it follows that when [FJF^, [-^r^J* K-^J 
vanish, then 

^_dp^ ^_dp^ M_^ 
dx^ dx^^ dx^ dx^^ dx^ dx^* 

will vanish also. 

The reasoning is general in its nature. If -F;, F^, ...F^ 
are independent with regard to^j,^^, ...p^, the determiuaut 



dF\ 
dp/ 



dF\ 
dp^ 



dK 
dp,- 



dK 
dp^ 



= A, 



00, 



does not vanish. This determinant is from its constitution 
as a determinant linear and homogeneous, not only with 
respect to any row or column of elements, but also with 
respect to the possible binary combinations which can be 
formed of two rows or columns, ternary out of three rows or 
columns, &c. provided that these combinations are themselves 
of the form ot determinants. In the language of the theorjr 
such combinations are called minor determinants. Hence if 
we construct the system of equations represented by (10), and 
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observe that the coefficients of any particular term of the form 

-P- — P- in the several equations form a system of such 

axj dXi ^ 

minors to the general determinant (11), it will be plain that 
the equations can by multiplication and addition be brought 
to a form in which the coefficient of that particular term will 
be A. At the same time the coefficients of all the other 

terms of the form -p — -^ will vanish. For a little atten- 
axj aXi 

iion will shew that they will be what the determinant A 
would become on making two of its columns or rows of ele- 
ments equal, and therefore will be identically equal to 0. 

Thus the Proposition is generally established for the case 
in which z does not explicitly appear in the functions 

F,,F,, F,. 

When z does appear in those functions the equations (6) 
will be replaced by 

dF ^.^^1. +^^=,0 
dXi ^^ dz dp^ dx^. dp* ^^i ' 

dXi -^' dz dpi dXi dp^ dx^ ' 

from which it is seen that the theorem above established will 
only need to be changed into the form employed in the state- 
ment of the general Proposition. 

As the above is one of the most important propositions in 
the entire theory of Differential Equations, it may be desire- 
able to illustrate it by examples. 

[There are no examples in the manuscript.] 

4. We resume the general theory. 

The integration of non-linear partial differential equations 
may be effected by two distinct methods, both resting upon 
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the ground of the above Proposition. The first of these 
methods, originally established by a different analysis from 
that which will here be employed, was discovered by Cauchy 
{Exercicea d^ Analyse) y and rediscovered by Jacobi {Crelle^s 
Journal). The second method, discovered by Jacobi at a 
later period, forms the subject of his posthumous memoir. 
Nova Methodvs. ..... 



Cauchy^a Method. 

We will, as before, begin with the case in which z does 
not appear explicitlv in the proposed partial differential equa- 
tion, which we shall represent in the form 

F,[x^, a?n. P^y i>0 = (1). 

We have seen that to find a complete primitive of the 
equation it is necessary and sufficient to construct a series of 
equations 

• • (2). 

F^[^X ^n, Pi, i>n)=«„l 

such that not only shall the conditions 

[^i^J = 0, [i^,i^J = (3), 

connecting the new functions i^, J^ with JP'j, be identi- 
cally satisfied, but also the series of conditions 

[i^«i^J = (4), 

F^ and Fj, representing any two of the new functions re- 
ferred to. 

The first of the above series of conditions amounts to 
this, that J^, F^ must be integrals of the partial differen- 
tial equation 

KF] = 0....:,... ;.(5). 



Digitized by 



Google 



106 NON-LIKEAR PARTIAL DIFPEBENTIAL [CH. XXVII. 

It is the peculiar aim of Cauchy's method to determine the 
integrals so as to cause the second series of conditions to be 
satisfied also. And it is shewn that this will be attained if 
the integrals of (5), which form the first members of (2), are 
such that the particular values which i?^, ....^^^ assume when 
qs^ is made to receive anj constant value, as 0, are difierential 

coefficients with respect to ajj, x^^ of any single function 

of those variables, the form of which may be arbitrarily 
assigned. 

The necessity of this condition is obvious. If the general 

values of ^j, p^ are difierential coefficients of a function z 

with respect to x^^ x^, then the particular forms which 

p^y j?^y- assume when x^ receives any constant value are 

simply differential coefficients with respect to a?-, x,^ of 

what z becomes under the same circumstances. To prove its 
sufficiency we must shew that when it is satisfied the condi- 
tions represented by (4) will be satisfied also. 

Since F^ and F^ are mtegrals of [F^P] = 0, 

[i^,i^J = 0, [F,F,] = (6). 

Also, since if in (1) and (2) we give to a;^ a particular con- 
stant value, as 0, and then in (2) regard j?„ as a function of 

^l> ••.••••^fi_l> J?i> Pn-l 

determined by (1), the system (2) will virtually contain only 

of which jpj, p^^ are differential coefficients of a single 

function with respect to aj^, a;^j, it follows from the pro- 
position of Art. 3, that any two functions F^ and F^ will 
satisfy mutually the condition 

2<=n-l /^a ^» _ ^^« ^>| ^ 

*~^ \clxi dpi dpi dxiJ * 
the differentiations having reference to 
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explicitly as they appear in F^ and jFi, and implicitly as 
involved in j?». Thus the developed form of the above equa- 
tion is 

^4=n-i [{iK . dK Ap^ (dF, dF, dp^\ 
^=' \\dXi "^ dp^ dxj \dp, "^ dp„ dp J 

\dpi dp^dpj\da:^ dp^ dxJ} * 
the forms of -^ and -^ being determined from (1). 

Performing the multiplications^ the above equations will be 
reduced to the form 

^"^ \dxi dpi dpi dXiJ 

dFa ^i^n^ifdpn dF,_d^ AF\ 
dp^ **^ \dxi dpi dpi dxJ 

_ dF^ ^i=rn-i (dp^ dF^^dp^ dF\ ^ ^ ^ .^. 
dp^ *^^ \dxi dpi dpi dxi) ""^ 

But from the form of the total differential of (1) we see that 

dF-i dFx 

dp^ dXi dp^ <S7 

^"■"■^' d^r 5/ 

dpn dp^ 

Hence 

2i=»-i fdpn dK_^ dll\ 
*^^ Vc^« dpi dpi dxJ 

\rf/?n/ ^''^ \da?4 dpi dpi dxJ' 
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Now since by Art. 3 

we have • 

"^^"^ \dxi dpi dpi dxj 

^_(df\dK_dF[dK\ 
\dx^ dp^ dp^ dxJ ' 

therefore 

^*^^ \dxi dpi dpi dxJ 

\dpj \dx^ dp^ dp^ dxJ ' 
In the same way 

-*'-^ \dx^ dpi dpi dxJ 

\dpj \dx^ dp^ dp^ dxj * 

The suhstitution of these values in (7) gives 

^j^n-ifdFq dB\^dF^ dlj\ , 
**^ \dxi dpi dpi dxJ 

fdj\y idF^ /dF\ dF,_ dt\ dFA 
\dpj \dp^ \dx^ dp^ dp^ ((xj ' 

^dF^rdF^dK^diF^d^i 
dpnXdx^ dp^ dp^ dxj] ' 
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or i4=»-i/'^. ^'_^.^A 

^dF.dF\_dKdF,^^ 

^^^\dXi dpi dpi dxi) ' 
which is precisely the equation 

We see, therefore that to solve the partial differential 
equation 

•^i(^» a?,»i?,> i^n) = 0, 

it is only necessary to construct the linear partial differential 
equation 



2i=n(^dP^dF\dP\^ 
*-^V^< dpi dpi dxJ ' 



and to obtain w — 1 independent integrals of this 



Fn{^,y x^,p,, p^^a^, 

such that if we determine from these conjoined with the given 

equation the values of p,,....jo^, then those of^^, p^. 

shall, when x^ is made constant, be the partial differential 

coeflScients of one and the same function of 0:^, x^^ with 

respect to these variables in succession. 

Now provided that we can find all the integrals of the 
above partial differential equation the particular determination 
required may be effected in the following manner. 

The Lagrangean auxiliary system consists of 2n — 1 ordi- 
nary differential equations 

^'■■■■■-^rk' "i <*'• 

dPi . dp^ dx^ dx^ 
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These admit of 2n — 1 integrals, one of which will be jP^ = Cj ; 
and this will agree with the given equation if we make c^ = 0. 
We have therefore, besides the particular integral jFi = 0, 
2n — 2 integrals of the form 

(9)v 

Now suppose it required to find a value of « as a function 
of ajj, ...a:„, which shall satisfy the given partial differential 
equation, and shall reduce when x^ = (or any numerical con- 
stant) to a particular given ftmction of a?^, ...a?n-i> which we 
will represent by -^/r (a?j, ... aj^^). Then on the assumption 
that x^ = 0, we have first the given equation 

« = f (a?,,...a:^J (10), 

secondly the derived equations 



_ dyfr{x^,...x„,,) 



"\ 






(11), 



Make in the 2n — 1 integrals a;„ = 0, and suppose at the 
same time a?^, ... a;^j, jp^y ...jp«.i to receive therem the same 
values as in the above derived equations. Then from the 
3n— 2 particular equations which we thus possess in the two 
systems united (particular because under the assumption that 
a?„=0), we can eliminate the 2n — 1 particular values of 
x^, ...a;^„ p...,p^^, and so obtain n — 1 equations among the 
constants. These express the conditions which are necessary 
and sufficient in order that the values of ^, , . . . ;p^. j thus derived 
from the integral equations may, when a;^ = 0, agree with the 
values assigned in (11). Accordingly if we substitute in these 
equations of condition for c^, . . .Ca^, the general values 0,, . . . <f>sn^ 
we shall obtain n — 1 equations between aj^, ... a?^, ^j, ...pn> 
which will at once be particular integrals of the system (8), 
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and will possess the property that the values ofp^ , . . . jOn-i whicli 
they in conjunction with jP, = give will when x„ = reduce 
to the values given in (11). Hence these values with that of 
j)^ derived from the same equation will make 

dz—p^dx^" —p^dx^—O 

an exact differential equation. In the integral of this it will 
only remain to determine the constant so as to make the value 
of z agree with that given in (10). All the conditions will 
then be satisfied. 

We may collect the results of the above investigation into 
the following Kule : 

To obtain an expression for as a function of the inde- 
pendent variables a?j, ...a;„, which shall satisfy the partial 
differential equation 

F{x^, ... a?„, Pj, ...j7j = 0, 

and shall when x^ is made equal to (or to any numerical 
value) reduce to a given function of 0;^, ... aj^j, which we will 
represent by '^ (ajj, ... a?,.i). 

BuLE. Construct the linear partial differential equation 

^^^\dxi dpi dpi dxj ' 

and forming its auxiliary Lagrangean system deduce its in- 
tegrals 

in addition to the known particular integral 2^= 0. 

Between the above integrals and the equations 

difr (a?!, . . . a?^i) df^s^^^^^^^^p^^) 

eliminate, after making a?« = 0, the quantities 
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Iq the resulting n — 1 equations replace 

and we shall have a system of equations which with F=0 
will determine values of^^, ...^„, which will render 

dz "Pidx^ — ... —pndx^ 

an exact differential. The integration of this will give the 
integral sought. 

In the case in which the given partial diflferential equation 
is of the form 

z being contained.explicitly, the linear equation to be solved is 
^.^^{fdF dF\dP dFfdP dP\) ^ 

and the argument by which it is shewn that the integrals of 
this to be employed in conjunction with i^= for the deter- 
mination of ;?j, ...^n need only be so conditioned as to make 
Pii '"Pn-i differential coefficients of one and the same function 
of a7j,...ir^j when x^=0 is in character the same as that 
already developed in the present Article. It is only necessary 

to substitute in its exposition j~+Pi~T' ^^^ 'J~ y *^^ ®^ 

for the other functions. 

' But as the auxiliary system 





dx. dx^ dz 








TF ■■' dF dF 


dp. 


! 




dp, dp^ 





(12) 




di'' it' — dV dF 

dx,-^P^dz dx^^'^dz 


irta 


illy inclades tlie equation 

dz —pidx^ — .. . —p^dx, = 0, 







the ultimate expression of the Eule will be as follows t 
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To obtain an expression for « as a function of ar^, ... x^ 
which shall satisfy the equation 

and shall when x„ is made e^ual to (or to any particular 
constant value) reduce to a given function '^ (i»i, ... x^^ of 
the independent variables x^, ... x^^. 

Rule. Let 

be the 2w — 1 integrals of the auxiliary system (12) which are 
additional to the particular integral F= 0. Make in these 2» 
equations a?,» = and forming the further equations 

^_ dyl^{x^, ... x^,^) 
^'' dx. 



dyjr{x^, ... a;^0 



eliminate the 2n quantities a?j, ... x„^i, z, p^ ... Pn» We thus 
obtain n equations among the constants c,,...Csm. 

Substitute in these equations (f>^ for c,, ... <l>a^ for c^, and' 
we have n equations connecting a?j, ... a?,, £;, ^j,...^^, from 
which with the aid of the given equation ^j, ... p«, may be 
eliminated, and there will result a single equation connecting 
x^j ... Xf^ with z. This is the integral sought. 

[It appears from the manuscript that an example was to 
have been supplied here.] 

5. Cauchy's ijaethod is evidently a general one. But its 
generality is not of the same kind as that which belongs 
to Lagrange's solution of linear partial diflferential equations. 
It conducts us, not to a form embracing every possible 

B.D.E. II. 8 
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solution, "but to a system of results from which every possible 
solution may he derived, by arbitrarily varying the form of 
the function which expresses the initial state of the dependent 
variable, that is the value of z when x^ = 0, and then per- 
forming certain eliminations. To obtain a complete primi- 
tive we should only have to assume as the form of z when 
a7„ = a function of the variables a?j, ... x^^ involving n 
independent constants. The form of this function is arbitrary. 
Each distinct determination of it under the conditions leads 
to a distinct complete primitive. The number of such com- 
plete primitives is infinite. 

/ There are some most important problems in which the 
knowledge of a single complete primitive is all that is re- 
quired. For this purpose the method of Jacobi which we 
shall now give may be employed. 



JacoMs Last Method. 

6. Supposing z to be not explicitly involved in the given 
partial differential equation 

which we shall as before represent by F^ = 0, the problem of 
the discovery of a complete primitive consists in the finding 
of 71 — 1 equations 

such that between any two functions FiFj the relation 
[i^i^i] = (1) 

shall be identically satisfied. The values of ^^ ,...??» deduced 
from the equations, by rendering 

dz —pjdx^ — ... —pjl^n^ ^ 

integrable lead us to the complete rrimitive expressed by its 
inteo:ral 



^o' 
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Now the idea upon which Jacobi's later methods rest is 
that of directly solving the diiFerent systems of linear partial 
differential equations flowing from the general condition (1), 
not of solving, as in Cauchy's method, one of those equations 
and then limiting that solution by conditions which virtually 
involve the satisfaction of the others. 

71 (tI ~- 1^ 

It is evident that the entire series of — ^^— — - conditions 

(1) will be satisfied if we determine F^ to satisfy the single 
equation 

then F^ to satisfy the system of two simultaneous partial 
differential equations 

iF,F^ = 0, IF,F^ = 0, 

then F^ to satisfy the system of three simultaneous partial 
differential equations 

and so on, until finally F^ is determined by the solution of 
the system of n — 1 partial differential equations 

[F,F^]=0, [iy^J = 0, ...... [J^^,FJ=0. 

Now all these are particular cases of the general problem 
of determining a function P which shall satisfy simultaneously 
the equations 

[2^,P] = 0, [P.P] = 0, [P.P|=0. (2) 

F^y Pj, ... F^ being given functions between each pair of 
which the equation 

[F,F2 = 

is identicallk satisfied. Here P will represent in succession 
the series i^, F^, ... F^. 

The given system is one of homogeneous linear partial 
differential equations. It belongs to the class of systems the 

8—2 
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general theory of which is discussed in Chap. xxvr. But it 
is not necessary to apply the theory in its general form. We 
need only a single integral ; for a single value of each of the 
functions -Fi, jP„ ... i^ suffices in combination with the given 
value of -Fj for the determination of a complete primitive. 
Now it may be shewn that the system is oi the class dis- 
cussed in Chapter xxvi. If expressed symbolically in the 
form 

A,P=0, A,P=0,...AJ? = 0, 

the condition 

will be identically satisfied. Hence Jacobi's method for the 
treatment of systems of this kind may be applied. 

That the system is of the kind asserted is a consequence of 
the following proposition. 

Proposition. If the equations 

[wP]=0, [vP] = 
are expressed in the symbolic form 

AP^O, AT=0, 
then the derived equation 

(AA'-A'A)P=0 (3), 

will be equivalent to 



r[wt;]pl = 0. 

V A — "S*"* (— — — — ^ \ 

~ *"^ \dxi dpi dpi dxj ' 

*"^ \dxi dpi dpi dxJ ' 
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Hence since j- is the coefficient of -j- in A'P, and -;- 

its coefficient in AP, its co efficient in the derivedequation (3) 
will be (Chap. xxv. Art. 5), 

A ^^ — A' ^ • 



or 

or 

d 
or 



^i^n fdu dy du^ d^v dv d^u dv d*u \ 

"^ \dxi djpidxj dpi dxidxj dxi dp^dx^ dpi dxidxj ' 

d ^i^nfdu dv du dv \ 
dxj *^^ \dxi dpi dpi dxj ' 

d r T 

In like manner the coefficient of -^— is 

dXj 

Hence (AA' - A'A) P = 2^? (^-M ^ _ ^^^JO ^ 
^ ^ \ dXj dpi dpi dXjJ 

= [mp], 

whence the Proposition is established. 

Applying this to the system (2) we see that any derived 
equation will be of the form 



[[i^.2^p]=0. 



But \FiFJI\ = by the conditions given ; hence the condi- 
tion (A^Ay — AyA<) P = 0, is identically satisfied. 

The results of Chapter xxvi. being thus directly applicable 
to the system under consideration, we see that a common 
integral of the system (2) may be found by a ^ries of alter- 
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nate processes of integration and derivation. We "begin hj 
seeking an integral of the first partial differential equation. 
By a process of derivation, always possible, followed by the 
integration of a differential equation between two variables, 
we arrive at a common integral of the first two partial difie- 
rential equations. Again, by a process of derivation followed 
by the solution of a differential equation we obtain a common 
integral of the first three partial differential equations. And 
so on, until a common integral of all is obtained. 

7. Another solution of the above problem has recently 
been given. Beginning as in Jacobi's method by finding an 
integral of the first partial differential equation, a process of 
derivation agreeing in principle with Jacobi's, only more 
extended, may lead us without further integration to a point 
at which the discovery of a common integral of the entire 
system will depend only upon the solution of a single difife- 
rential equation of the first order susceptible of being made 
integrable by a factor. Failing this, it will enable us to 
convert the given system of partial differential equations into 
a new system possessing the same general character, but con- 
taining one equation less. Upon this the same process may 
be tried with a similar final alternative — and so on till the 
required integral is discovered. {On the Differential Equa- 
tions of Dynamics, Philosophical Transactions, 1863). 



V 
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CHAPTER XXVIII. 



PARTIAL DIFFERENTIAL EQUATIONS OF THE SECOND ORDER. 



[This Chapter is a reconstruction on a larger scale of part 
of Chapter xv. At the end of the Chapter reference will 
be given to other writings of Professor Boole on the subject 
here discussed.] 

1. The general form of a partial differential equation of 
the second order is 



where 






F{x 


.y. 


^,P, 


q, r, 8,t)=0 




■w. 


p. 


dz 


r- 


dz 


r- 


d'z 
'da?' 


d'z 
"'dxdt,' 


d'z 
*~dy'- 





It is only in particular cases that the equation admits of 
integration, and the most important is that in which the 
differeutial coefficients of the second order present them- 
selves only in the first degree ; the equation thus assuming 
the form 

^r + /S?+ r<= r (2), 

in which J2, 8^ T, and F are functions of x, y, Zyjp and q. 

The most important part of the theory of the solution of 
this equation is due to Monge, and was extended by Ampfere 
to the more general. equation 

Br + 88+ Tt+ U{8'--rt)=V (3). 
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This equation, together with the particular equation of 
Monge, and the equation 

both which though falling under Ampfere's general form 
possess peculiarities demanding special notice, 1 propose to 
consider in this Chapter. I shall in conclusion make some 
observations on the theory of partial differential equations of 
the second order with more than two independent variables. 

Monge's method, and Ampfere's in so far as it is an exten- 
sion of Monge's, consists in a certain procedure for discovering 
either one or two first integrals of the form 



^o* 



«=/W (4). 

u and V being determinate functions of a;, y, z, p^ and q ; and 
f being an arbitrary functional symbol. From these first in- 
tegrals, singly or in combination, the second integral involving 
two arbitrary functions is obtained by a subsequent inte- 
gration. 

Now this procedure involves the assumption that the pro- 
posed equation admits of a first integral of the form (4). But 
such is not always the case. There exist primitive equations 
involving two arbitraiy functions, from wnich by proceeding 
to a second differentiation both functions may be eliminated 
and an equation of the form (2) obtained, but from which it 
is impossible to eliminate one function only so as to lead to an 
intermediate equation of the form (4). Especially this hap- 
pens if the primitive involve "an arbitrary function and its 
derived function together. Thus the primitive 

-2^ = ^ (y + ^) + '^ (y - ^) - ^ {0' (y + ^) - ^' (y - «) } — (5) , 

leads to the partial differential equation of the second order 

'-< = ? (6). . 

but not through an intermediate equation of the form (4). 

It is necessary therefore, not only to consider the case in 
which the assumed condition is satisfied, but also to notice 
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what has been done in those cases which do not at present 
fall under the dominion of any known method. 



Oeneau of the Equation^ 

2. Prop. I. A partial differential equation of the first 
order of the form u =/(v), or its symmetrical equivalent, 

F{u,v)^0, 

in which u and v are any functions of a?, y, z, p^ j, always 
leads to a partial differential equation of the form 

5r+ /&+ Tt+U{s^-rt) = F. 

For, differentiating the proposed first integral with respect to 
af, and with respect to y, we have 

dF fdu du du .du\ 
du \dx dz^ dp dq ) 

dF/dv dv dv ,^v\_ 
dv \dx dz^ dp dq J" * 

dF (du du du ,^\ 
du \dy dz^ dp dq J 

dF/dv dv dv , ^^ .^ ^ 
dv \dy dz^ dp dq J " ' 

For brevity, write 

/du\ J. du du , /du\ j, du ^ du 

and then eliminating 

dF dF 
du' dv* 
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we have 

((du\ . du du \ {(dv\ dv dv \ 

(fdu\ du du ) (fdv\ ^dv dv \ 

which, on effecting the multiplication, gives 

(du fdv\ ^ /du\ dv) 
[^ [dyj \dy) dp) ^ 

(fdu\ dv_du fdv\ _ /^\ d^ du /^\1 ^ 
'^%ix)di) dp\dx) \dt/)dq dq\di/J) 

(/diAdv^^du fdv\\ . 
"*" 1 W d^ dq \dxj) 

fdu dv du dv\ , , . 
\dq^ dp dp dq) ^ 

=(D(S-©(D <')• 

a result which, since u and v are by hypothesis given func- 
tions of X, y, z, p, J, is seen to be a particular case of the 
general form (3). 

We may hence deduce also the conditions under which 
particular forms included in the general form (3) arise. Thus, 
in order that the equation u =f{v) may give rise to a par- 
tial differential equation of the second order of Monge's form 

it is necessary that the condition 

du dv du dv ^ 
dq dp dp dq 

should be identically satisfied. This requires, by Chap. ii. 
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Art. 1, that u and v, considered as functions of p and j, should 
not be independent. 

3. The geometrical relations of the equation (3) are also 
remarkable. It may in particular be shewn that an equation 
of this form will be satisfied by the equation of any surface 
which constitutes the envelope of any system of surfaces 
formed by the variation of three parameters in subjection 
to two arbitrary conditions. For let the common equation of 
the enveloped surfaces be 

^=/(^.y»«»*j<') (S)» 

the parameters a, 5, c varying in subjection to the conditions 

conditions whiph, determining b and c as functions of a, may 
be reduced to the form 

5 = ^ (a), c^yfr{a) (9). 

Now the values of p and q being the same for any point 
in the envelope as for the same point in the generating surface, 
we have for all such points 

dfjx, y, a. h c) ^_ df[x,y,a,h,c) . . 

P' d^ ' ^ dy ^^'^^' 

These two equations in conjunction with (9) enable us to 
determine a, J, c as functions of a?, y, 0, p, q. Let these 
values be 

Then substituting in (9) we have 

equations which hold for all such points. These are then the 
partial diflferential equations of the first order of the envelope. 

Now each of these equations is of the general form (4) ; 
whence by Prop. i. the partial difierential equatioa of the 
second order is of the form (3), as was to be proved. 
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Let US actually construct this equation. 

Differentiating the first of the equations (10) with respect 
to X and to y, and regarding therein a as a function of those 
variables, and b and c as functions of a, we have 

ds? \dadx dbdx da dcdx da J dx ' 

^^ ^y I /^^y ^y ^ dy dc\da 

dxdy \dadx dbdx da dcdx da) dy ' 

from which we readily derive 

f d\f\da f ^y^^^o 
V da? ) dy \ dxdy) dx 

Proceeding in the same way with the second equation of 
the system (10) we have 

V dxdy) dy \ dy^ ) dx 
Hence, eliminating -5- and -^ , we have 

dy* dxdy dx* dj^ dy* \dxdy) ' 

the equation sought. 

Comparing this with the general form (3) we have the equa- 
tions 

dy* _ dxdy _^ _y _ tfe* dy' \dxdyl 

ii~ s T~U~ V 
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whence eliminating ^"^j -^, and ^^ we arrive at the 

equation, 

iS*+4(Z7F-5r)=0. 

This then is the condition which must be satisfied in order 
that the equation (3) may admit of an integral representing 
the envelope of a system of surfaces in which three parameters 
vary in subjection to two connecting conditions. It is only 
proved however to be a necessary^ not to be a sufficient, con- 
dition. 



Solution of the equation 5r + & + jfl + U{8* — rt) = F, when 
a first integral of the form F{u, v) = 0, exists. 

4. In the following sections we propose 

1st. To shew that when a first integral of the above form 
exists, its discovery depends upon the solution of two simul- 
taneous partial differential equations of the first order re- 
solvable into linear equations. 

2ndly. To shew how from such first integral or integrals 
the second integral is to be obtained. 

Prop. ii. If the equation 

Rr+Ss-\-Tt-\- V{s*-rt)^V 

admit of a first integral of the form F{u, v) =0, in which u 
and V are functions of x^ y, «,^, y, then will F{u, v) considered 
as a function ofx, y, z, p, q, and represented as such for brevity 
by F satisfy the two partial differential equations of the first 
order, 

\dx} dq \dy ) djp \dx)\dy} 
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\dq/ dq dp \dp / 



in 



which 



, j^ (fdF\ dF , (dF\ dF\ ^ 

(dF\_dF dF fdF\ ^ dF dF 
\dxJ'~dx -^ dz^ \dy ) dy ^ dz ' 



Regarding the function F in the proposed integral i^=0 
simply as a function of x, y, z, p, q, we have 



ah 



dp dq 



(11). 



Kdif J dp dq J 

On the other hand, regarding i^ as a function of x, y, z, p, q, 
mediately through u and v, we have the system 

dF(/du\ du du ] dF (fdv\ dv civ \ ^ 
du\\dx)^^''^dq''}'^d^\\di)^dp'"^di''\~^' I 

dF (/du\ du ^duX , dF ((dv\ dv dv \ ^\ ^^^^' 

and these systems are equivalent 

dF 
Now if from the second of these systems we eliminate -j- 

dF 

and -^ , we obtain (Art. 2), a result which must be equiva- 
lent to the proposed partial differential equation, 

Itr + 8s+Tt+U{s'''rt)^r ..(13). 

This equation then considered as a relation between r, s, t, 
must be an algebraical consequence of the relations (12), and 
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therefore of the equations (11). If then we determine alge- 
braically two of the quantities r, s, t, (we select r, t) from the 
system, and substitute their values in (13), that equation 
ought to be satisfied independently of the value of the re- 
maining quantity s. Now supposing^ and q to be both con- 
tained in F, so that neither -r~ nor ^— vanish, we have 
from (11), 



f^U^« 



,__ \dxj dq 



)^f 



dp 



t = 



(dF\ dF^ 
\dyJ dp 



dF 
dq 



dF 



substituting which in (13) there results 

\dx I dq \ay J dp \dx J \dy J dp dq 

+ \e f — Y -8~ — T (—]' 
\ \dq ) dq dp \d-p J 

Now as this equation is to be satisfied in virtue of the con- 
stitution of iJ, 8y r, Z7, F, and the function F^ and indepen- 
dently of 5, both the coefiicient of s and the absolute term not 
containing s must be separately equated to 0, Thus F con- 
sidered as a function of a?, y, 0, p, q^ and containing p, q, at 
least must satisfy the partial differential equations 

j^(dF\dF_^^/dF\dF 
\dx J dq \dy ) dp 

u(—\ (^^ + F— ^- 
\dx J \dy / dp dq" 

\dq ) dp dq \dp ) 

^-{(S)f-(f)f}-» J 



(15). 
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This result may also be established by forming tbe equa- 
tions of condition which express the proportionality of 
B^ /S, ...F, to the corresponding quantities m the constructed 
equation (7). From these equations of condition it is actually 
possible to eliminate in two distinct ways the quantities 

f-T-j, \~r]i ;y-> J-? tlie result being the formation of two 

partial dinerential equations for u agreeing in form with those 
above given for F. (See the memoir Ueber die partielle Differ 
rentialgleichung... GreWs Journal^ Vol. 61.) The actual 
transition from the former to the latter rests upon the con- 
sideration that the equation F (w, v) = 0, when F is arbitrary, 
is not really less general than the form <E) [F[u^ ^), v} = 0, in 
which the ^ is arbitrary. And here u has been replaced by 
F(u,v). 

The only condition respecting the application of the above 

equations is that we do not admit any relations which make 

• t OfF di • 1 

either ^r- or -r- to vanish. 
djp dg^ 

5. Prop. hi. The solution of the system of partial diffe- 
rential equations established in the last proposition may in all 
ca^es be made to depend upon that of simultaneous linear partial 
differential equations of the first order. 

In demonstrating this proposition we shall consider first 
the case in which P'= 0, then the case in which F= 0, lastly 
the case in which neither of these quantities vanishes. The 
ground of this division will appear in the investigation. 

Case 1. Suppose Z7= 0. The equation then is of Monge's 
form, 

Rr+8s-{-Tt^V. 
The second equation of the system (15) becomes 

\dq J dp dq \dp J ' 
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and therefore breaks up into the equations 

dq ^dp ^ dq ^ dp ' 

Wj and m^ being the roots of the quadratic equation 

Rm^-Sm+T^-O (16). 

As each of the above constituent equations is of the form 

dF^ dF 
dq dp ' 

the system (15) may be reduced to the form 

\dx) dp \dy ) dp dp dp ' 

which breaks up into the equations 

f=0, i?.(f) + T(f) + F».f =0. 
ap^ \axj \dy J dp 

The former of these we must reject (Art. 4). There re- 
mains for the determination of -Fthe system of linear partial 
differential equations 



dF_^dF^^ 
dq dp 



MS^^f)^"-?- 



.(17), 



and there will exist either one or two systems included under 
this form, according as the roots of the quadratic (16) are 
equal or unequal. 

B.D.E. II, 9 
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Case II. Let F = 0. The system (15) then becomes 

\dqj dq dp \dp J 

Eliminate Uhj maltipljing the first equation hy 

{dxjdp '^[dyjdq ' 

the second tj 

/dF\ fdF\ 
\dx)\dy)' 

and subtracting; we obtain, after rejection of the common 
. , dFdF 

We shall put this equation in the place of the second equa- 
tion of the system. This we are permitted to do under the 
restriction that in seeking to satisfy the system so changed 
we do not make use of any relations which would cause either 
of the two factors employed in the process of elimination to 
vanish or become infinite. 

The new e(juation reduces to one equation, or breaks up 
into two equations of the form 



©-»©- (-). 



\dy. 
m being determined by the quadratic equation 
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Making (—J^m[-^J in the first equation of the system 
(15), we get 

which breaks up into 

But if we combine the first of these with (18), we obtain 

and this combination causing both the factors employed in 
the elimination of U to vanish must be rejectea. There 
remains then the combination 

and this will represent either one or two systems of equations 
according as the quadratic determining m has equal or un- 
equal roots. 

Case III. Let neither ?7= nor F= 0. 

Multiply the second equation of the system (15) by an 
indeterminate quantity Z, and add to the first; then we have 

\dx) dp \dy) ^ 

^('^-«)f f- w- 

9—2 
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We shall enquire whether it is possible so to determine I 
as to resolve this into linear factors. . 

We might investigate this by resolving the equation as 
a quadratic with respect to -j- or -7—. But the form of 
the equation suggests what the forms of the linear factors 
must be if the resolution be possible. For as the squares of -t— 

J-nf 

and -T- both appear, and these squares alone, in the fimc- 

flW flW 

tion .to be resolved, it is clear that -=- and 3— will be the 

dq dp 

only differential coefficients of F which will appear in both 

linear factors in common. The most general supposition 

possible is then that one factor shall contain t— and -j- with 
( — J , the other the same with /^ J . 

Assuming then one factor to be of the form 
,^F. dF ^ (dF\ 

it is seen from the. form of the coefficients of the first three 
terms of (20) that the other factor must be of the form 

^^.TldF^-UfdF\ 
dg[ m dp n \dyj' 

and the resolved form of (20) must be 

I rfj m dp n \dy J) \ dq dp \dx)) " 

Multiplying out and equating coefficients, we obtain the 
conditions 
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m 

n 

mn 

m 



The third and fourth of these conditions are equivalent, 

and give n = l. The first and second are also equivalent, 

T 
and give m — Yr These values reduce the last equation of 

condition to 

so that Z is determined by a quadratic. The resolved form 
of equation (20) now becomes 

dF T dF , (dl 



To these results we may give a somewhat simpler form by 
making Ul=m; not the m used above. We have then as 
the quadratic for determining wi, 

m'+ 8m-\'BT- UV^O (21), 

and as the resolved form of (20), 
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Let m^ and m, be the values of m. Then we have from the 
last the two distinct equations 

(j.dF dF jj(dF\\{ dF dF jj(dF\\_ 
[j,dF^ dF^„(dF\\( dF^^dF^„fdF\) . 

and it is evident that these will be together equivalent to the 
equations (15) from which they were derived. 

Now to satisfy these equations simultaneously it is neces- 
sary that we should equate to one linear factor from each of 
their first members, if we equate to the first linear factors, 
we have 

^dF ^ dF^^^(dF\ ^ 

whence, by subtraction, 

/ xdF ^ 

This combination must therefore be rejected (Art. 4). For 
the same reason must the combination formed by equating to 
the second linear factors in the left-hand members of the 
above two equations be rejected. There remains then only 
the combinations formed by equating to the first factor of 
one of these members, and the second of the other. 

Thus we should have the combination 
dF , dF . ^^fdF\ 



^ dc[ dp \doii J 



.(22), 
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with the combination which would be obtained from this by- 
interchanging m^ and m,. ' 

6. It results from the foregoing investigations that the 
function F is in all cases to be determined by the solution of 
two simultaneous linear partial difierential equations with 
Jive independent variables. Now the theory developed in 
Chapter xxv. shews that the number of integrals of such a 
system cannot exceed three. That theory enables us both to 
determine what the number of integrals is, and to construct 
the system of ordinary differential equations, reducible to the 
exact form, upon which their discovery depends. 

We have seen that the knowledge of two integrals u = a, 
v = b o{ the system enables us to construct a general first 
integral 

of the partial differential equation (3). And the solution of 
this first integral would lead us to the second integral which 
is the final object sought. But the direct solving of a partial 
differential equation of the first order which is not linear and 
which involves in its actual expression an arbitrary fimction 
is difficult, and happily it may be avoided here. The fol- 
lowing propositions will enable us to accomplish the virtual 
solution by a different solution, founded however upon the 
same general principles. 

7. Prop. IV. ITie integrals of the respective systems of 
simultaneous linear partial differential equations upon which 
the determination of F depends are so relaied that if from tu>o 
such respective integrals the values of p and q are determined, 
they will render the equation 

dz ^pdx + qdy 

integrahle* And in the particular case in which the two systems 
become identical, any two integrals of the system stand in the 
same relation. 
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For, let ^ be an integral of the system (22), and "9^ an 
integral of the associated system obtained by interchanging 
m. and m, in the case in which these quantities are different. 
Tnen <E> satisfies the equations 

and ^ satisfies the equations 

But the necessary and sufficient condition in order that the 
values of p and q derived firom the equations 4> = 0, ^ = 0, 
may render dz -^pdx — qdy integrable, is 

(d^\^_d^(d^\ 
\dxj dp dp \dx) 

. -(Df-^O-" «• 

See Chap. XIV. Art. 11, Equation (36). 

Now if from the previous equations we determine the values 

(d^ fd^\ /d;9\ /rf*^ 
\dxj' \dy)' [dxj' \dy)' 



of 



and substitute them in the above equation of condition it will 
be identically satisfied. 
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The determination of (^J > ••• from the previous systems 

requires that U should not vanish. Hence the proposition is 
established except in the case of Z7=0, which is left doubtful. 



To examine this case let us revert to the system (17) which 
is proper to it. To that system since 

whence Bm^m^ = T, 

we may give the form 

dF dF ^ (dF\ , fdF\ V dF ^ 

or the form obtained from this by interchanging m^ and «i,* 

Substituting in these respective forms <E> and "9 in succes- 
sion for F, we find 

d^_ d^ fd^\ _ f^\ V d^ 

d^'^'^'dp' \dx)'"^'^'\dy)''B ^/ 

rf* dV /<^^_^ i^^Z^ 



^"^«^' \dx)^ ^"^'Kdy/ B dp\ 

and these values substituted in (23) reduce it to an identity. 
Thus the proposition is established generally. 

Lastly, as in the case in which the two roots of the quad- 
ratic for determining m are equal, the two systems of partial 
differential ec|uations for determining <E> and "9 become one, it 
follows that if from two integrals of that one system we can 
deduce values of^ and j these values will render the equation 

dz-^pdx — qdy^^O 
integrable* 
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8. Prop. v. When the system of simultaneous linear par- 
tial differential eqiuitions determining F admits of two integrals 
u = a, t; = i, it will admit or will not admit of a third inte- 
gral w = Cj according as the roots of the quadratic determining 
m are equal or unequal. 

The system in question, (22), becomes when we divide by 
UdLnd write for (--}-A and (-j-) their full expressions 

dF dF m^dF R dF _ 

dF dF TdF m^dF_ 
dx^P dz^Ul^'^'U dq^^' 



or A,jP=0, A,jP=0, 



in which 



*~dx'^^da^XJdp^ U dg 
Hence the equation 

(A.A.-A.A.)2r=0 
becomes 

In this expression the coeflScient of the first term only has 
been calculated. 

Now, by the theory developed in Chap. xxv. in order that 
the two simultaneous partial differential ^equations should 
have their full complement of integrals (three) it is necessary 
that the above equation should be satisfied identically. This 
involves three conditions, namely, 
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7n^ — Wj = 0, 

the first of which is the one aflfaned in the Proposition to be 
necessary. 

Secondly, it is to be shewn that if this condition be satis- 
fied and it* the system of given linear equations admit of two 
integrals w — a, v = J, it will admit of a third. 

Replacing m^ and m, by m the system becomes 
dF dF mdF B dF_ 
dy'^^dz'^ U dp'^ U dq"^' 

dx^P dz^ U dp'^ U dq^^' 

Now if we construct from this the corresponding system of 
ordinary differential equations, we shall find it to be 

dz —pdx — qdy = 0, 
dp''^dx--^dy = 0, 

dq--jjdx- ^rfy = 0. 

Now it is impossible that the first of these equations should 
be integrated without a previous determination of p and q as 
functions of a?, y, «, seeing that dx, dy, dz are the three ditfer- 
entiala entering into that equation. Such determination can 
only come from the integration of the second and third equa- 
tions of the system. But if these equations can be integrated 
in the forms t^ = a, t? = J, then %i and v being particular values 
of F satisfying the partial differential equations, it follows 
from the last Proposition that the values of y and q which 
they will yield will make the first equation integrable. 
Hence if the system admits of two integi*als it will admit of 
three ; as was to be shewn. On the basis of these Proposi- 
tions the theory of the second integration rests. 
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Theory of the Second Integration'. 
9. First suppose the values of m unequal. 

Then Wi = ai, Vi = Ji being the two integrals (and we have 
seen that there cannot be more than two) of one of the systems 
of linear partial differential equations, and u^=^a^^Vt — h% those 
of the other, the general first integrals of the given system 
will be 

The values of p and j determined from these will by 
Proposition iv. render 

dz -^pdx - qdy = 

integrable, and the integral of this will be the general integral 
of the proposed partial differential equation. For it will in- 
volve explicitly or implicitly two arbitrary functions derived 
from those, in the first integrals. 

It suffices however, following herein Charpit's method, to 
combine one general first integral derived from the one sys- 
tem with a particular first integral derived from the other 
system, e.g. the integrals 

^K> Vi) = 0» ^a = «- 

The values of p and q hence derived, and employed as 
before, will lead to a second integral involving one arbitrary 
function and containing two arbitrary constants. This con- 
stitutes a complete primitive from which the general solution 
will be obtained by converting one of the arbitrary constants 
into an arbitrary function of the other, and eliminating the 
latter between the equation and the one derived from it by 
differentiation with respect to that constant. 

Secondly, suppose the values of m equal. 

In this case we have but one system of partial differential 
equations so constituted however that if it admits of two inte- 
grals it will admit of three. 
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Let w=a, v = J, w=^c represent these integrals. Then if 
from these we eliminate^ and q we shall obtain a final inte- 
gral of the form 

jand this constitutes a complete primitive from which we shall 
deduce the general integral by making J = 0(a), c ='^(a), 
and eliminating a betwBcn the equations 

Q^ df{x,y,a,^{a),'rlr(a)} ^ 
da 

To prove thi^ let tis combine the general and particular 
first integrals 

v = ^(m), w = a. 

The values of ;> and q hence obtained make 

dz-'pdx-' qdy = 

integrable, and the result can be no other than the remaining 
integral «? = c, or rather what this would become on eliminat- 
ing j? and q from it. But since the equations by which this 
integration are to be effected are equivalent to 

w = a, v^^{<i)y 

w will become a function of a;, y, z, a and <^ (a). Also by 
Charplt's method c is to be treated as a function of a, so that 
ultimately we have the result above assigned. 

We have here supposed U not to vanish. If it do th<e 
theory assumes another but simpler form. Let 

be the two general first integrals. Then, since by the con- 
•dition at the close of Art. 2, if p be eliminated from these 
equations q will also disappear, it suffices to eliminate them 
togedier in order to obtain the general second integral. 
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10. Although the cases in which U- and F= have in 
the foregoing sections been treated for simplicity apart, their 
theoiy might have been deduced from that of the case in 
which neither [7 nor F vanishes. 

Thus to deduce the equations for the case of Z7= elimi- 
nate from the general system (22) ^ and -j- in succes- 
sion, and we find 

(^^-."^f -*.(!)+ ^(S=«. 
(«^— .-jf+j^(f)-*.(f)-»- 

But from (21) RT- i»j»w, = UV. 

« 

Substituting, and then dividing by [7 we find 






the equation determining m^, m^ being 
m^A-Bm + RT^O. 
This is equivalent to the results of Art. 5, Case i. 

11. We found it necessary (Art. 3) in order that the gene- 
ral partial differential equation of this Chapter should be satis- 
fied by the envelope of a system of surfaces the equations of 
which contain three parameters varying under two conditions 
that the relation 

ig" + 4(Z7F-5r) = 

should be satisfied. 

It appears from Art. 8 that this is but one of three condi- 
tions necessary and together sufficient for this purpose. The 
formal conditions for every form of ultimate solution con- 
sistent with the existence of a general first integral i^(w, t;) « 
can be deduced in the same way. 
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[In the Bulletin de VAcadimte Imjp^riale des Sciences de 
St Pftershourg^ Vol. IV. 1862, there is an article entitled Con- 
sidSrations sur la recherche des mtegrales premieres des Equa- 
tions diffSrentielles partielUs du second ordre^ par O. Boldt 
(Lu le 7 Juin 1861). 

The article occnpies pages 198 — 215 of the volume. Al- 
though the name does not quite correspond, I consider that to 
be a misprint, and I attribute the article to Professor Boole, 
partly from the nature of the contents, and partly because it 
is known by his friends that he was engaged at a time corre- 
sponding to the date here given in the preparation of a mathe- 
matical article in French. 

The object of the article is to determine the conditions 
necessary for the existejice of a first integral of the equation 

where 5, S, T, and TFare any functions of a?, y, «, -j- and --r- ; 

and also to determine the conditions which must hold in order 
that Ampere's method of integration may be employed. 

In Crelle's Journal, Vol. LXI. there is an article by Pro- 
fessor Boole, entitled Ueber die partielle Differentialgleichung 
zweiter Ordnung Rr+ 8s+ Tt+U{f'- rt) = T. 

The article is dated 1862 ; it occupies pages 309 — 333 of the 
volume. 

Among Professor Boole's manuscripts I found a memoir 
very closely resembling the article in Crelle's Journal; it 
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would appear that the memoir was drawn up with a view to 
publication in the Transactions of some English Scientific 
Society, and that this design was afterwards abandoned in 
favour of the article in Crelle's Journal. 

After some hesitation I have resolved to print this memoir. 
Even if the memoir had been identical with the article in 
Crelle's Journal it would have been convenient to the English 
reader to be able to avail himself of the investigations ; and 
the memoir contains remarks which do not occur in the article, 
and which are interesting in connexion with the history of the 
subject. There is some repetition of matter which has alreadj^ 
been given in Chapter xxviii.; but I was unwilling to impair 
the completeness of the memoir by abridgment or omission. 
Accordingly the memoir forms the next Chapter of the present 
volume. 

In Article 2 of the next Chapter will be found the pro- 
cess to which there is an allusion towards the end of Article 4 
of Chapter xxviii. 

It is obvious that the subject of partial differential equa- 
tions of the second order was much studied by Professor 
Boole, The chronological order of his writings on the sub- 
ject appears to be as follows : 

1. Chapter xv. of the first edition of his work. 

2. The article in the Bulletin of St Petersburg. 

3. The memoir which forms Chapter xxix. of the pre- 
sent volume. 

4. The article in Crelle's Journal. 

5. The Chapter xxviii. of the present volume.] 



Digitized by 



Google 



( 145 ) 



CHAPTER XXIX. 



ON THE SOLUTION OP THE PARTIAL DIFFERENTIAL EQUATION 

Br + 88+Tt+ U{/-rt) = F, in which J2, 8, T, U, V 

ARE GIVEN FUNCTIONS OF X^y^ «, p, j. 

1. The equation, the theory of the solution of which 
I propose to consider in this paper, is remarkable from its 
connexion with Geometry. If the equation of a surface 
contain three constants which vary as parameters in sub- 
jection to any two conditions connecting them, the gene- 
rated envelope will satisfy a partial differential equation of 
the above form. In other words any envelope of the surface 

i?'(aj, y, «, a, 5, c)=0 

formed by the variation of a, &, c in subjection to two con- 
necting conditions 

^i(a, J, c)=0, ^j,(a, J, c)=0 

is necessarily an integral of a partial differential equation of 
the form given above. 

Now this theorem is the more important, because it is 
only when three parameters in the equation of a surface 
vary in subjection to two relations that the envelope pos- 
sesses, irrespectively of the form of the connecting relations, 
any definite character. If there be but one connecting rela- 
tion it is possible to determine that relation so as to make 
the envelope assume the form of any surface whatever, and 
therefore tne possible system of envelopes is in such case 

B.D.K II. 10 
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unlimited. If there be three connecting relations the para- 
meters become absolutely constant and no envelope exists. 

The partial differential equation 

is remarkable also as including all the cases in which a 

Srtial differential equation of the second order admits a 
St integral of the form 

w=/(v), 

u and V being definite functions of a?, y, z^p^ y, and /(f) 
arbitrary in form. 

Neither of these statements is sufficiently general to con- 
stitute a theory of the genesis of the partial differential equa- 
tion under consideration, but the second one is more general 
than the first, and is indeed sufficiently so to serve as the 
ground of an investigation which connects the solution of 
the equation in all cases with the satisfaction of a system 
of simultaneous ordinary differential equations of the first 
order and degree. And this is the ground upon which the 
method of the papei: will rest. I propose to shew, 1st that 
the solution of the given equation on the assumption that 
a first integral of the form u=f{v) exists requires the satis- 
faction of a system of two partial differential equations of 
the first order and second degree ; 2ndly that this system may 
be resolved into four systems, each consisting of two partial 
differential equations of the first order and first degree, two 
of which systems are irrelevant and the other two relevant ; 
3rdly that the solution of the two relevant systems ulti- 
mately depends on the solution of a system of ordinary 
differential equations of the first order, and that from these 
ordinary differential equations the given equation of the 
second order may be deduced independently of the assump- 
tion above mentioned. 1 shall also discuss the theory of the 
second integration. And I shall exemplify another method 
of solution connected by a remarkable law of reciprocity with 
the above method. 
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First Investigation* 

2. Prop. I. lf'^=f{^) 5« a first integral of the equation 

Br + 8s+Tt+ U{/^rt)^V, (1), 

then will u and v, considered as Junctions of a?, y, z, p, q, each 
satisfy two partial differential equations of the form 

\dxj dq \dyj dp 



+KS©^'' 



du- du 
dp dq 



= 



.(2), 



tn 



, . 7 fdu\ J fdu\ ^ J ^ du ^ du , 
which ( -r- 1 and ( -r- 1 stand jor -^ — h^ -3- , and 



du du 
dy'^^Jz 
respectively. 

To demonstrate this proposition we shall form directly the 
partial differential equation of the second order of which 
u=f{v) is an integral and, comparing that equation with (1), 
deduce the conditions for the determination ot u and v. 

Differentiating u=^f{v), first with respect to x and secondly 
with respect to y, we have 

du du dz du dp du dq 
dx dz dx dp dx dq dx 



•' ^ ' [dx dz dx dp dx dq dx) ' 



du du dz du dp du dq 
dy dz dy dp dy dq dy 



=/'w{; 



6?v dz 
dz dy 



dv dp 
dp dy 



' "-' ^ dq dy) * 



10—2 
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.- ^ du ^ du ^ fdu\ du du , fdu\ 

or. if we represent ^+i> ^ by y , ^+y^ by y , 

|byr,|andgby«, and|by«, 

/rfw\ . du ^ du ^f. . {fdv\ , dv , ^dvl 

Eliminating /'(t?) we arrive at the partial differential equa- 
tion of the second order, 

{du (dv\ ^ dv f^W 
dp \dy) dp \dy)] 

(du fdv\^ dv /d^\ dv /du\ ^du fdv\) 
^[d^\dy) dq\dy)^dp\dx) dpKdxJV 

(dv fdu\ du (dv\\ (du dv du dv\ , , . 

-(!)©- ©(I) ("• 

It is seen that as respects the mode in which the quan- 
tities r, 5, t are involved this equation is of the same form 
as the given equation (1). That it may be equivalent, its 
coefficients must stand to those of (1) in a common ratio /k. 
This gives 

du fdv\ dv fdu\ ^ , . 

^ tej -^ \d^) ^f"^ ('*)' 

du fdv\ ^ dv fdu\ dv_ fdu\ du fdv\ _ ^ „. 
di \dy) dq_ \dy)^dp \dx)''dp\^) "^^ ^^^' 

dv /du\ du fdv\ rn / X 

dk u -^ w "'^^ ^'^' 
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du dv du dv ^ .j. 

di Tf-Tr If^- -■^^' 

©(D-d)©""^ W- 

As we have here five equations which are homogeneous with 
respect to the four differential coefficients of v and to /[t, it is 
clear that we can, by the elimination of these (]^uantitie8, 
obtain a relation connecting the differential coefficients of ti 
with J2, 8y Ty &c. But the peculiar cyclical form of the 
functions in the first members of the above system enables 
us to effect this •elimination so as to lead to two final equa* 
tions independent of v and fi. 

^^ ^y (^) (^) ' *°^ ('') ^y ^ ^' *"^ adding, we find, on 
rejecting tbe common factor /», 

MS) 1+^(1)1^^© (I) 

Again, multiplying (a) by (^) , (5) by (^) (^), (c) by 

(^j, and (e) by (^j^-H-^ ^, addmg, and again re- 
jecting^the common factor fi^ we have 
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Hence, u considered as a function of x, y, z, p, q satisfied 
the two partial differential equations (4), (5), both which are 
of the first order and second degree. 

As u and v enter symmetrically into the system (a), (J), 
&c., V will also satisfy two partial differential equations of 
the same form, viz. the equations 

\dx/ dq \dyj dp 

^(D"-<l)(|)-^(|)' ■ 



+ F 



{ldv\ dv fdv\ dv] _ . 



(6). 



Further, these two systems of equations constitute the 
complete system of equations resulting from the elimination of 
/jL from the five equations (a), (J), (c), &c. ; for in their deter- 
mination, no factor involving either the differential coefiicients 
of u and v, or the quantities jB, 8, T, &c. has been rejected 
directly or indirectly. 

I am not aware that the above results of elimination have 
been noticed before. 

3. Prop. ii. The system of partial differential equations 
above obtained for the determination ofuy viz. 

'du\ du rr/^^\ ^^ 



\dx) dq \dyj dp 



•m 
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admits of resolution into fimr systems^ each consisting of two 
linear partial differential equations of the first order. Of these 
systems two only are relevant to the solution of the problem. 

For, multiplying the second by an indeterminate quantity 
X, and adding the result to the first, we have 

p/<?M\ du ydu du m(^:U\du 
\dx) dq dp dq \dy) dp 

*^B)tAtM-'' '"• 

Now let us see if it is possible to determine X so as to 
make the first member of the equation resolvable into linear 
factors. We cannot say d, priori that such resolution is pos- 
sible as we should be able to do if that member were homo- 
geneous and of the second degree with respect to three instead 
of with respect to ih^four subject variables 

(du\ /du\ fdu\ fdu\ 
d^J' W' \d^)' W' 

Observing that the squares of ^ and -j- are wanting 

in the first member of (8) while those of (-^j and f -j-j 

appear, we are led to assume as the proposed equivalent of 
that member an expression of the form 

Multiplying the factors of this expression together and 
then equating the coefficients with those of the first member 
of (8) we have 
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Jtm+\m== U+8\ (a), 

mm' = T\ (J), 

\n=iW=mn! (c), 

Iin' = B {d), 

nm':=T (e), 

nn'=V. (/), 

From (J), (c), (cZ), we find 

T 
w=:F, w' = l, m = \V, m' = -p., 

values which will be found to satisfy (e) and (/) also, and 
which reduce (a) to the form 

FV-flfF\ + i?r-Z7y=0. 
Supposing \ thus determined, the equation (8) becomes 

The result is a little simplified if we retain m in place of \. 
We thus find as the resolved form of the given equation 

m being determined by the quadratic 

m*-8m+BT--UV=^0. 

If wjj, Wj be the values of m thus found, we have 
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and these two equations are manifestly together eqiial to the 
system (7). 

Now these equations can only be simultaneously satisfied 
by equating to 0, one factor in the first member of each ; and 
fhe different combinations which are thus possible give rise 
to four binary systems of linear equations. Let us examine 
these systems separately. 

If we simultaneously equate to the two first factors of 
the left-hand members of the last two equations, we have the 
systems 



^©-».(|)-4;-». 



a system which, when m^ and m^ are different, is reducible to 
the system 

It is clear that this cannot lead to a value of u satisfying 
the given differential equation (1), because it takes no 
account of the forms of 8j Z7, and T. Indeed if we actually 
eliminate 



/du\ /du\ du du 
\dx)' \d^)' d^' Tq^ 



from the above equations by means of the system 
/dv\ du , du 



/au\ du du 

(du\ du , du ^ 



.(10), 
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(derived from the assumed first integral u =f(v) hj making 
f{v) = c, and difierentiating the result first with respect to a?, 
then with respect to y), we find as the result 

Again, if we equate to the two last factors of the right- 
hand members of (10), we have 



which, if Wj and w», are different, reduce to 

And it is evident that neither are these equations consistent 
with the given equation (1), because they take no account of 
S, Uy and B. The equation of the second degree to which 
they actually lead is 

There remain then the two systems formed by combining 
the first factor of each one of the first members with the 
second factor of the other, viz. 






(11). 



.(14). 
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That these systems are relevant to the solution of the pro- 
blem under consideration may be shewn by eliminating from 
either of them by means of (10) the quantities 

(^\ /du\ du du 
diJ' \d^)' Tp' dq' 

The actual result will be 

F{i2r + &+31f+tr(s»-rO-F} = (13), 

which, except in the particular caae of F=0, reduces to the 
given equation. 

More generally, if in the equation 

u^f{v) 

u and V are any distinct solutions of the system (11), the 
same result of elimination may be deduced. For v by hypo- 
thesis satisfies the equations 

fdv\ . /dv\ rrdv 



■^(i)+^(l)^''l-»- 



Subtract these equations multiplied hjf'{v) from the corre- 
sponding equations of (11), and representing U'-f^v) by W, 
we have 

which being of the same form as (11) it follows that 

Tr=0 or w-/(v)=0 

also leads to thp partial differential equation of the second 
order (13). 
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4. Prop. hi. To reduce the determination of the first inte- 
grals of (1) to the solution of a system of ordinary differential 
equations. 

Each of the systemB (11), (12) presents u as satisfying 
simultaneously two linear partial differential equations of the 
first order. 

To deduce the value of u thus conditioned it will obviously 
suffice to multiply in each system one of the partial differen- 
tial equations W an indeterminate multiplier \, to add the 
result to the other equation so as to form a new equation 
which will, like those nrom which it is formed, be linear and 
of the first order, and which on account of the indeterminate 
character of \ will be equivalent to the two. From the 
auxiliary equations which we obtain in the process of solu- 
tion, \ must be eliminated. 

If in this way we combine the equations of the system (11), 
we have, on arranging the resulting equation according to the 
differential coefficients of Uy 



•|jZi, + »»^ + \(2'j-H».^)}J 



dp dq 

Hence we have the auxiliary equations 

dx _ dy _dp ^dq & 

jB + XWg"7Wi + \r"T'"'\F"jBp + wiig' + \(rj + w,p)' 

and it is to be remembered that m„ w, are the roots of the 
equation 
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15? 



.(I). 



Eliminating \ from the. first four of the above equations we 
have 

Udq^ + m^dx - Rdy = 0] 

Udp + m^dy-Tdx^O^ 

dz — pdx — qdy = o) 

This then is the system of ordinary differential equations 
deduced from (11) upon the integration of which the determi- 
nation of u will depend. 

A similar system, differing from the above only in the 
mutual transposition of m^^ and w,, is given by (12), viz. 

Udq + m^dx - Rdy = 0] 

Udp + m.dy-'Tdx^O) (II). 

dz — pdx — qdy =J)i 

If from either of these systems we can deduce two inte* 
grals of the forms 

it is obvious, from what precedes, that 

will constitute a first inte^al of the proposed (1), and there 
being two systems in question, two such first integrals, each 
invomng an arbitrary constant may coexist. 

6. Pbop. IV. To deduce the second integral of (1). 

It will be necessary to consider separately the cases in 
which m^ and m^ are equal and unequal. 

First let m^ and m^ be equal. 

Both the systems (I), (II) reduce to a single system which 



may be expressed in the form 



in th< 



dp^Yi^^-^Tr^y 



V 

m 



U 



dq = '—Yjdx+ Rdy 
dz = pdx 4- qdy ^ 



.(14). 
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Now, since the condition ;/ =^ is here satisfied, it is 

manifest that if from any two mtegrals of the above system 
of the forms w = a, v = b, simultaneous values of p ana q *be 
determined, these values vdll render the third equation oi the 
system integrable, and the efiect of its integration will be 
virtually to determine « as a function of a?, v, and three arbi- 
trary constants, viz. a, J, and a constant c introduced in the 
last integration. Let us represent the result in the form 

2 = <^(aj, y, a, h, c) (15). 

Now what relation will this result bear to the general solu- 
tion of the partial differential equation given, to the solution 
which we should obtain by integrating, not the particular 

Xations w = a, v = J, liut the general first integral w=/(v), 
ch includes them both. 

To integrate the equation u =f{v) it suffices to deduce any 
particular ecjuation involving an arbitrary constant b, which, 
m conjunction with u —f{v) will render 

dz —pdx — qdy — 

ntegrable, and to integrate the last equation regarding the 
arbitrary constant of integration as an arbitrary function of 
h. The result is a complete primitive in which, by the 
variation of & as a parameter the general integral is implicitly 
involved. 

Now either of the equations w = a, v = J will, in conjunc- 
tion with u =y W determine p and ^ so as to make 

dz ^pdx — qdy = 
integrable. Take the equation v = J, then u =/(v) reduces to 

IS, in pli 
1, we ha 



Thus, in place of the equations w = a, t? = J, of the previotls 
section, we have 



Digitized by 



Google 



ART. 6.] OP THE SECOND ORDER, 159 

for the determination of jp and q. The constant c introduced 
in the final integration becomes also, according to the above 
theory, a function of J, and the complete primitive is of the 
form 

z^<l>{x, y, 5, /(J), yjr (J)} (16), 

while the general integral is found by eliminating b between 
this equation and its differential with respect to b. 

The general integral therefore represents the envelope of 
the surface represented by (15), a, J, c being parameters sub- 
ject to any two connecting conditions. 

As Wj, m^ are supposed equal, a necessary condition of the 
possibility of this species of integration is that 

)g»-4(JKr-i7F)=0 (17), 

or 

the value of m is --- , and the system (14) reduces to 

Udq + jdx-'Jidt/^-O ^ (^^)' 

dz — pdx — qdy = . 



We conclude therefore that if (17) be satisfied and we can 
from (18) deduce a value of z in terms of a?, y, and three 
arbitrary constantly the equation expressing that valvs will be 
a complete primitive^ and the general integral will be found 
by making the constants vary in subjection to two arbitrary 
conditions* 

Ex. Let the given equation be 

xgr + ypt + xy{s^— rt) =pq* 
Here R — xq^ S=Oy T^ypj U-xy, V^pq. 
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The condition (17) is satisfied, and (18) becomes 
xydp—ypdx^O, 
xydq — xqdy = 0, 
dz "pdx — qdy = 0. 

From the two first of these we find 
p=^ax, q==by, 

whence from the third, 

This is the complete primitive, and the general primitive 
consists of all possible equations derived from this by making • 
a, b, c vary in subjection to two conditions. 

Ex. 2. Given 



(l + 2»)r~2p25 + (l+/)<-- 



= -(!+/ + 2^'- I 

Here the equation for m reduces to 

m^-{-2pqm+pY-0y 

I 
whence wi = — 2^, and the system (18) gives I 

(l+y + j';* I 
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Subtracting the upper equation multiplied by pq from the 
lower one multiplied by l + 2^ and dividing by l+i?' + J^ 
we. haye 



whence 



*+V(i+/+?^""' 



i h. 



In like manner, 

Hence determining p and j, 

, _ {x-a)dx + {y-l)d y 

vii-(^-ar-(y-*)r 

Therefore (a? - a)'* + (y - S)' + (« - c)' = 1. 

^rom this form of the complete primitive it is evident that 
the general integral will represent all possible tubular surfaces 
formed by the motion through space of a sphere of constant 
radius unity. 

Secondly, let m^ and m^ be unequal. 

Then since, in neither of the systems (I) and (II) is the 

condition -f- = -r satisfied, from neither system separately 

can values oip and q be obtained which make dz=pdx + qdy 
integrable. 

But, as will be shewn, any two integrals obtained, the 
one from the one system and the other from the other, will 
give values of p and q which will render dz=:pdx + qdy in- 

B.D.E. II. 11 
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tegrable, and the general solution will consist of all possible 
integrals of the latter equation thus obtained. 

Or if the complete first integral of either i^stem be com- 
bined with any particular integral involving an arbitrary- 
constant obtained from the other, the two will fiimish values 
of ^ and J which render dz=pdx-\-^dy integrable, and its 
integral will be a complete primitive involving one arbitrary- 
function in, its expressed form, another in the connexion of 
its two constants ; the general primitive being found in the 
usual way by making the constants vary as parameters in 
subjection to a. single arbitrary connecting condition. 

In fact it may be shewn that if we attempt by the process 
of Charpit or Lagrange to integrate the partial differential 
equation of the first order u=f{v), deduced we will here sup- 
pose from the system (i), we virtually construct the system 
(ii) in the auxiliary equations upon which the process of 
solution turns. I have obtained a direct proof of this proposi- 
tion, but I think it preferable and at the same time suflScient, 
to direct attention to the prior ground upon which it rests in 
the relations of the systems of partial differential equations 
(11), (12) from which the systems of ordinary differential 
equations (i), (il) are derived. 

Let P=0 represent any integral of the system (11), and 
Q = any integral of the system (12). Then we have 



-f--(S--(f) 
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/Hence we deduce . "^s^ 

\dp \dx) dp \dx) dq \dy) dq \dy)) 

={*©-. (D)(S)-{-©--(f)}© 
M-(S-Kf))(f)-WS)-Kf)}(f)- 

The second mem'ber of this equation is identically 0, Hence 
dividing by F we have 

fm.^(f)^fff).^Q(f.yo (19). 

dp\dx) dp\dxj dq\dyj dq\dyj ^ ' 

But this is the known condition under which the values 
of p and J deduced from the equations P = 0> Q = make 
dz —pdx-\-qdy integrable; see Chap. xiv. Art. 13, Equation 
(36). 

We conclude then that if from the systems (i), (ii) we can 
dedtice two corresponding sy sterns of integrals 

then will the first integrals of (1) he 

while the second integral will consist of all possible relations oh^ 
tained either \st hy specifying the forms off , f and obtaining p 
and q a;s functions of x and y and integrating dz =pdi: + qdy^ 
or 2ndly, by specifying one of the functions jf[y f^, leaving the 
other arbitrary y determining p, q, integrating dz = pdx + qdy^ 
and regarding the final constant of integration as an arbitrary 
parameter. 

11—2 
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.*^ Ex. Given ar + bs + ct + e{8^-rt) = h, the coefficiefJ® 
being constant. 

Here-B=a, 6^=5, T==c, U=e, 7= A. 

Hence m,, m^ are the roots of 

m' — iw + ac — eA = 0, 

and the systems (i), (ii) give 

edq + m^dx — ady = 0] 
edj[> + m^dy — cdx=sO) 

edq^ + m^dx — arfy = 0) 
edp + mj(% — cdx = J * 

Whence the first integrals are 

eq + mjc- ay =^f^{ep-\- m^y "- ex), 

from which all possible second integrals are to be derived in 
the modes above explained. 

Let us take the second of those modes and give to the 
second of the above first integrals the particular form 

6p + mjf --cx= Gf 

G being an arbitrary constant. From this, and from the 
other integral, left in its complete form, we have 

cx-m^ + G ay- m^x+f^^ {(m, - m^) y + C} 

P- e ' *"■ e ' 

whence, substituting in the formula dz^pdx-^-qdy, inte- 
grating, replacing the arbitrary form 

r 
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^id introducing an* arbitrary function of G for the arbitrary- 
constant, we have 

for a complete primitive. The general primitive . consists of 
all possible relations obtained by eliminating G between the 
above equation and 

0=i[a; + f{K-m.)y+(7}] + ^'((7), 

when the forms of ^ and -^ are specified. 



Second Investigation. 

6. If from the equation 

Rr-\-88-\-Tt'\-U{8^-rt) = V. (20), 

we eliminate r and t by means of the equations 
dp = rdx + «dy, 
dq = sdx + tdy^ 
the result will be 

[i?d/ - Sdxdy + Td^ - U(dpdx + dqdy)\ s 

^Bdpdy+Tdqdx-Udpdq-Vdxdy (21). 

There are different considerations (all of them however in- 
volving, as I have been led to think, a more or less explicit 
reference to some theory of the genesis of the given partial 
differential equation) which indicate that its solution depends 
upon that of the equations obtained by equating to the part 
affected and the part not affected by «, viz. upon the solution 
of the equations 

Bdy' - Sdxdy + Tdx"" - U {dpdx + dqdy) = (22), 

Edpdy-^r Tdjdx- Udpdq - Vdxdy =^ (23). 
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Without entering into these considerations let us inquire 
what consequences may be deduced from these equations 
assumed to be true. 

It is seen that these equations are connected by a remark- 
able reciprocity with the partial differential equations (7). 
They will in fact be converted into these equations if we 
change 

dx, dyy dpy dq, UIV, S (24), 

into 

du du du du ^ ^ ^ . 

respectively. From ihis formal connexion it follows that if 
we multiply (22) by \ and add to the result (23), we shall be 
able to determine \ so as to permit the resolution of the equa- 
tion thus formed into linear factors. Ultimately we shall, as 
appears from Art. 3, reduce the system (22), (23) to an equi- 
valent system of the form 

{-Rdy-m^dx + Udq) (- m^dy- Tdx + Udp) = 0, 

{-Rdy-m^dx'\' Udq) {-m^dy- Tdx+ TJdp) = 0, 

m, and m^ being determined by the equation 

or, changing the sign of w, 

(- Edy + m^dx + TJdq) {m^dy - Tdx + Udp) = 1 . . 

(- Rdy + mjdx + Udq) {m^dy - Tdx + Udp) = ol 

m^ and m^ being as in the former investigation roots of 

m^-Sm'\-RT'-UV=0. 

Equating to the corresponding factors of the first mem- 
bers we have 



'-Rdy-{'m^dx-\- Udq^O 
— Rdy + m^dx + Udq 
m^dy — Tdx + Udp = 



= 0) 

= oi 



m^dy — Tdx + Udp = 
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The first of these, m^ and m^ being different, is resolvable 
into 

Udq-Bdy^O, dx = 0; 
the second into 

Udp-Tdx = 0, dy = 0, 

and it is obvious that neither of these can lead to the given 
partial differential equation (1). The first of them combined 
with the equations 

dp=rdx-\-8dy^ dq = sdx -\- tdi/ (27), 

leads in fact to the partial differential equation 

B- m=:0 '. (28), 

the second in like manner leads to 

T-Ur = (29). 

But equating to the non-^corresponding factors of the first 
members of (26) we have 



• Rdy + m^dx + Udq = 
m^dy — Tdx + Udp 



m^dy - Tdx+ Udp = 
— Bdy + m^dx + Udq^ 






Now these systems when completed by the equation 
dz = pdx + qdy agree with the systems (i), (ii) deduced in the 
previous investigation. 

It remains to shew that these systems actually lead to the 
given partial differential equation (1) directly. Eliminating 
from either of them, combined with the system (27) the dif- 
ferentials dxj dy^ dp, dq, we shall have as the result 

UlEwSs+Tt+UfZ-rt) - F1 = (30), 

which, rejecting the factor CT, as from (13) we rejected F, is 
the differential equation proposed. 
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Ground of the Reciprocity above noticed, 

7. The reciprocity above noticed is not of a primary cha- 
racter, but is foundea upon two prior laws which I shall pro- 
ceed to demonstrate. 

If from the partial differential equations of the system (7) 
we eliminate Fand substitute the resulting equation in the 
place of the first equation of the system we shall obtain the 
equivalent system 

H (^ _ S— —+T (—\ 
\dq) dq dp \dp) 

-"{(Sl+d)!)-"] 

These equations are both symmetrical and it will be ob- 
served that they are convertible the one into the other by 
changing 



..(31). 



into 



du du du du ~ ^ -^ 
dq' 'dp' ^' Ty' ^' ^' ^ 



du du du du ^r a TT 

Xr» :^> :^i ^9 "^> ^ 



(32), 



dx' dy' dg' dp' 



(33), 



respectively. This is a law of reciprocity which connects 
solely the differential coefficients of u and the coefficients 
TJ^ S, V of the original equation. 

Again w =0 is by hypothesis a solution of the given partial 
diiferential equation. Regarding it however simply as an 
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equation which is true and the truth of which is consistent 
with that of the equations 

dp = rd^ + sdy\ ^^^^^ 

dq = sdx + tdy) 

and differentiating it first with respect to a?, secondly with 
respect to y, we have 

du du dz du dp du da 
dx dz dx dp dx dq dx ' 

du du dz^.du dp du dq ^ 
dy dz dy dp dy dq dy^ ' 

equations to which we may give the form 



_(diA_ 
[dx)" 

fdu\_ 
KdyJ" 



du du 
dp dq 



du\ du du 
dp dq* 

Now this system is of the same form as the system (27) 
and will agree with it if we change 

_du ^du du du . . 

1^' Ty' dp' dq ^^^^' 

into 

dp^ dq^ dxy dy (36), 

respectively — a change which does not affect the coefficients 
of the given equation, and which is therefore the expression 
of a law of reciprocity distinct from that last noted. The 



of the given equation, and which is therefore the expression 

" reciprocity distinct from that last noted, 
combination of these two laws does however lead to the 
law exemplified in the researches of the previous Article ; 
see (24) and (25). 

The question here arises whether it would not have "been 
tetter to employ from the first the symmetrical foims (31) of 
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the partial differential equations of the first order and second 
degree upon which u depends, than the unsymmetrical forms 
(7). It was indeed from the symmetrical forms that the 
chief results of this paper were originally obtained, but the 
unsymmetrical forms lead to the same end in a simpler way, 
and therefore they have been made use of in the present 
memoir. 

It may be proper to notice, in concluding this section, that 
the symmetrical forms in ordinary differentials would have 
emerged in place of the unsymmetrical ones of (22) and (23), 
if the quantity s* — rt had been retained instead of 8. The 
equations 

dp = rdx + sdy^ dq=sdx + tdt/j 

enable us in fact to reduce the given equation (20) to the 
form 

Bdp" + Sdpdq + Tdf - Vidpdx + dqdy) 

= {/--rt) {Edf- 8dxdy+ Tds? - U[dpdx + dq^dy)]. 

Hence arises the symmetrical system 

Rdf + Sdpdq + Tdg" - V {dpdx -^dqdy) = 0, 

Bdy^ - Sdxdy + Tdx^ - JJ {dpdx + dqdy) = 0, 

which is connected with the system (31) by the single law of 
reciprocity expressed in (35) and (36). 



Postscript, 

8. At the time when the above investigations engaged my 
attention I was totally unaware that the subject of them had 
been discussed by Ampfere {Journal de VEcoh Folytechnique, 
Tom. XI.) and recently by Professor De Morgan ( Cambridge 
Philosophical Transa^tionSy Vol. IX. Pt. IV.). I feel it there- 
fore incumbent upon me to state why after acquainting my- 
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f5elf with the results of their labours, I offer this paper for 
publication. 

The method of Professor de Morgan so far resembles the 
first method of this paper, and that of Ampfere the second, 
that while the former makes the solution of the problem 
depend directly upon that of simultaneous partial differential 
equations of the first order, the latter makes it to depend 
directly upon the solution of simultaneous ordinary differ- 
ential equations of the first order. The formal connexion of 
these methods by the law of reciprocity is, I believe, esta- 
blished for the first time in this paper. The system of partial 
differential equations of the second degree (7) has not, so far 
as I am aware, been given before. 

But a point which I think of deep importance is the follow- 
ing. By connecting, as in this paper, the differential equa- 
tions of the second degree, whether ordinary or partial, by an 
indefinite multiplier which is afterwards determined so as to 
admit of the resolution of the system into its component linear 
elements, we assure ourselves that each step of the solution 
offers a complete sequence to that which has gone before, and 
it only remains then to separate the different elements and 
determine whether they are relevant or irrelevant to the end 
in view. That any such distinction exists has not, so far as 
I am aware, been noticed before. And it seems to me the 
more important that it should be noticed because the solution 
of partial differential .equations in cases far more general than 
those above considered seems to depend upon the satisfaction 
of simultaneous differential equations of a degree higher than 
the first. I have in fact by an application of the Calculus of 
Variations arrived at the conclusion that the theory of the 
solution of all partial differential equations of the second 
order, whatever the number of variables may be, is very inti- 
mately connected with the satisfaction of a system of differ- 
ential equations of the type 

dF j^ dF . , dF ^^ ^ 
-df--^dxdy + -^dx'^0, 

F= representing the given partial differential equation, x 
and y any two of the independent variables, and r, 5, t the 
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second differential coefficients of the dependent variable with 
respect to x and y. 

I may perhaps at some future day resume the subject, to- \ 
gether with an inquiry into the theory of the solution of the 
partial differential equation of this paper, when the conditions 
under which the auxiliary equations (i), (ii) are supposed to 
be integrable are not satisfied. 

9. Note. It may be desirable to establish directly the 
converse form of one of the results of Proposition iv. For 
this object we shall shew that the equation of the envelope of 

2; = ^ (a?, y, a, b, c) (1), 

where a, b, c are connected by any two conditions of the 
forms 

will satisfy a partial differential equation of the form 

Br+83+Tt+U{s'--rt)=r (2), 

in which also 

Differentiating (1) we have 

___ d(f> d^ da d<f} dh d(f> dc 
^ clx da dx db dx dc dx 

_ d<f> d<f> da d4> dh d^ dc 
^ dy da dy db dy dc dy 

and by the nature of an envelope these reduce to 

P-% ?=f : (^)- 

Again differentiating these equations with respect to x and 
y, and writing for simplicity 

J^^A .^=5 ^^-=(7 
dadx ' dbdx ' dcdx ' 

dady ' dbdy ' dcdy ' 
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-we have 

cfxc/y dy dly dy ' 

dxdy dx dx ' dx^ 

t^^ + A-p + Bf+C'^. 

w «y «y ^y 

Hence we find 

\ ay ay , ayj\ ax dx dxj 

\ dx dx dx) \ dy dy dyj 

— ( ATf A*7T\ (^^ ^^ ^ ^^\ 
"" ^ ^ \dy dx dx dy) 

4- (OA'— CA) (— — — ~ — ^ 
^ ^ \dy dx dx dy) * 

Now since a, J, c are connected by two conditions, so that 
h and c are functions of x and y only as being functions of a, 
we have 

da dh ^da^ ^ __ ^^_^ ^ — o 

dy dx dx dy'' ^ dy dx dx dy ' 

dc da ^de_ da _ 
dy dx dx dy~^ ' 
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Thus the above equation reduces to 

(-^y-('-^-)('-f)=». 



or 



df"" ^ dxdy''^ da? ^^* "^"cb? df VdbV"^ ' 

This equation is of the general form (2). Its coefficients 

-j^ , &c, are determinable as functions of a?, y, z^ jp, j when 

the form of the complete primitive (1) is given. For this 
purpose the complete primitive with the two derived equations 
(3) suffice. 

Again, comparing (4) with (2) we have as the conditions 
of their equivalence 

B 8 T __^ V 

dt^ "" dxdy dm? dot? dy* \dxdy) 



dy* dxdy dm? dot? dy* \dxdy/ 

suppose R, 8, Ty U, V 
8'--4:{BT^UV) = 0. 



conditions which suppose J2, 8, T, U, V connected by the 
relation 
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CHAPTER XXX. 



ADDITIONS TO CHAPTER XVII. 



[The present Chapter consists of additions to Chapter xvii. 
Art. 1 was intended to follow Chap. xvii. Art. 1.] 

1. The theory of the solution of linear differential equa- 
tions in a series flows very beautjfuUy from their symbolical 
expression. It is usual«in treating this subject to assume the 
form of the series, and deduce from the differential equation 
the law of its coeflScients ; but the symbolical form of the dif- 
fierential equation determines in reality the form of the solu- 
tion as well as the law of derivation of its successive terms. 

Let us begin with the binomial equation 
Operating on both sides with {/^(Z))}"*, we have 

iJ which *(^)=|(i- 

Hence {1 - ^(2)) €^}t. = {/o(2>)r 0. 

Now {j^(2))]'"^0 will be determined by the solution of a 
linear differential equation with constant coefficients, and will 
be necessarily of the form 

^P+jB(2+CB+..., 

in which .4, jB, C^, ... are arbitrary constants, and P, ^, J8, ... 
are functions of the independent variable. 
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We have then 

therefore u^{l--if>{D)€'^p{AP+JBQ+ GB+ ...). • 

Now let us represent ^{D) e*"^ by p ; then 
u={l^pr{AP+BQ + GB + ...) 
= {l + p + p' + p'+...){AP+BQ+CB+...) 
=^A{l + p + p' + p'+...)P 
+ B{l+p + p' + p'+...)Q 
+ C(l + p + p' + p' + ...)B 
+ ... 

Eepresent the first line of the above expression by w,, 
then since 

p^ = <^(Z>) €'^<f>{D) €*•* ... m times 

= €'^'^ff>{D + mr)(l>{D + mr'-r) <A(i? + r), 

we have 

u^^A{P-]-e'^(l>{D + r)P+e*'^<l>{D + 2r)<P[D + r)P 

+ ^'^<l>{D + Sr)<l>{D + 2r)<l>{D-\-r)P+...}, 

in which it only remains to perform the operations indicated 
by (l}{D + r), by ^{D -{■2r)4>{D + r), ... on the function P. 

Let us in the first place suppose the symbolic ftmction 
f^{D) to be of the form (D-o) {D-b) ...; then 

{f,{D)rO^Ae^' + Be''+.... 

Here P= e"^. Hence substituting in the above expression 
for Uj and observing that/(i>) €^=f{n) €***, we find 

u,=:A^^ {1 + ^(a + r) e**^ + ^ (a + 2r) ^(a + r) 6*'^+ ...}, 

or, since €^ = a?, 

Uj^ = Ax"* {I -f <l>{a + r) a?**+ ^(a + 2r) ^(a + r) 0^+ ...} ; 
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and 

u==Ax'* {l + if>{a-yr) x''+<f>{a+2r) 4>{a + r) aj*'+ ...} 
+ Ba? {1 + if>{b + r) aj'+ (f>{h + 2r)(f>{b + r) iB»^+ ...} 
+ ..., 
the solution sought 

Consider now the general equation 

Here we have, representing "^Ty^ by ^«(i>), 

{l + ^,(2>)6* + .... + ^.(2))6«*}« = {/,(2>)}-0; 
therefore 

«= {1 + <p,io)^+ <A.(i?) e^v {/.(^)r 0. 

Here we have first to determine {y^(i>)}"*0, then to deter- 
mine the effect of the operation represented by 

[1 + 4>,{D)^+ ^.(Z))€»»}-^ 

upon this. 

Now {/o(-D)}*^0 is given by the solution of a linear diffe- 
rential equation with constant coefficients, and will therefore 
be of the form 

AP+BQ+CB+...J 

Ay B, Of ... being arbitrary constants, and P, Q, B, ... func- 
tions of 0. 

Again, since 

{i+^,(i})€*+ +<i>,{D)^r 

1 

"■l+^,(i>)e* + +<f>,{I))e^'' 

it may be shewn by a process of actual symbolical division, 

B.D.E. II. 12 
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attending to the laws of combination of symbols, that ti. 
expression may be expanded in the fonn 

To determine the functions F^{D), F^{D), we maj 

proceed as follows. From the equation 

(1 + ^,(2)) €» + + ^,(Z))^r 

= F,{D)+F,{D)^-{.F^{D)^+ ... 
we have 

l^{l+4>,{Dy+..,+<l,^{D)^][F,{D)+F,{Dy+FlD)^+...] 

= F,{D) + {F,{D) + ,l>,iD)F,{l)-l)}^+ (1). 

Hence F^{D) = 1, 

F,{D) + ,f>,{D)F,{D-l) = 0; 
therefore F^{D) = -4>,{I>)F^{D-1), 

and so on. Hence FJD), F^{D), are determined in suc- 
cession. The general law is as follows: the coefficient of 
c"* in the second member of (1), when m is greater than 1, is 

F^iD) + <l>,{D)F^,iD-l) + <f,,{D)F^{D-2) + ...i2), 

whence 

F„{D)=^-<l,,{D)F^,{l)-l)-4>,iD)F^iD-2)-.... 

By this formula the successive values of F^{D) can be 
deduced from those of F^^{D), F^{D), .... 

Combining the above results we obtain thus for u the ex- 
pression 

u^{l + F,{D)^ + F,{D)e^+....}{AP+BQ + ...} 

=^A{P-]-F,{D)^F+F^{D)e^P+....} 

+ B{Q + F,{D)^Q-{-F,{D)e^Q + ....} 

+ .... 

Let us in applying this expression first suppose that the 
factors oif^{D) are real and unequal, so that fjfi) is of the 
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/ 

'' form {D — a) {D — h) {D — c).... Further, let us suppose that 
no two of the quantities a, J, c, .... differ bj an integer. 

Then {f,{D)Y^(i=A^ + Bi>»+...., 

whence we may assume 

P=e«», <2 = €**,.... 

Thus the expression for u becomes 

+ B{i» + F^{D) 6<"»)»+ J?;(2)) e'»«»+ ....} 

+ ....; 

or, since F{D)^ = F{m)^, 

u =A{^+ F^{a + 1) e'"*''»+i?;(a + 2) e"^''»+ ....} 

+ B{^+ i?;(J + l)e'**'>»+J?;(6 + 2)6'"*|»+ ....} 

+ .... 

Hence, replacing e* bj x, 

u=A{3f + F^{a + l)x'^' + F^(a+2)af^+....] 

+ B{a? + F,{h + 1) oT' + F,{h + 2) ai^ +....} 

+ .... 

In (2) replace in like manner D by a + i and we have, 
putting i for m, 

Fi{a + {) + (f>^{a + i)Fi_/^a + t-l) 

+ ^«(« + **) ^i-2(« + * " 2) + = 0, 

or, if Fi{a + {) be represented by u^j 

Put m for a + 1, thus 

«*«+ ^iH «*«^i + <^, W u^+ = 0. 

This agrees with the law established in [there is no refer- 
ence in the manuscript, but the law intended appears to be 
that given in Chap. xvii. Art. 9.] 

12—2 



Digitized by 



Google 



180 SYMBOLICAL METHODS. [CH. XXli: 

Secondly, suppose that r of the factors oif^[D) are equal 
and of the form B-^a. 

Then {^(i>)}"*0 contains a term of the form 

€^(Co + q^ + c,^+ +Cr.,0- 

Hence the corresponding portion of u is of the form 

where v stands for c^ + cfi + c^ +.... + c^.^ ^"*. 
Now JF;(i) + a + i> 

= k(a + l) + i^/(a + l)i) + i^/'(a + l)^^+....|i;, 

which on performing the differentiations becomes a polyno- 
mial of the form 

-4o + A^+ +-4,.,^\ 

"We see thus that (3) will assume the form of a series of terms 
€"^, e^"'*'*)^ .... each multiplied by a polynomial of the (r — l)"" 
degree in Q. Or arranging the terms otherwise it will con- 
sist of a series of terms of the form 

A+A^+ ^Br_jsr^, 

in which 5^, jB,, B^_^ are series involving e"^, 6^"^^^^, 

e^'^^ Or lastly, changing ^ to a;, the portion of w in 

question is of .the form 

5,+ 5, (log a:) + +5^, aog^rs 

B^y B^,..„B^^ being polynomials in each of which the lowest 
power of X is aj*, and the successive powers increase by unity* 

This establishes the assumption in [there is no reference 
in the manuscript; probably Chap. xvii. Art. 10 is to be 
supplied.] 
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Thirdly, let f^ {D) contain r factors D — a^^ D — a^y.^.D — a^ 
in which a^, a^, ... a^ differ from each other by integers, toge- 
ther with other factors. 

The portion of u corresponding to the factor D^a^ will be 
{1 + jp;(2?) €^ + i^.(i)) €2^ + ....}^€«^, 
in which 

^« (i>) *-• — {^, (i>) e*2^^. (2>) €<-"• 

Thus F^[D) €**^ consists of terms of the form 

i being one of the numbers 1, 2, .... n. Hence i^(2>)€<**'^* 
will consist of terms of the form 

j being one of the numbers 1, 2, .,. n. Continuing this until 
t +y + &+... = w, we see that F^{D) e"^ will ultimately con- 
sist of terms of the form 

^.(2>)6«^,(D)e*^»(2>)e«..., 

i,J, k,... receiving arbitrarily any of the values 1, 2, ... n, 
and i+j+k+ ... being eqnel to m. 

Thus the portion of u derived from A^ will consist of 
all possible terms of the form 

AU^) ^^AD) e*^(i>) ^ .... {^ 
AMD)fAD-t)MD-i-J) ... ^^, 
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Let t = a, t+j=0, i +J + A; + excluding the last 

term =/*; and let the symbolical numerator which involves 
only direct functions be represented by/(i)), and we have 

mm 



f,{D)f,{D-a)f,{D-fi) ...../.(i)-/*) 



g(m+at)fl 



in which a, ^, ....fjb,m are integers ascending by differences 
not exceeding w. 

[A few lines of the manuscript here are obscure, and I 
venture to express in other words the idea which seems to be 
involved. 

Let JD — a, denote one factor oi f^iJD)^ then the correspond- 
ing factors in the denominator of 



are . (-D-a.) (i>-a,-a) (D-a. — fj) (5). 

Now if a, is not greater than a,, then a, + fM is less than 
a^+m; hence no factor in the expression (5) can be identical 
with JD — m — af. But if a, is greater than a,, then one fac- 
tor in the expression (5) may be identical with D — m — a^. 

Hence it follows that the denominator of the expression (4) 
may contain D — m — a^ to the power r — 1, but not to a 
higher power.] 

And, since 

we see that u will contain r sets of terms together of the form 

A + Bi^ogx) + C(loga:)»+ + K(^ogxY'\ 

Ay J?, (7.... being polynomials in x. 

This establishes the rule in [there is no reference in the 
manuscript; probably Chap. XVII. xlrt. 10 is to be supplied.] 
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[There is no hint in the manuscript as to the position 
which Article 2 was intended to occupy ; and the reasoning 
does not seem fully developed.] 

2. Prop. The solution of the equation 

MD)u+f^{I))e'>u+ +f,{D)^u = 

being expressed in the form 

it is not necessary to introduce new constants in interpreting 
J^iiD), .... ; it suflSces to interpret particularly if only uni- 
formly and consistently. 

For let 

{/bC^ro=^P+5(2 + ; 

and in interpreting 

F^{D)e^{AP+BQ+ ) 

let a new constant be introduced which was not in the inter* 
pretation of 

F^,{D)e^^')'{AP+BQ+ ). 

Now F^{D)eJ + <f>,{D)e^F^^{D)e^'-'^' 

+ <l>,{I))e^F^{D)^^^'+ = 0, 

therefore 

F^{D) €^^ = - <I>,{D) €'F^,{D) €(«-^)^- .... 

hence the new constant comes from {^(^)}'^0, and the term 
containing it must be -4T, or B'Q,...., where -4', 5* .... are 
constants. Suppose it A'P; 

then as derived from this, 

i^«.x(i?)e<"^"{/.(i>)}-0 = -^x(i>)e*^'i', 
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Thus F^„{D)^^ M^)^^^..^'''*"'' ' 



The law of derivation is exactly the same as in the deri- 
vation of Fi{D) e^ from Fi_^{D) €(*-^)^ 

[Art. 3 seems intended for a reconstruction on an extended 
scale of part of Chapter xvii. Art. 3.] 

3. We proceed to consider more fully the theory of the 
"binomial equation 

Now the possibility of solving the equation depends upon 
the nature of the symbolic function <I>{D). It is perhaps the 
most general account of the present state of the theory to say 
that there exist certain primary forms of this function which 
render the equation solvable, and that to each of these pri- 
mary forms an infinite number of the forms are reducible by 
general theorems of transformation. As these theorems 
admit of a statement which is independent of the form of the 
function ^(i?), we shall establish them first. 

Pkop. !!• The function <f> {D) in the equation 

can without otherwise changing the first member of that equa- 
tion he lat affected with any constant factor, or 2ndly con- 
verted into <l){D + a),or Srdly converted into {^(— jD)}"*. 

First. Let U=f{e^), and in the equation 
u + <f>{D)e^'u=f{^, 

1 d d 

let ^ = a'^. Then jz = jZ>» ^^^ *^^ equation becomes 
u + a4>{D)e^u=f{a^^, 
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in which i) = -i2j7 . Thus ^{fi) has teen aflfected by a con- 



stant factor a. 

Secondly. In the same system let u — ^v. Then 

or ^v + ^j> {D + a) ^v =/(e^), 

therefore t? + ^ (2? + o) e*^!? = e'^^fi^ . 

Here <f) (Z>) has been changed into ^ (2> + a). 

The result of this transformation may be conveniently ex- 
pressed by the following theorem. 

The equation 
will he converted into 
hy the relations 

Thirdly. In the same equation let ^ = — d' ;, then 

dd^ dff' 
and we have 

w + ^ (- jD) e-'^w =f(€'^, 

in which D^-^. Hence 

u^^'^^{•-^I)^r)]u^f{€^)^ 
therefore e^w + ^ (r — J9) w = e**'/(€"'^, 

whence w + {^(r-2>)r€^w = {<^(r-i))r6*^/(€-^'). 
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In this equation let u = ^v. Then bj the last theorem, 

Thu8 
is converted into 

V + {<!>{- D)re<'v = {<f> (- i?)r/(e-^, 

in which J9 = -^^ , &y assumtng 

The above transformations leave the index r in the first 
member unchanged. If however we assume ^ = - , whence 

^ = a^, we should have 

w + ^(a2))€«'=/(e^). 

By combining this with the previous results we see 
that it is possible to convert 6(i)) into 6(aD + b)f and into 
{<l>{aD + b)r\ 

But the most important transformation of the function <f> (2?) 
is that which is established in the following proposition. 

[The proposition referred to is Prop. iii. of Chap. xvii. 
Art. 3.] 

[Article 4 was intended to follow the words "or subse- 
quently in the derivation of w" in Chap. xvii. Art. 4.] 

4. It becomes therefore important to establish rules for 
the treatment of the constants which in these diflferent ways 
arise. 
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Now the entire process of solution consists of three stages, 
namely : 

1st, the determination of V hj the equation 

2iidl7, the solution of the transformed equation 
Srdly, the determination of u hj the relation 

Let us consider these separately, supposing <f> (D) to con- 
tain a single factor y: — j- which is made to disappear in the 

generation of '^{D), so that a and J differ by a multiple of r. 
Thus the given equation is of the form 

u-^-^{D)€^u=U (6). 

The transformed equation is of the form 

V-^(i))€'^17=F, 

in which u = P,^-^^v, r^Pr^^^U. 

First, suppose a — b = nr, where n is positive. 
Thus 

u = {D + a){D + a-r) ... {D + a-nr+r)v, 

V='l{B+a){I> + a-r)...{D + a-nr + r)Y U. 
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Hence 

= {D + a)...{I) + a-nr + r){l-'^jr[D)i^r{u^+C^<c^ 

+ C,e-''*-^« + + C,e-<''--^»l....(7), 

where l^ is a particalar value of 

|(2)+ o) .... {D + a-nr + r)\ V. 

The part containing the constants will consist of terms of 
the form 

(D + a) ..." {D + a-nr + r) {l-'^(2>)e^}-'(7e-<'^» 
= (2> + a) (D + a-r) ... (D + a-nr-^r) |l +'>lr{D) e'* 

+ 1 (2>) ^^|r{D)^+.. X Ce-<-*f 

= C(i) + o)(i) + a-r)... (i) + o-nr + r)|€-("*->» 

4-^(i)) €-(-*-^»+'^(i>)f (Z) - r)e-C-*-*i»+ ...|. 

Now all these terms vanish up to the one containing 
g-(a-«r)» . therefore we have to perform the operation 

C{D + a){D + a-r) ... (D + o-nr + r) on 

U- iP) ^{D-r)...^{D -jr) €-<-")• 

+ f (2))f (i) - 1) :.. f (i) -jr -r) €-<"^)» + ....k 

where J = n — { — l', that is, we have to perform the operation 
C{D + a){B + a-r)... {D + a-nr + r) on 

U{D) ■<^{D-r)...^{D-jr) e-(-^» 
+ f (D) t»^(2)) ... >^(2> -jr) 6-<«--'5»+ ....1. 
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Now '>^{P)^{P - r) ... ir{D^jr) e'^-^' 

^'\lr{nr — a) yfr^nr^r-- a) ...'^^{nr'^jr — a) e"**""^* 

therefore we obtain 

Thus this expression is the same in form for all values of 
t. Therefore dl the terms containing an arbitrary constant 
in (7) are equivalent to only one term. 

Secondly, suppose a — J == — wr. 

Then w = |(i> + 6) (i) + &-r)...(2> + 5-wr + r)l 'v, 

F= (i) + J) (i) + J - r) ... {D + b-nr + r) U. 

Here there are no constants in F. But u contains n arbitrary 
constants not in v, and as there is no subsequent process in 
the method for destroying these or reducing them to mutual 
dependence, it is necessary that the relations connecting them 
should be sought by comparing the solution with that given 
by the method of development in series. 

Note. It would be better to reduce (6) to the form 

before the demonstration, 

[Article 5 was intended to follow Chap. xvii. Art. 7. 

There is a memoir by Professor Boole on the subject of this 
Article, entitled On the Differential Equations which deter- 
mine the form of the Roots of Algebraic Equations, The 
memoir occupies pages 733 — 755 of the Philosophical Trans- 
actions for 1864.] 
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5. If we agree to regard 9A primary those forms of bino- 
mial equations which are integrable but not through any 
reduction effected by the Propositions of Art. 3, and to which 
equations through the application of those propositions other 
equations are reducible and so made integrable, it becomes 
very important to enquire what these primary integrable 
forms are. It does not appear at present possible to give a 
general answer to this question, but so far as is known, such 
forms if belonging to differential equations of a degree higher 
than the first stand in a remarkable connexion with the 
theory of algebraical equations. By the study of this theory 
Mr Harley was led to the conclusion that y defined as an 
implicit function of x by the algebraical equation 

/-ny + (w~l)a? = (8), 

n being greater than 2, satisfies the binomial differential 
equation 

in which €^ = a?. In this expression the factors of thj nume- 
rator are equidifferent, as of the denominator, their common 

w — 1 
difference being , but the equation is not resolvable by 

Propositions ii. and iii. into forms, the integrability which 
had before been recognised. 

The above result first reached by induction was confirmed 
by Mr Cayley by the aid of Lagrange's theorem. 

To the form (8) all algebraic equations of the third, fourth, 
and fifth degrees are known to be reducible. 

Mr Harley has subsequently found that y considered as a 
function of x defined implicitly by the equation 

y"-wy*"'+ (n-l)a? = 
jsatisfies the sjrmbolical differential equation 
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the factorial notation according to which 

[mY = m{m-l){m-2) ...{m-n + 1) 

being here adopted. 

These results are implicitly involved in a more general 
theorem which I shall now demonstrate. 

Theorem. If J/i^y^ Vn ^^^ *^^ w roots of the algebraic 

equation 

and if the m^ power of any one of these roots he represented 
hy u, and hgaoy 6^ then u as a function of 6 satisfies the dif 
ferential equation 

And the complete integral of ike above differential equation 
will he 

* u^Gj,,^+Cj,,^ + C7^.^ 

Let y** = Zy then the given equation may be expressed in 
the form 

2-1 

in which 6 = - 1. Hence, by Lagrange's theorem, 
the general term being 
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[CH. XXX. 



whicli on effecting the differentiations and adopting the fac- 
torial notation becomes 



m 



1 



-l)r 



-1 



r-1 tw-r 



and this expression will be found to represent the first term 
as well as tne others of Lagrange's expansion provided that 
we interpret the form 



M-'byl. and|>rb7 3^. 



Further, the above general development includes the n 
particular developments of u or y* arising from the giving 



to J* its n particular algebraic values. In this way it 
sents the mr power of each of the n roots y^, y^, .... 



repre- 
• y» in 



succession. 



Now representing the above general term hj u/f^ we shall 
have 



^^ ^ 1 

n 



m 



Wr_» = 



m -f (n — 1) r 

— — - — n 

n 

n [^— n]*" 
Therefore, after reduction and replacing 6 by — 1, 






n \ [n n"^ J 



.(9). 



It follows therefore that the complete series of which the 
general term is tt,a*' will if represented by u satisfy the diffe- 
rential equation 
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^ n n \ \n n J 



[nf' -^u = o...{i). 

If we integrate the equation in a series (Cliap. xvii. Art. 9), 
ihe initial terms of the value of u will be 

C,+ C,a+C,a' + C^,ar-\ 

the succeeding terms being formed from these by the law 

(9). Hence, if the arbitrary constants C^C^ <7«-.i be so 

determined as to make the above initial terms agree with 
the first n terms of the Lagrangean expansion in any of its 
particular forms, the succeeding terms will also agree, and 
the Lagrangean expansion will thus become a particular inte- 
gral of the equation (I). The aggregate of such particular 
integrals, each afiected by an arbitrary constant, will therefore 
also be an integral of the differential etjuation, and will, in 
fact, constitute its general integral, subject to exception only 
in the case in which for a. particular value of m the integrals 
2/ 1^9 yr^ y»*" cease to be independent. 

For instance, if 7n= — 1, and we reduce the equation to 
the form 

it is seen that except when ti = 2, we have 

yr'+yr' +y«-^=o. 

Here then the solution 

u = 6V+ Cj,,^ + (7«y/ (10) 

ceases to be general for it becomes • 

« = (C,-Cf.)y-+((7.- a„)2,-... + (CL.- C,)y-^^^, 
and virtually involves but n — 1 arbitrary constants. ^ 
If, however, we give to the integral the form 

B.D.E. IL 13 
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the last term of which becomes a vanishing fraction when 
?n = — 1, we find for the general value of u in this case 

^ = Oj/^-' + C^," + CLi^^-i 

+ dnij/r^^ogy^ + yf'logy^ +y»"'logy^), 

and in this way we may proceed in failing cases generally. 

Lastly, it may be observed that in certain cases the diflferen- 
tial equation (I) admits of reduction to an order lower by 
unity than its own. And in particular this happens in the 
failing cases above noticed. Thus, if in (I) we make wi = — 1 
the equation will be expressible in the form 

D(D-l)...(D-w + l)w 

n\^ n w J ^ ' ' 

whence, operating on both members with (i) — w 4- 1)~*, we 
have 

^ ^ n\_ n ^ J 

The general integral of this equation will be expressed by 
(10) provided that a proper relation be established between C 

and the constants C^, C^, C^. If we choose to determine 

C so as to give to the integral the particular form y~\ we shall 
find on substituting for u its Lagrangean development making 
m = — 1, J = — 1, and calculating the coefficient of a""^ or 
^n-i)fl ijj ^i^Q fij.gt member of the difierential, 



a. 



[n - 3]" 



Hence, if n be greater than 2, we have (7=0. It follows 
therefore that if n be greater than 2, the equation 

^^r'u-l\^i>-l-iY^u=o (ID, 
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in which e* =a has for its general integral 

y^, 2^2, y^-i being any n — 1 roots of the equation 

It may he useful to notice the forms which the ahove 
results assume when is changed into — 0, and therefore D 
into — jD ; see Art. 3. 

It will be found that (I) becomes 

L 7i nj \n nj 

of which the integral is therefore 

^ = c^i^r + 0,y^ .....+ a,y,« 

2/1,^2^ 1/n being the roots of the equation 

y-_l2,-i + l=0 (11); 

and log a being denoted by 0; 

while as the equivalent of (II) we have 

r n _ 9"|«-i 

""^ Tn-l iY^ ^^ = <> (IV), 

[_ n n_ 

of which, supposing n greater than 2, the integral is 

j/iyj/gj j/n-i being any w — 1 roots of the same algebraic 

equation. 

13—2 
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Mr Harley's results may^ readily be deduced from the 
above. Thus it will be found that the equation (11) re- 
duces to 

r-n< + (n-l)a; = 

if we make 



1 1 



y-»=(«-l)-a,-«, a= ^"-^j" '^ 



Hence, making x=^ and representing ■^, by D\ we have 



for the transformation of (IV) 



n ' 



2) ^2>', 

W— 1 



1 -f 



t^ = (n — 1) " € * t. 
Substituting and multiplying the result by e-, we find 



r n J., 2n-r 
^ /n-l\*^ M^-l w-1 

'-(-irj ^ — TW — 



.e^^')H=Oy 



which is Mr Harley*s first equation. 



m 



If in (I) and (III) we make 1 = a, whence m^^n-^nay 



n 



and at the same time change a into ah ", and y into yJ *, we 
shall obtain the following somewhat more general statement 
of their united import. 
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The differential equations 
Tn-l 



1 
u 



- — -^ — '-'""'• 



b 



._i i^nir 



L=^^— r(f— ) 



.6^tt = 0, 



are both satisfied by the general integral 

; ^ c,yrr+ G^^-^ .-.. + o^yr^^, 

when y^, y,, .... y^ are the roots of the algebraic equation 

y*-ay*^ + J = 0, 

provided that for the first equation a = €*, and for the second 
« = €-*. 

' If n<=: 2, the above equations assume the forms 

"-& — ^w^ — "^^^' 

f-(a-l)' 

6. [The two principal papers by Mr Harley on the dif- 
ferential equations exhibitea on page 190 are the following : 



(1) On the Theory of the Transcendental Solution of 
gebraic Equations, f 
Vol. V. pages 337... 360. 



Algebraic Equations, Qtuirterly Journal of Mathematics^ 
"ol. 



(2) On a certain class of Linear' Diflferential Equations. 
Manchester Memoirs.. Third Series. Vol. ii. pages 232... 245. 
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In a letter bearing date January 13, 1864, Professor Boole 
pointed out to Mr Biarley that his second equation might also 
oe deduced from the general theorem discussed in Art. 5. 
Employing the above notation the deduction may be pre- 
sented in the following form. 

The equation (11) will reduce to 

r-n^ + (n-l)a; = 0, 
if we make 



and for the transformation of (III) we have 

1 

>-^I 



^^Un^lfe, B^nU, 



t« = (n — 1) " € • u. 

These substitutions being effected we arrive, after some 
slight reductions, at the following equation, 

n"[(n-l)2y-m]'*-^i)V-(n-l)[ni)'-m-l]VM'=0, 

which, making w = 1 and u = <, gives 

n*[(n-l)D'-l]*-'I>'«-(7i-l)[niy-2]V« = 0, 

an equation which admits of reduction. In fact, operating on 
both members with {1/ — 1)"^, and determining the constant, 
as in the former case, by the aid of the Lagrangean expansion, 
we find' 

n^\(ji-\)VYH- [niy - n - 1) [nU-^il^eU = [n- 1]»-^€^, 

which is Mr Harley's second equation. 

The references and deduction here given were to have 
been added to the memoir which is cited in page 189, ac- 
cording to Professor Boole's desire; but by some accident 
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they were not printed, and the omission was not discovered 
until after his death. 

Mr Harley has lately succeeded in obtaining the fol- 
lowing extension of Professor Boole's theorem. 

The differential equation 



«' [^^]"« 



Vn-r d ^ m .1 [r d m SY r a 
\_ n ax n \ \n ax n J 

is satisfied by the wi*^ power of any root of the equation 

^* — a?^'*"'' + a = 0, 

u being considered as a function of a?. 

From this he deduces the following ; the differential equa- 
tion 

L r dx r\ L ^J 

is satisfied by the m*^ power of any root of the equation 
y*-w^*"'*'-f (n- 1) aj = 0. 

For the materials of this Article I am indebted to Mr 
Harley.] 
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CHAPTER XXXI. 



THE JACOBIAN THEORY OP THE LAST MULTIPLIER. 



1. A SYSTEM of n differential equations of the first order 
and degree containing w + 1 variables admits of n integrals 
of the form 

W^ = Cj, Wj = Cj, ••• W„ = Cn^ 

Wj, Wj, ... w„ being independent functions of the original vari- 
ables. When w — 1 of these integrals have been found they 
enable us to eliminate n — 1 variables, with their differentials, 
from the given system of equations, and so to obtain a single 
final differential equation of the first order between the two 
remaining variables. The final equation admits of being 
made integrable by a factor, and its solution so found would 
constitute the n^^ and last integral of the system. We pro- 
pose in this Chapter to develope the theory of the above 
integrating factor as established by Jacobi. The term * prin- 
ciple of the last multiplier,' which is more usually employed, 
seems objectionable; for the essence of Jacobi's discovery 
consisted not in demonstrating the existence or the nature of 
the last integrating factor, but in the peculiar form of the 
method which he gave for its determination, and in the rela- 
tions which are implied in that form. The discovery may be 
briefly said to consist in this; viz. that instead of forming by 
means of the n — 1 known integrals the final differential 
equation between two variables and applying methods analo- 
gous to those of Chap, v., to determine its integrating factor, 
we construct antecedently to all integration a linear partial dif- 
ferential equation of the first order, any one integral of which 
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will enable us to assign an integrating factor of the final difi*er- 
ential equation, whatever the order of the previotcs integrations 
may have been. Again, this partial differential equation de- 
pending for its construction only upon the form of the system 
given, we can often by examining it affirm beforehand that if 
all the integrals but one of the system be in any way found, 
the final integral will be deducible by quadratures. This 
happens in the case of the most important of all systems of 
differential equations — that of Dynamics. 

Further, an ordinary differential, equation of the n^ order 
being reducible to a system of n differential equations of the 
first order, Jacobi's theory may here also enable us to pre- 
dicate the possibility of the last integration when the previous 
integrations have been effected. 

Beginning with a single differential equation of the first 
order reduced to the form 

dx ^dy 

in which X and Y are functions of the two variables x and y, 
we know by Chap. v. that the inte^ating factor fi will be 
given by the solution of the partial differential equation 

^^T-o «• 

the form of which should be carefully noticed. 

Consider next a system of two differential equations of the 
first order expressed in the general form 

dx ^dy ^dz y^v 

X""T""Z ^ ^' 

X, F, and Z being functions of the three variables oj, y, z^ and 
suppose one integral, represented by 

4>{x,y,z)=^c .....(3), 

to be known. The function <f) {x, y, z)y or, as we shall express 
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it for brevity, <p, will obviously satisfy the partial differential 
equation 

^f+^f+^s-" (*)■ 

of which indeed the given equations form the Lagrangean 
auxiliary system ; see Chap. xiv. 

If from the given integral we determine 2? as a function of 
X, y and c, and substitute its value in the first of the given 
dilierential equations, viz. 

dx ^dy 

the latter will be converted into a differential equation be- 
tween X and y. But we may leave to the equation its prior 
form, provided that we regard X and Y as functions of the 
variables x and y, both explicitly as they appear therein, and 
implicitly as they are involved in z. And this being so, the 
equation (1) will become 

djfiX) djfiX) dz d(fiY) difiY) dz^ 
dx dz dx dy dz dy 

The values of -j- and -j- in this equation must be found 



from the known integral (3) ; they are 

dz d(f> d(j> dz d^ d<f> 
dx^ dx ' dz^ dy"^ dy * dz^ 

substituting which we have 

d{tiX) d^ d(jjLX)d(f> d(fiY) d4> d{fiY)d^^ .^. 
dx dz dz dx dy dz dz dy " 

This then is the partial differential equation for determin- 
ing fjL. But the construction of this equation supposes 6 to 
be known. We propose to shew that fi can be determined by 
a process in which the only partial differential equation to be 
solved can be constructed without the knowledge of if>. 
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Since by actual differentiation 

dx \ dz) dz \ dx) " dx dz dz dx* 

it follows, writing fjbXtov A, that 

d (fiX) d^ d {fjbX) d<f>_ d / Y^\^ — f X^^ 
dx dz dz dx^ dx\ dz) dz \ dx) ' 

Similarly 

d{fiY) d<f> djfiY) d^^ d / y^^—f Y^\ 
dy dz dz dy dy V d^J dz \^ dy) * 

Lastly, we have 

Now adding the last three equations together we see that 
the first member of the result vanishes by (5) : we have thus 

-s('^S)-l(''^f)-l('^t)=«- 

The second line of the first member is equal to 
and therefore vanishes by (4). There remains then 
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it 4lence if we put 



"S-^. 



we have 



d(MX) ^djMY) ^djMZ)^^ ,g, 

dx dy dz 



If then by the solution of this equation a value of M dis^ 
tinctfrom be founds the function •r? will be an integrating 

Tz 
factor of that final differential equation which remains when z 
has been eliniincUed from the system (2) by means of any 
knovm integral ^ = c. 

It will be observed that the equation for M is analogous in 
form to the equation for fi in the previous system. And this 
suggests the form of the general theorem. 

Thus proceeding to the case of a system of three equations 
dx ^dy ^dz __dt 

we see that if 

be a known integral, -^ therefore satisfying the equation 

'^s-^f+^t^^t=» «■ 

then the system 

dx ^dy __^ 

will virtually, involve only the variable^ a?, y, «, since t 
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through the known integral becomes a function of x, y, z. 'x^ 
equation (6) now becomes 

d{MX) djMX) dt djMY) djMY) dt 
dx dt dx dy dt dy 



or putting 



d{MZ) djMZ) dt_ 
^ dz '^ dt dz^^' 



dt _ d'^ ^ dyfr 
dx'~ dx ' dt '"' 



djMX) dyfr d{MX) dyjr d{MY) dyjr d{MY) dyjr 
dx dt dt dx dy dt dt dy 

djMY) d'^ d(MY) d^|r^^ 
dz dt dt dz ^ 

and this is equivalent to 

dx dy dz dt 

and therefore becomes on rejecting the term in the second 
line by (7), and putting 

d{Nr^ d{NY) d{NZ) d(NT) _ 
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If from ttis eq uation a value of N distinct from be ob- 

N 
tained, then Jf=^, and therefore 

N 
^ df d^ 
dt dz 

This is the final multiplier, i.e. the integrating factor of 
the final differential equation between x and y which remains 
when z and t have been eliminated from the given system by 
means of the two known integrals. In calculating /l6 from 
the above formula we must proceed as follows. The value of 

-^ must be found from any given integral -^ = c ; but that of 

•=^ must be found from another integral from which by means 

of the former one t has been eliminated. Thus the general 
forms of the integrals will be 

^(a:,y, z,c)=c\ 

Lastly, the values of -^ > -^ found as above, and that of 

N given by any solution (distinct from 0) of the partial dif- 
ferential equation (8) having been substituted in the expres- 
sion for /i, we must eliminate z and t from that expression by 
means of the two known integrals. The resulting function 
of a?, y, c and c' will be the integrating factor sought. 

The reasoning above employed is in its nature quite inde- 
pendent of the number of tlie equations of the original sys- 
tem. The general theorem to which it leads may be thus 
stated. 

Theorem. The system of n differential equations 

:^ — ^1 — ^2 — ^ 

X r;-r,--i; 
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being given, if a system of w — 1 integrals 

"be so reduced hj elimination that the variable y^ shall not 
appear in ^^, the variables y^, y^ shall not appear in ^3, and 
so on, then the integrating factor fi of that final differential 
equation between x and y^ will be given by the formula 

M 

in which M represents any integral distinct from of the par- 
tial differential equation 

d[MT) ^ djMY:) ^ ^^ djMY: )^^ 
dx dy^ **" dy^ 

In applying this theorem the expression for fi must be 
freed from all the variables except x and y^ by means of the 
given integrals. 

This is Jacobi's theorem. On account of its great importance 
I propose to give another demonstration of it founded upon 
the Calculus of Variations. 

2. Second demonstration founded upon the Calculus of 
Variations, 

It will be most convenient to present the proposed system 
of differential equations under the symmetrical form 

dx^ dx^ _ dx^ 

the independent variables being x^, x^, a?„ of which 

Xj, Xa, .... X„ are any functions. We have thus w — 1 diffe- 
rential equations, and we are to seek the integrating factor of 
the differential equation which remains when by means of 
w — 2 known integrals n — 2 of the variables with their diffe- 
rentials have been eliminated. 
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Suppose F—c to be any integral of the system, then P 
satisfies, and it suffices that it satisfies, the partial differential 
equation 



y dP ^^dP 



.+^??-«. 



dx, • ~'dx/'"^^dx,~^ (''^• 



Now if in place of a?„ »,, ....oj, we introduce a new sys- 
tem of independent variables u^, t^, .... u, which are functions 
of the former, then we shall have 

dP dP ,y-dP 



£/■,, Ut, .... U^ being functions of «,, «, m„. And by the 

theory of the transformation of multiple integrals, . 






where 



^= 



du^ 



du^ 



dx^'' 



du^ 
dx^ 



The foregoing equation we may express in the form 
^ I -2i~T- dx^ dx^ ... da?^=S / -jf T~^^i ^^« •••• ^*^«i« 

Hence, representing by S an operation of differentiation 
which affects only the form of P as a function of a?j, ar,, ... xr^ 
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or of Wj, Wa, ....a^, and not the independent variables them- 
selves, we have 

and therefore integrating by parts and equating the portions 
on each side which remain under the sign of n-fotd inte- 
gration, 

Whence again transforming the integral in the first member 

^ p dXi ^p du^ du^ .... du^ 
J dxi H 

and this being true quite irrespectively of the form of P, 
we have 



H^ dx,'^^ du\H) 



In this equation Jacobi's theorem is virtually contained. 
For let the given equation be multiplied by any factor. Then 
changing in the above Xi into MX^y and U^ into JfC^, we 
have 

1 ^ d{Mx:f _^ i.(^\ 

H^ dx, ^^duXs)* 
Hence, if If be determined to satisfy the equation 

aXi 
B.D.E. IL 14 
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we shall have 

^i(^-» (*)■ 

This is wholly independent of the relations connecting 
Wj, w,,...w„ with a-j, a?,,....a;„. Now choose the n — 2 variables 
t*j, Wj, .... w^ so that Wi = Cj, Wa = Cj, .... w,^ = c,^ shall be 
integrals of the given partial differential equation (a). Then 
that equation transformed becomes 

dP dp_ 

of which the auxiliaiy ordinary equation is 

At the same time the equation (b) becomes 
d fM „ \ . d IM . 



d^n-t 



(f^.^)-4(l^.)=»- 



M 
Hence j> is the integrating factor of the preceding diffe- 
rential equation between u^^ and u^. 

Jacobi's theorem in ita most general form is thus seen to 
be the following 

Theorem. If the system of differential equations 

be transformed by the introduction of a new system of vari- 
ables Mj, w,, ....w„, so chosen that 

t/j = Ci, t^j = c„ Uf^-C^ 

shall be integrals of the given system, then the final differen- 
tial equation between u^^ and u^ shall have for its integrating 
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M 



factor -g:, in which M is any function satisfying the partial 
differential equation 

djMX,) djMX;) ^ d{MX:) _^ 



dx. 



dx^ 



dx^ 



and J? stands for the determinant 

du. du^ 



dx,''- 


-dx„ 


dx^'"' 


'"dx. 



The form of Jacobi's theorem obtained by the previous 
demonstration may be deduced from the above by choosing 
for u,^^u^ two of the original variables, for example a?,^, a?„, 
and transforming the integrals m^, w^, .... w,^ so that w, shall 
contain only a?,, ... a?„, W3 shall contain only a?3...a?„, and so on. 



Examples. 

3. Jacobi has established by means of the above theorem 
the very remarkable theorem that in any ordinary dynamical 
problem the forces depending not upon the time but upon the 
material constitution of the system, if all the integrals but 
two of the dynamical equations are found, the two remaining 
integrals can be found by quadratures. 

1st. In a dynamical system of free points the forces act- 
ing upon which depend only upon the position of the points, 
we have if we represent the entire system of rectangular co- 
ordinates taken in any order by aj, y, «,... and the correspond- 
ing resolved forces divided each by the corresponding mass 
by Xj Yj Z^,,, the system of equations 



de 



= X, 



de 



— Y 



14-2 
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or putting 

dx ^ f dy ^ f 



X y X Y 

Now as X, Y... do not contain t we may consider first the 
system 

dx _^dy _ dx' _ dy 
X y XI 

and it is evident that if we can find all the integrals of this 
system, t will be given by the equation 

fdx 



-/$+. 



x' having been first converted by means of the supposed in- 
tegrals into a function of x. 

To determine the last multiplier of the system last writt^ 
we have first the equation 

d{Mx') d{My') d{MX) , d{MY) 

dx '^ dy ••••^ dixf "^ dy' "••""' 

which since X, F... do not contain x\y\.. is satisfied by 
Jf = a constant Giving to the constant the particular value 1, 
we see that if 

are w — 2 integrals of the system, and if by means of these 
we eliminate w — 2 of the variables and construct the differen- 
tial equation between the two remaining variables, the inte- 
grating factor of that equation will be ^, in which H is the 
functional determinant of W|, Wg, .... w«v 
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2ndly. Suppose the system subject to a material connex- 
ion which establishes an equation of condition among some 
or all of the co-ordinates. If we represent the co-ordinates 
taken in any order and multiplied each by the square root of 



the corresponding mass by a;, y, . . . the corresponding resolved 
forces by X, F, ... and the equation of condition expressed by 
means of the above modified co-ordinates by ^ = 0, the diffe- 
rential equations will be 

the transformation above employed reducing all the equations 
to the same type. [See the next Chapter.] 

Making 

the system becomes 

j^ dx dy dx dy' 

dx ay 

and the Jacobian equation for M becomes 

dx dy "" dx' 

dy 

Now <^ does not contain a?', y'.... Let us inquire whether 
it is possible to determine M also as a function of a?, y .... 
without x\ y.... so as to satisfy the above differential equa- 
tion. 
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The equation woiild become 

, dM . , dM , , ,, fd\ dd> , dk d6 , \ ^ 

or if we write 

, d . t d ^ 

and from this we must eliminate X. 

Now since ^ = 0, we have by differentiating and putting 
dx _^ , 
dt^"^'"' 

^'^ + 2^1+- = ^' 
and again differentiating 

Mda^d^ dy _ 
"^dx dt^dy di"^ "' 



or since 



'df-'^'^^d^''- 
* da?^^ df^""^^'"^dxdy^"" 

+x^+r^ + .... 

ax ay 

-{©■HIT--}-. 
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and differentiating with respect to a?', 



or if we make 



Similarly 



Therefore 






=*l^l «=«• 



or 



dx dx dy dy 

fd\ d(f> d\ d^ \ n ^ f\ 

and now elimmatmg ^, £ + ^^+ ... 

we obtain QSif - 3/8^ = 0, 

which is satisfied by if = Q. 

4. [Among Professor Boole's manuscripts I found five 
pages in German, forming part of a memoir, which was pro- 
bably intended for Crelles Mathematical Journal, The 
memoir was to have discussed two applications of the Calcu- 
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]u8 of Variations; one to the Jacobian Theory of the Last 
Multiplier, and the other to the Solution of Pfaff's equation 

X^dxi + X^dx^+ +X^dx^=^0, 

But there is only a single paragraph relating to the second 
application. 

The manuscript contains the same demonstration of the 
Jacobian Theory of the Last Multiplier as in Art. 2 of the 
present Chapter; after this demonstration some remarks occur 
of which the substance will now be given.] 

It is worthy of notice, that Jacobi in the 36th volume of 
Crelle's Journal, deduced by the aid of the Calculus of Varia- 
tions the result on which the preceding demonstration of the 
Theory of the Last Multiplier depends. In fact, he shewed 
that if V denotes any function of 



^«i ^j 



dz dz 



and Fbe transformed by the introduction of a new*system of 
independent variables Wj, w,, ... u^, then the following rela- 
tion holds. 



where 



^/dV d dV 
I dz dx. J dz 


d dV\ 
dx,^dz\ 


d{LV) d d{LV) 
dz du, jdz 
du^ 


d d{LV) 
du, T dz 




dx. 


"du. 




A — 










dx. 


dx. 






du,' 


du^ 
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JacobI applies this result to the transformation of the ex- 
pression 

d'V ^F ^ 
dx^ dy^ dz^ * 

But neither Jacobi himself, noi^ any other person, so far as I 
know, has drawn attention to the application of the result 
which I have given here. 

[The substance of the single paragraph relating to the 
second application of the Calculus of Variations will now be 
given.] 

Clebsch has earned the thanks of all who are interested In 
the higher parts of the Theory of Differential Equations, since 
he has performed the same service for Pfaff s problem as 
Jacobi did for the Theory of Partial Differential Equations of 
the first order, and thereby for the equations of Dynamics. 
But while I recognise the great importance of the results, I 
consider it desirable to give a simpler deduction of the system 
of partial differential equations therein involved, and on which 
the other results depend. 
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CHAPTER XXXII. 

THE DIFFERENTIAL EQUATIONS OP DYNAMICS. 



[It will be seen that this is only a fragment of the Chapter 
which was to have appeared under this title.] 

I do not propose in this Chapter to discuss the origin and 
interpretation of the differential equations of motion or to enter 
into those details of their application which are found in all or- 
dinary treatises on Dynamics. But they constitute a system 
analytically so remarkable from the forms in which it is 
capable of being expressed, and from the general methods of 
integration which emerge out of those forms, that they are 
well deserving of a special attention. 

Eeferred to rectangular co-ordinates the differential equa- 
tions for the motion of a system of points free or connected 
are 
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Here m is the mass at the point (a?, y, z), m! that at (x\ y\ «'), 
X, F, Z the resolved forces at (ar, y, z) tending, severally to 
increase those co-ordinates, and so on. Lastly 

«^=0, i/r = 0,... 

are the equations of condition each of which may involve all 
the co-ordinates, and X, /x... are indeterminate multipliers. 

The above is usually termed the first Lagrangean form of the 
difibrential equations. In applying it we must either elimi- 
nate \, /i... from the given equations, and then by the equa- 
tions of condition just so many of the co-ordinates with their 
differentials, or we must retain X,/i,. . . as variables so conditioned 

that the values of -^ , -^ ,. . . in the system shall satisfy iden- 
tically the differential equations involving -^, -^j-** de- 
rived from <^ = 0, 'i|r = 0,... viz. the equations 

^ = ^ = 

de^' df "'••• 

The first Lagrangean sjrstem may by a slight transforma- 
tion be reduced to a form in which all the equations are of 
one type, viz. of the type which they would have if all the 
masses were equal to unity. 

For taking the first equation of the system and dividing 
by m^ we may express the result in the form 






from which we see that if a:, y ... had been taken to represent 
the entire system of co-ordinates taken in any order and mul- 
tiplied each by the square root of the corresponding mass, and 
X, F... the corresponding resolved forces taken in the same 
order and divided each by the square root of the correspond- 
ing mass, the system of equations would have been 
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df dx <?a3 **' 

df^ dy '^ dy '" 



all being of one type. In general investigations this form 
is to be preferred. 

From the first Lagrangean form another known as the 
second Lagrangean, and from this again a, third known as 
the Hamiltonian are derived. The second Lagrangean form 
is properly speaking an expression for the effect of a trans- 
formation of co-ordinates in the most general sense upon the 
original system, i.e. of a transformation which in place of 
a;, y, ... the entire system of given co-ordinates substitutes 
a new system of variables f, i;, ... the expressions of which as 
functions of a;, y, ... are known. It is not necessary that this 
new system of variables should be co-ordinates in the proper 
sense of that term, determining three by three the positions 
of the several masses; it suffices that they should in their en- 
tirety determine and be determined by the co-ordinates given. 

The second Lagrangean form may be established as 
follows : 

Differentiating the equations <^ = 0, ^Irrs 0,... with respect to 
any one of the new variables f we have 



d<f) dx d(f> dy _ 
dx d^ dy d^ '" "" ' 

dyjrdxdjrdy ^ 
dx d^ dy d^'" ' 



whence if we multiply the equations of the given system by 

dx dy 

di'dl' 



dx dy 111 1 

jfc , -ii ,... and add, we nave 



dx d^x dy d^y __ ^.dx ydy 
dlW^d^W ^dl"^ ^ dl" 
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CHAPTER XXXIIL 

ON THE PROJECTION OF A SURFACE ON A PLANE. 



[The following memoir was found among Professor Boole's 
manuscripts; a Title and Introductory Remarks were to have 
been prefixed, but with this exception the memoir appears to 
be finished for publication. It is sufficiently connected with 
the subject of Differential Equations to find a place in the 
present volume. 

The memoir by Sir John Herschel to which allusion is 
made is entitled, On a new Projection of the Sphere; this was 
read before the Royal Geographical Society of London on 
the 11th of April, 1859, and was printed as part of the Journal 
of the Society, Vol. xxx. 1860, pages 100... 106. A chart of 
the World on Sir John flerschel's projection has been pub- 
lished by A. and C. Black of Edinburgh. 

The history of the subject will be found in Chapter xxiil. 
of the Coup d'oeil historique sur la Projection des Cartes de 
OSographie... Par M. D'Avezac, Paris, 1863. 

For the materials of this introductory notice I am indebted 
to Sir John Herschel.] 

1. Let Xy y, e be the rectangular co-ordinates of any point 
on the given surface; x', y* the co-ordinates of the correspond- 
ing point on the plane of projection. Let the equation of the 
. given surface be 

or, for simplicity, 
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The condition of projection upon which Sir John Herschel's 
investigations are founded, and which we shall adopt here, is 
that of the similarity of corresponding infinitesimal areas on 
the surface and on the plane. The object of the problem then 
in general is the discovery of the mode in which a?', y' depend 
upon a?, y, and z in accordance with the above condition; its 
object in any particular case is the determination of a?', y' as 
functions of a?, y, z. 

Regarding then a?', y' as ultimately functions of a?, y^ z we 
have 

, , doti , . doi ^ . da?' , 
cto = -7- oa; + —r- ay-\'—rdZi 
ax dy ^ dz 

in which dxy dy^ dz ^^e not independent, but are connected 
by the condition 

dF , dF , . dF , 
-j-tto + — - ay +-T- efo == 0. 
dx dy ^ dz 

Now for brevity write 

doi _ dai ^ - cfo' _ 
dx "^ ^ dy ^ ^ dz " ^ 

— =u4 — =jB — =(7- 
rfa? ' fl?y ' dz ' 

then 

&j' = ae&H-Jdy + c<fe (1), 

dy :=adx + b'dy-\-ddz (2), 

O^Adx-hBdy-^- Cdz (3). 
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Now the condition of the similarity of infinitesimal cor- 
responding areas may be resolved into the two following 
conditions, viz.: 

1st. The equality of their corresponding angles. 

2ndly. The proportionality of their corresponding sides. 

And these conditions we shall introduce separately. 

1st. Assuming any point a?', y on the plane of projection, 
let oi alone vary, and the infinitesimal line generated is daj', 
while (since dy' = 0) (2) and (3) become 

ddx 4- Vdy + ddz = 0, 

Adx + Bdy-\- Cdz = 0, 

whence, if we write 

L^Bc-CV, M=^Ca'^Ac\ N^AV-'Ba;, 

_ dx dy dz , . 

we have T^M^N ^^^' 

so that the direction cosines of the infinitesimal line on the 
surface F corresponding to the line dx' on the plane (a;', y') 
will be 

L M N 

(L»H-if+i^)*' (X^ + if'+i^)*' (i' + Jf^+i^)* ^ ^' 

In like manner, if y' alone vary, we shall find for the 
direction cosines of the infinitesimal line on the surface F 
which corresponds to dy' on the plane 

r ^ N' 

{L'*+M"+N'*)^' (i'^+Jkr^+i^)*' {L'*+M'*+N"'}i ^ ''' 

T^here L'^Bc- Cb. M^ Ca-Ae, N' = Ah -Ba. 

Bj the first of .the conditions of similarity the (ingle be- 
tween these lines on the surface z^nst he a right angle since 
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dai and dy' are at right angles. Hence we have, from 
(5) and (6), 

LL + MM' + NN'=^0 (7). 

2ndly. The ratio of the length of the element dx' to the 
corresponding element on the surface is 

dx' 



Jdx'-^dy' + dz^' 
or, by (1), 

adx + hdy 4- cdz 

4d3?^-dy'^dz' ' 
and therefore by (4) 

equating which to the corresponding expression for the ratio 
of the length of drj to that of its projection on the surface, 
we have 

^TJTW+W ^/IrTW^+lr ^ ^' 

Now if we substitute for X, M, Nj L\ M\ N' their values, 
we shall find 

aL + bM+cN=^A{Vc-hc')+B{c'a-'Ca') + C{a'b - ah% 

a'r+VM'+ c'N' = A {he'- h'c) + B{ca'-' c'a) + G {aV-a'h), 

and the second members of these equations differ only in sign. 

Thus (8) may be expressed in the form 

iA{Vc-ld) +B{c'a^ca') + a(ai'-a'J)l 

\{L' + M' + N')^ {L" + M"+N")V '"^^' 
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But the first factor of the first member of this equation 
being the determinant of the system 

adx + hdy + cdz = 0, 

Adx^-Bdy + Cdz = 0, 

expresses when equated to zero the condition that if in the 
system (1), (2), (3) dy* vanishes dai shall also vanish; and 
doi and dy being independent, this condition cannot be satis- 
fied, so that (9) reduces to 



[ = 0, 



whence 

i'»+jir'+JV^-.z*-2r-.^«o (10), 

and this, with (7), will fiolly express^ the conditions of simi- 
larity, 

2. If we multiply (7) by 2 V^, and add and subtract 
the result from (10), we obtain the equivalent system 



!i' + i V- 1)« + (if' + if V:^i)» + (J^ + -N:V-i)* = oi 
!;'- i V^/ + (if '- if V=l)« + (JV- JvV=i)« = 0^ ' 



(11). 



-vr Tt ^ T t — ; dW doa dF dx* 

dff dt^ dz dy 



± 



\S^ dz. dz -dy) 



dF d{x'±^W'-l) dF d{x'±y'^Pl) 
dy dz dz dy ' 

B.D. E. II. 15 
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Writing then 



[CH. XXXIII. 



we have 



a?' +y' V- 1 = u, a?'- yV- 1 = 1;, 

dy dz dz dy^ 
dy dz dz dy* 



In the same way 



dz dx dx dz^ 

dz dx dx dz^ 

dx dy dy dx'* 

dx dy dy dx^ 



Substituting which in the system (11) there result 

\dy dz dz dy) '^ \dz dx dx dz) 

^ /dF du_dF du\\ 
\dx dy dy dx) " 

(dF dv _dF do"^ /dF d^ yjdv^ 
\dy dz dz dy) \dz dx dx S/ 

\dx dy dy dx) " 



.(12), 
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to wkich we maj give the somewliat more convenient form 
\dx dx dy dy dz dz) ^ '* 

{(S)'-(f)^(f)}{©'-(|T-(S)} 

\dx dx dy dy dz dz) "~ '"^ 

These are partial differential equations of the first order, 
serving to determine u and v as functions of a;, y^ z. 

But it is not necessary to solve the equations in their 
general form. For, a?, y, and z being connected \yy the equa- 
tion of the surface, the above equations msj always be so 
reduced as to involve only two independent variables. As 
latitude and longitude determine the position of a point on 
the earth, so two co-ordinates of any given species will deter- 
mine the position of a point on the given surface, and these 
co-ordinates, when fixea upon, become the independent varia- 
bles of the problem* 

Let 8 and t represent such co-ordinates, and let their ex- 
pressions in terms of a;, y, z give 

which equations combined with that of the given surface will 
reciprocally determine a:, y, z as functions of 8 and t Then 
Ist the differential coefficients of F which in the equations 
(I), (II), are functions of a?, y, z may be transformed into func- 
tions of 8 and t ; 2ndly , we have 

du ^du ds du dl 

dx^ ds dx dt dx^ 

du __du ds du dt 

dy" ds dy dt dy'* 

du _ du ds du dt 

dz" ds dz dt dz ' 
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and as -T- , t"*** *^® known functions of a?, y\ «, they also 

are expressible in terms of 8 and t. The result of these sub- 
stitutions will then be to convert (I) into a partial differential 
equation in which u is the dependent and a and t the inde- 
pendent variables, and this equation being, like (I), of the 
first order and second degree in the differential coefficients of 
w, will be of the form 



'(sy^eis+^dT'o- 



For t^ we shall have an exactly similar equation with the 
same coefficients. 

The above equation is, by the solution of a quadratic, 
resolvable into two equations of the form 

du ^ du ^ du ^ du ^ 

To these correspond the respective auxiliary equations 

dt+\ds = Oy dt + \ds^O (13). 

If the integrals of these are 
respectively, then we have 

Now V beinff determinable by an equation of the same 
form as u, it follows that of the above two values of u one 
must be assigned to v, so that the solution of the problem will 
be contained in the system 

or In the system 
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The particular forms of the arbitrary fdnctions ^ afe 
will depend solely upon the nature of the problem under < 
sideration. 

One other point remains to be noticed. The first mem- 
bers of (12) are essentially positive, being composed of 
squares ; so are then the first members of (I) , (II) ; and so, 
if the intermediate transformations are real, is the first 
member of the eg^uation whose coefficients are P, Q, R, 
Hence the quadratic determining X^, \ will have imaginary 
roots of the form a±pj—l. Ultimately therefore it will suf- 
fice to integrate one equation of the system (13) and then to 
deduce the solution of the other by changing nf^ into 

3. Application of the above formuke when the given sur- 
face is an oblate spheroid, such as the earth. 

Let the plane of the equator be that of projection, the 
centre being the origin. Let the co-ordinates x, y pass 
through the meridians of and of 90** respectively, and z 
through the poles. The equation of the surface will be 

a^ ^b^ ^^*^' 

where a is the earth's equatorial, b its polar radius. Let also 
the latitude of the point a?, y, z be represented by s, the 
longitude by t. We have 

^ a^ ^b' ^' 

dF_2x dF_2y dF _2z 

cfo;~Z' dy" a^' dz'^W' 

and substituting in (I), 

fx du y^ du z duV _ 
~V dbs^a^l^'^b^'dz) ~ ' 
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xi we represent ^ by A*, 

(...,..*.^{(|)V(|)V(g)] 

Now as Xf y are rectangalar co-ordinates in tibe plane of 
the equator^ and x passes through the first meridian, we have 

^ = tan ^. 

X 

Again, representing in the annexed figure the meridian of 
the point P^ or {x, y, z) 
touched by the straight line 
QB in the same plane, we 
have g3f=V ^+y, MP=-z. 
(Therefore if *Jx' + y^=r, the 
equation of the meridian is 

--4---=:1 




that of the tangent 



rr zz _ 



r , z being current rectangular co-ordinates of the tangent. 
Hence 

^ ^-^ c?z Wz 

Vr V5n^ 



But CQR = latitude. Therefore finally 



8 = tan"' 



Vx'+y*' 



X 



(16), 



and we must now transform (15) so as to make s and t the 
independent variables. 
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From the above equations combined with (14) we find ^ 

_ aAcos^ _ aAsin^ _ atan^ , . 

^"^h' + Un'? ^"VA" + tan«i' ^~AVF+^^'"^ ^' 

and substituting in (15), 

-(co8*^ + 8m«^ + tan«^j=0 (18). 



sec 



Again, 



du ^du da du dt 
dx^ da dx dt dx^ 



du 
dy' 


du da 
'da dy^ 


du dt 
dt dy' 




du 
da' 


^du da du dt 
"^ da dzlidz' 




- 


k'zx 


- 


sin« cos« costJS 



Now 3- = -p=== =- 

^ ^/^4y(i»"+y + A*«') «^ 

where jEr= A' + tan'«. In like manner 

cfe _ — sing cosg sintt/S 
dy'^ ah ' 

da _ h cos'g ^.g 
cfe"" a ' 

dt —BJrLtJS , 
dx'~ ah ' 

dt co^tjE 
ay' ah ' 



*«0 
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Hence 

du JH ( . du . du\ 

Substituting these values in (18), and dividing by the com- 
mon factor -jYT we have on deduction 

(f y+ C08«* {1 + (A'- 1) C08M'@'=0, 

which is resolvable into 

^^ - 7=1 C08« {1 + (*•- 1) cos^} ^ =0, 

2^ + Npl C0S5 {1 + (A*- 1) C08*«}^ « 0, 

partial differential equations of which the integrals are in- 
cluded in the comnion formula 

^*"' ic08»{l+(A'-l)c08*«} - 
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« t^l' , . a I (since e'=s — 5— 1 

« log tan (-7 + 5 + ^ log -. — ► 

Hence 

«-*h{Cri^)''-(M)}*'^]' 

or, changing ^ (f) into ^ (e*), 

Let r and d be the polar co-ordinates of that point in the 

Elane of projection which corresponds to the point whose 
ititade and longitude on the surface are $ and t ; and let 

/l-esin»\l^/2r + i), 
\l + eam$J U 2/* 

then the complete solution assumes the very simple form 

r6«vir= <l> (fife*«>/=Oi re-«v=J= -^ (56*«v=0 (HI). 

Of particular deductions the most interesting is that which 
arises from the supposition that the parallels of latitude are 
projected into circles round the pole. This requires that r 

B.D.E. lu 16 
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should be independent of .<; a condition which is satisfied in 
the most g^nei:^! manner by assuming . 

we then find 

whence, on multiplication and division, 

u 

whence, A and S being new arbitrary constants derived from 
Oand C • *' 

If we observe that and t should vanish together, we have 
5 = 0, and tjie equation 0^ ±nt shews that the surface of 
the sphere will be projected into a sector of a circle, the arc 
of which is to the circumference of the circle as n : 1. Thus, 

if n = -, thQ sphere is projected upon a quadrant^ and so 



on. 



• The other equation gives 

\ V* 2/J \1 + e sin 8/ 

li 8 = we find r = A, whence A is the distance of the 
equator from the pole in the plane of projection, and if that 
distance, which is arbitrary, be assumed as the unit, we have' 



-tan{(^+|U7;'-^^!"^)^ 
[Vl 2/j XX+esmsJ. 



for the distance from the pole of that parallel whose latitude 
is s. We may give to this expression a better form by 
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assuming jp^-^ +8, and introducing an auxiliary quantity q 
determined by the equation 

ecosp = cosj. 
We have then 

r.(unf)-(oo.|)-. 

The following table gives the values of r for the sphere and 
for the spheroid whose eccentricity is '08 (which is about that 
of the earth), for each ten degrees of polar distance, for the 

v-alues n = 1, and ti = -r • 
4 



Polar 


n= 


= 1 


n 


1 

■"4 


Distance. 










Sphere. 


Spheroid. 


Sphere. 


Spheroid. 


10« 


•0875 


•0880 


•5439 


-5447 


200 


•1763 


•1774 


•6480 


•6490 


30° 


•2679 


•2694 


•7195 


•7205 


40' 


•3640 


•3658 


•7767 


•7777 


SO^' 


•4663 


•4682 


•8264 


•8272 


60° 


•6774 


•5792 


•8717 


•8724 


W 


•7002 


•7017 


•9148 


•9153 


80» 


•8391 


•8400 


•9571 


•9574 


W 


1-0000 


1-0000 


1-0000 


1-0000 


100® 


1-1918 


1-1904 


1-0448 


10445 


\W 


1-4281 


1-4250 


1-0932 


1-0926 


120<> 


1-7321 


1-7265 


1-1472 


11463 


130^ 


21445 


2-1357 


1-2101 


1-2089 


140<> 


2-7475 


2-7340 


1-2875 


1-2859 


1500 


3-7321 


3-7114 


1-3899 


1-3880 


1600 


6-6713 


5-6372 


1-5432 


1-5409 


1700 


11-4301 


11-3581 


1-8387 


1-8358 



THE END. 
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